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INVESTIGATION  OF  THE  LOW  SPEED 


FIXED  GEOMETRY  SCRAMJET 

"INLET  DESIGN  PRACTICE  MANUAL" 


Edited  by  James  Johnson 


FOREWORD 


This  report  was  prepared  by  the  General  Applied  Science 
Laboratories,  Inc-,  Westbury,  New  York,  a  wholly  owned  sub¬ 
sidiary  of  The  Marquardt  Corporation,  on  Contract  F-33(615)- 
67-C-1084,  Task  3012,  "Investigation  of  the  Low  Speed  Fixed 
Geometry  Scramjet-"  The  work  was  administered  under  the 
direction  of  the  Air  Force  Aero  Propulsion  Laboratory,  Air 
Force  Systems  Command.  Mr.  R.  Canny,  (APRP)  was  the  project 
engineer  for  the  laboratory. 

The  studies  presented  began  in  September  1966  and  were 
concluded  in  August  1967.  Mr.  James  Johnson  was  responsible 
for  the  technical  direction  of  this  program  at  GASL.  The 
technical  efforts  were  performed  by  the  following  scientists; 
Messrs.  M.  Abbett  and  G.  Bleich  for  the  two-dimensional  flow 
analysis  and  programming,  and  Drs.  S.  Angelucci  and  C.  Ruger 
for  the  thr je-dimensional  flow  analysis  and  programming. 

This  report,  GASL  TR-667,  was  submitted  November  1967. 


ABSTRACT 


Computer  progr-:-jns  have  been  written  to  aid  the  design  and 
analysis  of  two-dimensional,  axisymmetric,  and  three-dimensional 
supersonic  inlets.  This  report  presents  the  fundamental  analytic 
techniques,  the  use  and  operation  of  the  computer  programs,  and 
the  application  to  the  design  of  supersonic  inlets.  The  programs 
are  written  in  Fortran  IV  for  use  on  the  7094  high  speed  digital 
computer . 

The  two-dimensional  and  axisymmetric  programs  presented 
herein  are  written  for  generalized  inviscid  supersonic  internal 
flow  problems  with  uniform  or  non-uniform  inlet  entry  conditions 
for  entropy,  total  enthalpy,  pressure,  Mach  number  and  flow 
direction.  The  program  capabilities  include  the  intersections 
and  reflections  of  both  family  waves,  the  formation  of  shocks  or 
expansions  at  corners,  the  formation  of  shocks  by  coalescence  of 
waves  from  smooth  walls,  and  the  formation  of  contact  discon¬ 
tinuities  . 

The  three-dimensional  programs  presented  herein  are  written 
for  the  calculation  of  the  inviscid  supersonic  flow  fields  associ¬ 
ated  with  basic  elements  of  three-dimensional  supersonic  inlets. 

The  methods  utilized  exact  and  linearized  supersonic  flow  theory 
with  engineering  approximations  to  yield  solutions  to  the  following 
unit  problems;  delta  wing  flow,  conical  shock  interacting  with  a 
plane  surface,  plane  shock  interacting  with  a  conical  surface,  and 
the  interaction  of  two  different  conical  flow  fields. 
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NOMENCLATURE  FOR  SECTION  IV 
speed  of  sound,  a8  =  y  P/p 
J  Ma-1  /  p  u3 

coefficients  of  specific  heat 


specific  static  enthalpy, 


specific  total  enthalpy,  H  =  h  +  q3/2 


Mach  number ,  M  =  q/a 


pressure 

magnitude  of  the  velocity  vector ,  q8  =  u3  +  v8 

gas  constant 

entropy 

velocity  component  in  x  direction 
velocity  component  in  y  direction 

abscissa  (horizontal  axis  =  z  in  cylindrical  polar 
frame  (z,  r,  9)  for  axisymmetric  flow) 

ordinate  (vertical  axis  =  4  in  cylindrical  polar 
frame  for  axisymmetric  flow) 


angle  flow  turns  when  going  through  oblique  shock 


c  /c  ,  the  ratio  of  specific  heats 
p  v 

characteristic  slope 


viscosity  (=0  in  this  report 


V 


Prandtl-Meyer  turning  angle  (see  equation  15) 


p  density 

or  angle  oblique  shock  makes  with  incident  flow 

r  streamline  slope,  r  =  v/u 

A  angle  velocity  vector  makes  with  x  axis 

Subscripts 


1,  2  conditions  in  front  of  and  behind  shock,  respectively 

A,B,C,etc.  variable  is  to  be  evaluated  at  the  indicated  point 
Superscripts 


j  0  for  two-dimensional  flow 

1  for  ax i symmetric  flow 

- -  averaged  quantity,  also  indicates  a  line  between  two 

points 

I,  II  down  running  and  up  running  characteristic' 

respectively 

The  following  is  standard  notation  on  the  figures 

_  down  running  characteristic 

-  up  running  characteristic 


oblique  shock 

contact  discontinuity  (interface,  slip  stream) 
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INTRODUCTION 


The  use  of  three-dimensional  flow  fields  in  the  design  of 
supersonic  and  hypersonic  inlets  has  been  under  experimental 
and  analytical  investigation  at  GASL  for  several  years  under 
Air  Force  sponsorship  as  reported  in  References  1  through  6. 
Three-dimensional  designs  can  be  used  for  either  fixed  or  vari¬ 
able  geometry  inlet  configurations  over  a  wide  flight  Mach  number 
range  to  obtain,  for  example,  high  inlet  contraction  ratios  and 
capture  mass  flows,  to  prevent  boundary  layer  separation,  and  to 
reduce  heating  loads  to  the  inlet  structure.  General  use  of 
three-dimensional  designs  however,  has  been  limited  by  the  devel¬ 
opment  of  practical  analytic  techniques.  In  most  cases  where 
three-dimensional  analytic  techniques  are  available,  they  are 
either  so  highly  restricted,  or  complicated  and  time  consuming, 
that  they  are  impractical  for  use  in  a  typical  inlet  design  effort. 

Practical  methods  the  design  and  analysis  of  three- 
dimensional  supersonic  inlets  have  been  under  development  at 
GASL  for  the  Air  Force  Aero  Propulsion  Laboratory  as  part  of 
the  overall  investigations  of  the  low  speed  range  supersonic 
combustion  ramjet  (SCRAMJET) .  These  methods  combine  exact 
and  linearized  three-dimensional  theory  with  engineering  ap¬ 
proximations  to  yield  solutions  to  unit  problems  which  are  basic 
elements  of  three-dimensional  supersonic  internal  flow  fields, 
i.e.,  delta  wing  flow,  coniual  shock  interacting  with  a  plane 
surface,  plane  shock  interacting  with  a  conical  surface,  and 
the  interaction  of  two  different  conical  flow  fields.  A  com¬ 
bination  of  these  unit  problems  can  be  used  to  design  and 
analyze  a  large  variety  of  supersonic  and  hypersonic  three- 
dimensional  inlets. 

These  methods  for  two  and  three  dimensional  inlet  design 
and  analysis, have  been  programmed  in  Fortran  IV  for  use  on 
the  IBM  7094  digital  computer.  This  report  presents  the  com¬ 
plete  description  of  the  fundamental  analysis,  the  program  use 
and  operation,  and  the  application  to  the  design  of  three  and 
two  dimensional  supersonic  inlets.  This  work  was  performed 
for  the  Wright-Patterson  Air  Force  Base  Aero  Propulsion  Labora¬ 
tory  under  contract  F33  (615) -67-C-1084. 


II 


PROGRAM  OBJECTIVES  AND  APPROACH 


The  overall  objective  of  this  effort  was  to  generate  digital 
computer  programs  for  aid  in  the  design  and  analysis  of  two  and 
three  dimensional  supersonic  inlets.  The  analytic  techniques 
which  form  the  basis  of  this  work  have  been  developed  at  GASL 
during  the  past  several,  years  under  Air  Force  sponsorship.  The 
purpose  of  this  report  is  to  present  these  fundamental  analytic 
techniques,  the  use  and  operation  of  the  computer  programs,  and 
the  application  to  the  design  of  supersonic  inlets.  The  aim  is 
to  provide  a  practical  tool  for  two  and  three  dimensional  inlet 
design,  in  addition  to  a  detailed  description  of  the  analysis. 

To  this  end,  the  following  sections  describe; 

Section  III  -  A  discussion  of  the  methodology  used  for  three- 
dimensional  inl^t  design  and  the  use  of  the 
computer  programs  in  these  efforts. 

Section  IV  -  A  complete  description  of  the  two-dimensional 
and  axisymmetric  internal  flow  field  analysis 
and  computer  programs. 

Section  V  -  A  complete  description  of  the  three-dimensional 
analysis  and  computer  program  elements. 

It  is  seen  then,  that  this  report  can  be  used  as  a  "Design 
Manual"  to  aid  in  three-dimensional  inlet  design,  or  as  an  initial 
step  in  developing  a  more  sophisticated  or  modified  analysis.  In 
either  case,  it  must  be  noted  that  the  current  work  presented  here 
has  limitations  which  should  be  the  object  of  continuing  study. 
These  limitations  will  be  clear  as  they  arise  in  the  following 
sections  in  this  report.  A  summary  of  the  overall  limitations 
however,  is  given  here. 

(1)  Viscous  effects  are  neglected 

(2)  Completely  supersonic  flow  required 

(3)  Completely  supersonic  leading  edgss  required 

(4)  Three-dimensional  computations  only  for  unit  problems. 
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METHODS  FOR  THREE-DIMENSIONAL  INLET  DESIGN 


The  detailed  three-dimensional  and  two-dimensional  analytic 
techniques  and  computer  programs  will  be  described  in  the  follow¬ 
ing  sections.  It  is  necessary  here,  however,  to  describe  the  over¬ 
all  methodology  for  using  these  techniques  before  going  into  the 
analysis. 

The  overall  inlet  design  process  is  schematically  shown 
in  Figure  1.  It  is  seen  that  several  steps  are  required  in 
the  design  process: 

a)  choice  or  definition  of  inlet  design  performance 
characteristics 

b)  choice  or  definition  of  inlet  design  concept,  e.g., 
two-dimensional,  three-dimensional,  with  or  without 
centerbody ,  etc . 

c)  two-dimensional  calculations  to  define  essential  per¬ 
formance  properties 

d)  three-dimensional  corrections  to  define  exact  perfor¬ 
mance  propetties 

e)  evaluation  of  performance  with  respect  to  initial  re¬ 
quired  performance  levels 

f)  repeat  steps  (b)  to  (e)  until  required  performance  is 
obtained . 

A  summary  of  each  of  the  above  design  process  steps  is  given 
in  the  following  discussion. 

The  initial  steps  (a)  and  (b)  are  closely  interrelated  since 
the  required  inlet  design  performance  characteristics  generally 
dictate  the  choice  of  the  inlet  design  concept.  The  requirement 
of  fixed  geometry  operation  over  a  wide  Mach  number  range,  for 
example,  might  require  a  more  highly  three-dimensional  inlet 
design  (with  possibly  one  or  more  inlet  central  bodies)  than  a 
fixed  geometry  inlet  design  for  a  more  restricted  Mach  number 
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FIGURE  1  -  INLET  DESIGN  PROCESS 
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variable  geometry  inlet  design  to  attain  performance  goals.  In 
any  case,  these  initial  steps  of  the  design  process  require 
decisions  based  upon  inlet  design  experience. 

The  next  step  (c)  utilizes  the  two-dimensional  and/or  axi- 
symmetric  computer  programs  described  in  Section  IV.  This  step 
would  be  the  final  step  for  a  completely  two-dimensional  or 
axisymmetric  design.  For  the  case  of  a  three-dimensional  design, 
however,  it  forms  the  first  step.  Several  overall  inlet  design 
geometries  are  roughed  out  by  the  inlet  designer  and  the  computer 
program  provides  a  detailed  analysis  of  each  one.  The  overall 
properties  of  the  various  inlet  flow  fields  are  then  evaluated 
with  respect  to  the  initial  design  goals  and  the  most  promising 
designs  are  chosen  for  further  consideration  by  including  the 
effects  of  three-dimensionality. 

The  next  step  (d)  utilizes  the  three-dimensional  computer 
program  elements  described  in  Section  V.  The  basic  elements  of 
a  large  variety  of  supersonic  three-dimensional  inlet  flow  fields 
are  programmed  for  machine  calculation  in  this  phase  of  the  design 
process.  These  elements  as  shown  in  Figure  2  are  the  delta  wing 
flow  field,  the  interaction  of  a  conical  shock  with  a  plane  sur¬ 
face,  the  interaction  of  a  plane  shock  with  a  conical  surface,,  and 
the  interaction  of  two  different  conical  shocks.  The  techniques 
used  for  these  machine  computations  however,  can  also  be  used 
manually  to  compute  the  flow  field  of  almost  any  supersonic  three- 
dimensional  configuration.  The  level  of  three  dimensionality 
designed  into  the  inlet  geometry  in  this  step  of  course  depends 
again  upon  the  initial  design  goals.  The  overall  procedure  here 
again  is  to  investigate  various  geometries  until  the  desired 
performance  conditions  are  achieved.  However,  here  the  overall 
configuration  is  known  from  the  two-dimensional  and  axi symmetric 
analysis  of  step  (c) .  The  procedures  in  this  step  are  to  incor¬ 
porate,  for  example,  various  leading  edge  sweep  angles,  cone 
angles,  etc. 

An  inlet  configuration  is  now  defined  and  analyzed.  Step  (e) 
calls  for  an  evaluation  of  the  resulting  performance  based  upon 
the  initial  design  performance  goals.  At  this  point,  either  these 
goals  are  satisfied  or  steps  (c)  and  (d)  are  repeated.  In  this 
manner  then  the  required  inlet  performance  characteristics  are 
achieved . 


> 


IV 


TWO-DIMENSIONAL  AND  AXISYMMETRIC  SUPERSONIC  INTERNAL 
FLOW  COMPUTER  PROGRAM 


A.  Summary  and  Introduction 

A  computer  program  has  been  written  for  the  computation 
of  two-dimensional  and  axisymmetric  inviscid  supersonic  internal 
flows  for  configurations  typical  of  supersonic  and  hypersonic 
inlets.  The  program  is  designed  for  use  either  for  the  design 
and  analysis  of  completely  two-dimensional  or  axisymmetric 
inlets,  or  as  part  of  the  design  effort  of  a  three-dimensional 
inlet.  as  described  in  Section  III.  The  program  is  written  in 
Fortran  IV  for  use  on  the  IBM  7094  and  the  CDC  6600  high  speed 
digital  computers. 

The  overall  computer  program  can  be  summarized  as 

follows: 


General  -  Calculation  of  complete  inviscid  two-dimensional 
and  axisymmetric  supersonic  internal  flow  fields. 

Input  -  Uniform  or  non-uniform  entry  conditions  for 

entropy,  total  enthalpy,  pressure,  Mach  number 
and  flow  direction. 

Capabilities 


•  Intersection  and  reflections  of  waves  of  same 
and  opposite  families  -  no  limit  to  number. 

•  Formation  of  shocks  at  compression  corners. 

•  For-  ation  of  shocks  by  coalesence  of  waves 
from  smooth  walls. 

•  Formation  of  contact  discontinuities. 

•  Reflection  of  waves  from  walls. 

•  Perfect  gas. 
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Output  -  Fluid  mechanic  and  thermodynamic  parameters,  wave 
locations  and  contact  discontinuity  location  at 
each  point  of  the  characteristic  mesh. 

A  complete  detailed  description  of  this  computer  program  is 
given  in  this  section.  Part  B  describes  the  types  of  flews  con¬ 
sidered  and  summarizes  the  assumptions  and  simplifications  in  the 
analysis  and  computer  program.  In  Part  C  the  Euler  equations  of 
motion  are  presented  and  non-dimensionalized .  In  addition,  the 
special  forms  of  the  governing  equations  are  given  here  for  charac¬ 
teristic,  oblique  shock,  and  Prandtl  Meyer  expansion  computations. 

Part  D  describes  the  manner  in  which  the  computation  is  performed 
and  gives  outlines  of  the  various  types  of  computations.  Appendix 
I  gives  most  of  the  information  required  by  the  user  of  the  program. 

In  Appendix  III  some  sample  calculations  are  presented. 

A  person  who  wishes  to  use  the  program  with  a  minimum 
amount  of  effort  should  skim  Part  B  and  Part  D.  He  should  care¬ 
fully  read,  in  order.  Part  C-a  and  Appendices  I-III . 

B.  Outline  of  the  Program 

a)  Type  of  Flows  Considered 

The  program  computes  internal  flows  typical  of 
supersonic  and  hypersonic  inlets.  Included  in  the  flow  field  com¬ 
putation  are  weak  shocks  of  both  families,  their  formation  at  sharp 
corners  and  by  the  coalesence  of  characteristics  (envelope  shocks) , 
the  intersection  of  shocks  of  the  same  and  opposite  families,  the  reflec¬ 
tion  of  shocks  on  the  inlet  walls,  the  intersection  of  weak  shocks  and 
slip  streams,  and  centered  corner  expansions.  The  gas  is  considered 
to  be  an  inviscid  perfect  gas  of  constant  ratio  of  specific  heats. 

Since  the  method  of  characteristics  is  used  to  perform  the  calculation, 
the  flow  must  remain  supersonic  throughout.  Strong  shocks  are  not 
considered . 

b)  General  Comments  -  Approximations  and  Simplifications 

A  schematic  representation  of  a  typical  inlet  is  given 
in  Figure  3.  The  flow  impinges  on  the  inlet  at  some  specified  angle 
of  aetack.  There  may  be  an  expansion,  a  shock,  or  neither,  at  the 
leading  edge  of  the  lower  wall  and  at  the  cowl.  Due  to  smooth  or 
abrupt  compressions  along  the  two  walls,  shocks  may  form.  In  addition, 
expansion  waves  may  be  generated  by  smooth  or  sharp  expansions  of 
the  walls.  The  shocks  and  expansions  will  then  interact  in  complicated 
ways. 
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FIGURE  3  -  TYPICAL  INLET  CONFIGURATION  (AXI SYMMETRIC) 


Because  of  storage  and  time  considerations  and  the 
desire  to  keep  the  program  as  simple  as  possible  while  retaining 
sufficient  accuracy,  various  assumptions  and  simplifications  have 
been  made.  These  are  outlined  here  with  their  justification. 

1.  Intersection  of  Two  Weak  Shocks  of  the 

Same  Family 

In  general,  such  an  intersection  will  result 
in  the  two  shocks  continuing  as  either  two  shocks  and  a  contact 
di  -continuity  or  as  one  shock,  an  expansion  fan,  and  a  contact 
discontinuity  (Figure  4).  However,  in  cases  of  interest  here 
the  second  shock  or  expansion  fan  (OS  of  Figure  4)  and  the  con¬ 
tact  disconuity  (OD)  usually  have  a  much  smaller  effect  on  the 
flow  field  than  does  the  primary  shock  (OC) .  Hence,  in  the 
program,  two  shocks  of  the  same  family  intersecting  are  continued 
as  only  one  shock  of  that  family. 

2.  Intersection  of  Two  Shocks  of  Opposite  Families 

When  two  weak  shocks  of  opposite  families  inter¬ 
sect  they  will  in  general  continue  as  rwo  shocks  of  opposite  families 
and  a  contact  discontinuity  (Figure  5) .  In  the  program  such  an  inter¬ 
section  is  treated  with  no  simplifications. 

3.  Intersection  of  a  Weak  Shock  and  a  Contact 

Discontinuity 

Ordinarily  such  an  intersection  will  continue 
as  a  shock,  a  centered  expansion,  and  a  contact  discontinuity 
(Figure  6) .  Fresently,  the  contact  discontinuity  is  dropped  when 
it  intersects  an  oblique  shock. 

4.  Formation  of  Attached  Oblique  Shocks 

By  Discontinuous  Compression 

This  is  done  in  a  straightforward  fashion 
(Figure  7)  with  no  simplifying  assumptions. 

5.  Reflection  of  Oblique  Shocks  from  Walls 

This  is  done  in  a  staight forward  fashion,  the 
shock  being  reflected  with  strength  such  that  the  reflected  shock 
has  a  turning  angle  which  is  minus  that  of  the  incident  shock. 

(Figure  8)  . 
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FIGURE  4  -  INTERSECTION  OF. TWO  SHOCKS  OF  THE  SAME  FAMILY 


FIGURE  5  -  INTERSECTION  OF  TWO  SHOCKS  OF  OPPOSITE  FAMILIES 


Expansion 


FIGURE  6  -  INTERSECTION  OF  AN  OBLIQUE  SHOCK 
AND  A  CONTACT  DISCONTINUITY 


Shock 


FIGURE  8  -  REFLECTION  OF  AN  OBLIQUE  SHOCK  FROM  A  WALL 
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6.  Formation  of  Detached  Oblique  Shocks  1 

By  Coalescence  of  Characteristics  j 

\ 

After  each  point  of  the  characteristic  mesh  i 

has  been  computed,  tests  are  made  to  detect  the  crossing  of  character-  '  | 

istics  of  the  same  family,  indicating  the  formation  of  an  oblique 
shock.  When  such  a  crossing  occurs,  a  shock  is  entered  in  the  flow 
field  at  that  point  and  its  subsequent  development  is  followed. 

(Figure  9) 

j 

7.  Centered  Expansions  Due  to  Discontinuous 

Expansion  Corners  5 

These  are  performed  by  dividing  the  turning  angle 
of  the  corner  into  a  finite  number  of  segments  and  then  making  ] 

regular  characteristic  computations.  (Figure  10) 

8.  Reflection  of  an  Oblique  Shock  from  the  Axis 

of  Symmetry 

In  the  case  of  axisymmetric  flow,  the  possible 
flow  configurations  arising  when  an  oblique  shock  reflects  on  the 
axis  of  symmetry,  are  several  and  complicated.  There  is  no  provision 
in  the  program  to  treat  this  problem.  One  can  expect  the  program 
to  usually  blow  up  as  an  oblique  shock  approaches  the  symmetry  axis. 

However,  it  may  be  that  the  shock  will  reflect,  as  from  a  wall,  away 
from  the  symmetry  axis. 

9.  Compression  Corner  on  the  Symmetry  Axis  in 

Axisymmetric  Flow 

I 

( 

The  only  compression  corner  allowed  on  the 
symmetry  axis  in  axisymmetric  flow  is  at  the  leading  edge  of  the  , 

lower  wall.  In  addition,  the  initial  data  line  must  go  through  > 

the  leading  edge  of  the  lowe.  wall.  (See  Part  D-a)  ' 


FIGURE  9  -  FORMATION  OF  A  DETACHED  OBLIQUE  SHOCK 
BY  COALESCENCE  OF  COMPRESSION  WAVES 


\ 

a 


FIGURE  10  -  CENTERED  EXPANSIONS  DUE  TO 

DISCONTINUOUS  EXPANSION  CORNERS 


where  x  and  y  are  axes  of  a  right  handed  cartesian  frame 
(cylindrical  frame)  with  velocity  components  u  and  v  re¬ 
spectively.  The  x  axis  is  positive  in  the  general  direction 
of  the  incident  flow.  j=0  for  two  dimensional  flow  and  j=l 
for  axi symmetrical  flow,  p  is  the  density,  p  the  pressure,  S 
the  entropy,  and  the  bars  indicate  dimensional  quantities. 
Now  we  divide  all  pressures  by  reference  pressure  p  , 

densities  by  p  ,  temperatures  by  T  ,  lengths  by  1  ,  veloc- 
.  .  ,  Q.  /=-= — i  -  o  ^  o 

lties  by  VPo/po'=  T0  where  R  is  the  gas 

constant  and  pQ  -=  pQ  RT0  ,  and  we  subtract  S,  from  all  en- 

J'.--ide  the  difference  by  c  . 

v 

u,  V,  x,  y, 

S  =  (S-C  /c 


tropie"  •!■>. 
dimension- 
ar  e ,  (w 


o 

the 
•-  itities 

'  P  = 


P,  D 
P/P  o 


V 


o 

Then,  in  non- 
S,  the  equations 
etc . 


(puyJ)  + 


x 


(pvy  )  = 

y 


o 


(4) 


puux  +  pvuy  +  Px  =  0 


puv  +  pvv  +  p  =  0 
x  y  y 


uS  +  vS  —  0 
x  y 


(5a) 

(5b) 

(b) 


where 


S  =  In  (p)  -  y  In  (p) 


(7) 


The  user  of  this  program  is  free  to  pick  any  reference 
length,  I  ,  he  chooses.  If  the  incident  flow  is  uniform  and  the 
simplified  input  is  used  (see  Appendix  I) ,  the  reference 
pressure  and  temperature  are  chosen  to  be  the  free  stream  (i.e. 
incident)  static  pressure  and  temperature.  The  non-dimensional 
specific  static  enthalpy  is  given  by 


h 


£ 

p 


and  the  specific  total  enthalpy  by 


H  =  h 


+ 


2  2 
u  +v 


2 


(8) 


(9) 
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7  2  2  2 

u  +  v  )  is  the  modulus  of  the  velocity  and  a  (=yp/p  ) 

is  the  square  of  the  speed  of  souuu. 

Since  the  flow  is  inviscid  and  steady,  H  and  S  are  con¬ 
stant  along  a  streamline.  When  there  are  shocks,  the  entropy,  s, 
suffers  a  jump  increase. 


For  completeness,  the  remaining  parts  of  this  section  v/ill 
present  well  known  special  forms  of  the  governing  equations  to  be 
used  when  we  make  certain  computations,  e.g.,  characteristic  compu¬ 
tations.  All  variables  in  the  equations  are  non-dimensional.  Those 
readers  wishing  more  detailed  explanations  of  the  mathematical 
character  of  supersonic, inviscid, perfect  gas  flows  should  look  in 
any  standard  text  book  on  gas  dynamics,  e.g.,  References  7,8,9. 

b)  Characteristic  Equations 


Along  characteristic,  or  Mach,  lines,  in  the  absence 
of  shock  waves,  the  Euler  equations  can  be  combined  so  that  we  have 
only  total  derivatives  of  certain  quantities  in  the  direction  of  the 
characteristic  line.  It  is  well  known  that  for  our  cases  the 
chair  act  eristic  directions  are  given  by 


/  2 

J  u  +  v 


.1  ch£  uv  -  a 

K  ~  dx  ~  2  2 

u  -  a 


2  2 
•  a 


(10a) 


.II  dy  uv 
X  =-£  =  ~~  2  2 


/  2 

+  a  J  u + 


2  2 
v  -  a 


(10b) 


u  -  a 


X111  =  =  Z  =  T 

A  dx  u 


(10c) 


Henceforth,  X1  will  be  referred  to  as  a  down-running  characteristic 

II  I  II 

and  X  as  an  up-running  characteristic.  Along  X  '  the  Euler 

equations  can  be  combined  to  give 


dr  +  b  = 

dx  ”  dx 


ll. 


X1'11)  t 


(ID 
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where 


and  M  (=q/a  )  is  the  Mach  number 


III 


b  =  J M2  -  1  /pu2 
Along  the  streamline,  we  have 

dS  =  dH  =  0  (no  snocks)  (12) 

We  can  write  equations  (11a)  in  differential  form  as: 


At  -  bAp  = 


(t  -  X1)  T 


along  a  down  running  characteristic  and 


At  +  bAp  = 


(t  ~  X11)  t 


j  Ax  =  G  Ax  *j  (13a) 


Ax  -j  =  G11  Ax  -j  (13b) 


along  an  up  running  characteristic. 

c)  Oblique  Shocks  -  Rankine  Hugoniot  Equations 


Let  subscript  1  refer  to  the  conditions  in  front  of 
(upstream)  the  shock  and  subscript  2  refer  to  the  conditions  behind 
(downstream)  the  shock  and  let  a  be  the  acute  angle  the  incident 
flow  makes  with  the  shock  wave. 
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A  is  the  turning  angle  of  the  flow.  The  Rankine  Hugoniot  equations 
are 


2  2 

P2  2y  Ml  sin  a  -  (y-1) 

P1  y  +  1 


(14a) 


p2  (y+D  p2  +  (y-1)  Pj_ 
pj_  ~  (y+i)  px  +  (y-i)  p2 


(14b) 


+  In 


2  (p2/ Pi  “  1) 

tan  6  =  tan  t  - - - 

(tan  a  +  P2/p1) 


M2  =  - 

2  .  2 


2  2 

(y-1)  M^  sin  a  +  2 

2  2 
/  C\  r  rx  w  •  ..  ** 


2  2  2 
sin  (cr-6)  [2y  M^  sin  a  -  (y-1) ] 


(14c) 


(14d) 


2y  m2  sin2  a-  (y-1) 
y  +  1 


2  2  T 

(y+1)  m.  sin  a 

-  y  In  - - - - —  (14e) 

f »/_  1  \  M  ci'n  ‘  rr  +  2 


H2  H1  (14f ) 

tan  1  (t2)  =  tan-1  (t^)  +  6  (14g) 

(14h) 

d)  Prandtl-Meyer  Expansion 

The  Prandtl-Meyer  relation  is 


v 


(15) 


where  v  is  the  angle  turned  from  Mach  =  1  to  reach  the  Mach  number 
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M.  If  we  have  a  corner  with  other  than  sonic  velocity  just  prior 
to  the  corner,  then  denoting  the  conditions  just  prior  to  the  ex¬ 
pansion  by  subscript  1  and  just  after  the  expansion  by  subscript  2, 
we  have 


1  + 


1  + 


2 


y-1 


V 

Y-l 


,  y-l  2 
1  +  — Mi 


<y—  1  2 

1  +  VM2 


1 

nr-1 


y~l 


(16b) 


(16c) 


H 


2 


(16d) 


(16e) 

(16f) 


s2  =  sL 


t  =  tan  S 


(16g) 


D.  The  Program  -  General  Outline 

a)  Initial  Conditions  -  Boundary  Conditions  - 

Coordinate  Systems 

As  indicated  previously,  the  purpose  of  the  program 
is  to  compute  supersonic  and  hypersonic  inviscid  inlet  flows.  The 
coordinate  system  and  body  geometry  are  represented  schematically 
in  Figure  11.  The  x  axis  is  aligned  in  the  general  direction  of 
the  free  stream  flow  (it  is  the  axis  of  symmetry  in  the  case  of  axi- 
symmetric  flow)  but  need  not  be  exactly  parallel  to  the  incident 
flow.  It  is  positive  in  the  general  direction  of  the  incident  flow. 
The  y  axis  is  normal  to  the  x  axis  and  is  positive  in  the  direction 
determined  by  rotating  the  positive  x  axis  90°  counterclockwise.  It 
corresponds  to  the  radial  axis  in  a  cylindrical  (x,r,  0)  frame  in 
the  case  of  axisymmetric  flow.  It  is  assumed  that  the  tips  of  both 
the  cowl  and  the  lower  wall  are  pointed,  and  of  small  enough  included 
angle  relative  to  the  incident  flow  so  that  any  tip  shocks  are 
attached  and  there  are  no  regions  of  subsonic  flow.  Furthermore, 
for  axisymmetric  flow  with  a  lower  leading  edge  on  the  axis  of  sym¬ 
metry,  it  is  assumed  that,  in  the  coordinate  system  used,  the  x 
coordinate  of  the  tip  of  the  cowl  is  greater  than  or  equal  to  the 
x  coordinate  of  the  tip  of  the  lower  wall. 

The  line 'parallel' to  the  y  axis  through  the  tip  of 
the  lower  wall,  AA' ,  is  considered  to  be  the  "data  line"  of  the 
incident  flow  (see  Appendix  I  ). 

The  geometry  of  the  upper  and  lower  walls  must  be 

specified . 

b)  Order  of  Progression  of  The  Computation 

1.  In  The  Absence  of  Shocks 

The  computation  proceeds  along  down  running 
characteristics  issuing  either  from  points  on  the  initial  data 
line,  (AA 1 ,  Figure  11)  or  from  the  upper  wall.  Referring  to  Figure 
12,  let  us  follow  the  outline  of  a  typical  computation,  forgetting 
for  the  moment  the  possibility  of  shocks  or  centered  expansions. 

Mesh  points  are  indicated  by  circles  and  the  flow  is  completely 
specified  initially  along  the  initial  data  line,  AA'.  The  flow 
field  at  point  C  can  be  computed  for  the  known  conditions  at 
points  A  and  B  (see  Part  D-c-1) ,  where  C  is  the  point  of  inter¬ 
section  of  the  up  running  char acteri ^  .ic  from  A  and  the  down 
running  characteristic  from  B.  The  down  running  characteristic 
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FIGURE  11  -  SCHEMATIC  OF  COORDINATE  SYSTEM  AND  BODY  GEOMETRY 


from  C  intersects  the  wall  at  D,  and  the  flow  conditions  at  D 
can  then  be  obtained  from  the  known  conditions  at  C  and  the 
specified  wall  slope  at  D.  (See  D-c-2) 

The  computation  proceeds  in  this  manner 
working  along  down  running  characteristics  from  the  initial 
data  line  to  the  bottom  wall.  Eventually,  however,  one  of 
the  down  running  characteristics  will  cross  the  cowl  as  illus¬ 
trated  in  Figure  13a.  C,  the  point  formed  by  the  intersection 
of  the  characteristics  from  D  and  B,  is  computed  as  an  ordinary 
point.  Then  the  conditions  at  E,  the  tip  of  the  cowl,  are 
obtained  by  interpolation  between  points  A,  B,  C,  and  D.  The 
computation  is  continued  by  working  along  the  down  running 
characteristic  issuing  from  E  (forget  for  the  moment  the 
possibility  of  a  shock  at  E) ,  Figure  13b.  When  that  character¬ 
istic  is  completed,  the  intersection,  with  the  upper  wall  of  the 
up  running  characteristic  FG,  from  point  F,  can  be  determined. 

The  flow  conditions  at  G  can  then  be  computed  using  the  known 
conditions  at  F  and  the  specified  wall  slope  at  G.  Then  the 
computation  can  be  continued  in  a  normal  fashion  working  along 
the  down  running  characteristic  issuing  from  G. 

In  the  absence  of  shocks  and  centered 
corner  expansions,  the  entire  inlet  flow  field  could  be  com¬ 
puted  in  this  fashion.  In  actual  fact  however,  there  will 
often  be  not  only  one,  but  several  oblique  shocks  (of  either 
or  both  families)  which  complicate  matters  somewhat. 

2.  In  the  Presence  of  Attached  Shocks 

Although  up  and  down  running  characteristics 
have  equal  stature  in  the  mathematical  theory  of  supersonic  in- 
viscid  flow,  the  situation  is  somewhat  different  in  the  program. 
That  is  to  say,  since  the  computation  proceeds  along  down  running 
characteristics,  we  have  implicitly  created  a  favoritism  between 
the  characteristics  of  opposite  families.  This  difference  is 
important  when  shocks  and  contact  discontinuities  are  present. 

This  unequal  treatment  of  characteristics 
of  different  families  does  not  affect  the  basic  order  of  the 
computation  when  up  running  shocks  are  present.  (See  Figure  14) 
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FIGURE  13  -  MODIFICATION  OF  CHARACTERISTIC  MESH  AT  THE  COWL 


FIGURE  14  -  SKETCH  OF  DOWN  RUNNING  AND  UP  RUNNING  SHOCKS 


FIGURE  15  -  PATTERN  OF  DOWN  RUNNING  CHARACTERISTICS 

NEAR  DOWN  RUNNING  SHOCK 
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If  there  are  such  shocks,  the  intersection  of  each  shock  with 
each  down  running  characteristic  defines  a  mesh  point  (actually 
two,  one  in  front  of  and  one  behind  the  shock,  see  Part  D-e-1) . 

The  only  complication  in  trie  program  is  that,  as  it  proceeds 
along  a  down  running  characteristic,  the  program  must  be  instruc¬ 
ted  to  look  for  the  intersection  of  the  down  running  character¬ 
istic  with  any  up  running  shocks  from  the  previous  down  running 
characteristics.  When  such  an  intersection  occurs,  a  special  com¬ 
putation  must  be  made. 

However,  the  situation  is  not  quite  so  simple 
when  there  are  down  running  shocks.  For,  the  down  running  char¬ 
acteristics  down  stream  of  the  shock  terminate  at  the  shock  (see 
Figure  15),  i.e.  characteristics  terminate  at  a  shock  of  the  same 
family.  Thus,  as  the  program  proceeds  along  a  down  running 
characteristic,  it  must  make  a  special  provision  upon  encountering 
a  down  running  shock  since  the  down  running  characteristic  does 
not  continue  on  the  other  side  (up  stream)  of  the  shock. 

The  proceedure  that  is  followed  in  such  cases 
is  illustrated  in  Figure  16  for  the  case  of  a  compression  corner 
on  the  upper  wall.  When  such  a  corner  is  detected,  the  con¬ 
ditions  upstream  of  the  shock  are  obtained  by  a  regular  character¬ 
istic  wall  computation  between  the  corner  and  A,  the  data  at  A 
being  computed  by  interpolation  on  the  previous  down  running 
characteristic.  The  program  then  computes  the  down  running 
characteristic,  BC,  that  would  issue  from  B  if  no  compressive 
corner  were  present.  All  the  data  computed  on  characteristics 
upstream  of  and  including  BC  is  referred  to  as  region  X  and  is 
stored  for  later  use  by  the  program. 

With  the  specified  turning  angle  and  the  con¬ 
ditions  in  front  of  the  shock  at  B,  the  shock  angle  and  the  con¬ 
ditions  behind  the  shock  can  be  determined  (see  Part  C-c) .  The 
shock  will  always  have  a  slope  less  than  the  local  slope  of 
the  down  running  characteristics  in  region  I,  hence,  the  shock 
will  always  be  in  region  I.  The  intersection,  E,  of  the  shock 
and  the  next  up  running  characteristic  of  region  I  is  determined 
and  the  conditions  upstream  of  the  shock  are  determined  by  inter¬ 
polation  between  points,  0  andP.  In  region  1+1,  the  region  behind 
the  shock,  FE  is  the  down  running  characteristic  from  the  upper  wall 
to  the  point  E.  Jt  is  assumed  that  the  conditions  at  F  are  equal 
to  the  conditions  behind  the  shock  at  B.  Then,  the  flow  behind  the 
shock  and  the  shock  slope  at  E  can  be  determined  by  simultaneously 
solving  the  compatibility  equation  (Eq.  13)  along  FE  and  the  Rankine- 
Hugoniot  equations  across  the  shock. 
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I  th  Strip 


Down  Running  Characteristic  (Region  I) 
Up  Running  Characteristic  (Region  I) 

Down  Running  Characteristic  (Region  II) 
Up  Running  Characteristic  (Region  II) 


FIGURE  16  -  CHARACTERISTIC  MESH  AT  COMPRESSION  CORNER 

ON  THE  UPPER  WALL 


With  the  known  conditions  behind  the  shock 
at  E,  the  intersection,  H,  of  the  up  running  characteristic  from 
E  with  the  upper  wall  can  be  determined.  The  flow  field  at  H 
can  be  determined  by  solving  the  compatibility  equation  (13b) 
along  EH  using  the  known  wall  slope  at  H» 

___  How  we  work  along  the  down  running  character¬ 

istic,  HJ,  in  region  I+l„  if  the  shock,  BD,  crosses  the  next  up 
running  characteristic  of  region  I  at  a  point  K  before  it  inter¬ 
sects  HJ,  then  the  conditions  in  front  of  the  shock  at  K  are 
determined  by  interpolation  in  the  characteristic  mesh  of  region 
I.  By  interpolation,  the  point  M,  on  EH,  whose  down  running 
characteristic  passes  through  K,  and  the  data  at  M  are  determined. 
The  flow*  behind  the  shock  at  K  is  determined  by  simultaneously 
solving  the  compatibility  equation  along  MrC  with  the  Rankine 
Hugoniot  jump  conditions  across  the  shock  at  K. 

Then  the  location  and  conditions  at  L  are 
determined  by  a  regular  characteristic  computation  between  H  and 
K.  This  process  continues  until  the  down  running  characteristic, 
HJ,  intersects  the  shock,  BD,  at  a  point  J.  The  conditions  in 
front  of  the  shock  at  J  are  determined  by  interpolation  as  above; 
the  conditions  behind  by  using  the  compatibility  equation  along 
LJ  and  the  jump  conditions  at  J. 

The  result  is  that  the  flow  field  is  divided 
into  a  number  of  regions,  or  strips,  the  boundaries  of  which  are 
down  running  shocks. 


3.  In  the  Presence  of  Detached  Shocks 

It  may  happen  that  shocks  are  formed  by  the 
coalesence  of  compression  waves  due  to  a  smooth  turning  of  either 
of  the  walls  (see  Part  D-f-2) .  The  formation  of  up  running  shocks 
by  such  a  process  in  no  way  affects  the  progression  of  the  com¬ 
putation  on  the  gross  scale  now  under  discussion. 

However,  if  a  down  running  shock  is  formed  by 
the  coalescence  of  compression  waves,  we  again  divide  the  flow 
field  into  two  strips  as  illustrated  in  Figure  17. 
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Down  Characteristic 
\  of  i+l  strip 


Last  Down  Characteristic 
of  I  Strip 


/  s 


Down  Running  Shock  ^ 

Formed  by  the  Coalesence  of 
Compression  Waves 


FIGURE  17  -  DOWN  SHOCK  FORMED  BY  THE  COALESENCE 

OF  COMPRESSION  WAVES 
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4.  The  Presence  of  Centered  Expansions  on  the 
upper  or  lower  walls  causes  no  essential  change  in  the  pro¬ 
gression  of  the  computation  (see  Part  D-d). 

c)  Computation  of  Regular  Points 

Now  that  we  have  indicated  the  gross  order  of  pro¬ 
gression  of  the  computer  program,  we  can  look  into  the  various 
specific  computations  in  a  little  more  detail. 


1. 


Interior  Flow  Field  Points 


As  mentioned  in  Part  I  of  this,  section,  regular 
interior  points  are  determined  by  the  intersection  of  a  down  running 
characteristic  with  the  up  running  characteristic  issuing  from  a 
mesh  point  on  the  previous  down  running  characteristic.  This  is 
illustrated  in  Figure  18  where  C  is  the  point  to  be  computed.  All 
of  the  points,  including  B,  on  the  ( I— 1 )  st  down  running  character¬ 
istic  have  been  computed,  as  have  all  the  points  through  A  on  the 
i  th  characteristic.  Point  C  is  the  next  point  to  be  computed. 


First 

Equations 

13a  and  13b 

are 

written 

in 

difference  form: 

(Tc  -  ta> 

<PC  -  PA> 

-I 

V 

j 

(17a; 

-  Bac 

=  gac  (xc  - 

<Tc  •  V 

+  V 

(Pc  -  PB> 

=  5"  <*„  - 

BC  C 

V 

j 

(17b) 

where  the  barred  quantities  are  averaged  quantities  on  the 
appropriate  characteristic.  Equations  17a  and  17b  are  solved 
simultaneously  to  give 


b  p+b  p-t+t- 


pc  = 


AC  A  BC  B 


G*  (x  — x  )  —  (x  -X  ) 


B  l  AC  '  C  A'  BC  '  C  B  J 


b,_  +  b  _ 
AC  BC 


-I 


Tc  =  TA  +  bAC  <PC  '  PA>  +  GACtXC-XA>J 
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(13a) 

(18b) 


i  tft  Down  Runnine 
Characteristic 


(i-l)st  Down  Runninc 
Characteristic 


FIGURE  18  -  REGULAR  INTERIOR  POINT  COMPUTATION 


In  addition, we  have  the  following  definitions  and  equations 
available? 


sc  =  ln  pc  -  r  ln  pC'  or  pc  =  e 


ln(pc)-Sc 


(19a) 


Hc  =  A  Pc/pc  +  “2  '  °r  qc 


PC 

2<V^r^> 


(19b) 


2  2  2  2 

qc  =  uc  +  vc  ‘  or  uc 


v  2 

_fc_ 

1+T2, 


(19c) 


T 


c 


or  vc 


U  T 

c  c 


(19d) 


Now,  it  is  fairly  common  to  perform  a  regular  computation  in  some¬ 
what  the  following  form. 

1)  Get  the  location  of  C  by  the  intersection 
of  characteristics  from  A  and  B. 


2)  Compute  pc  and  tc  using  (18a)  and  (18b) 
where,  in  the  first  pass,  the  averaged  quantities  are  taken  as  those 
at  A  or  B,  whichever  is  appropriate. 

3)  Project  the  streamline  of  slope  tc  back 
from  C  to  intersect  the  line  AB  c-  meeting  A  and  B  (not  shown  in 
Figure  18) . 


4)  Since  S  and  H  are  constant  along  streamlines 
(in  the  absence  of  shocks)  obtain  Hc  and  Sc  by  linear  interpolation 
along  AB. 

5)  Use  equations  19  to  determine  uc,  v^,  pc 
from  the  known  PC'  Tc'  S 6C'  Hc* 

I  II 

6) '  Compute  the  coefficients  br.,  G_  ,  G_  ,  etc., 

average  the  coefficients  along  each  characteristic,  and  iterate 

until  the  characteristic  slopes  at  C  do  not  change  within  some  specified 

tolerance. 
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However,  if  one  proceeds  in  this  manner,  and  if 
the  flow  is  nonisentropic,  he  will  find  that  the  linear  inter¬ 
polation  in  step  4  introduces  errors  that  can  become  quite  large 
(Reference  4) .  These  errors  are  quite  evident  when,  after  computing 
a  number  of  characteristics,  one  checks  the  mass  flow,  where  the 
error  can  exceed  40%. 


A  much  better  way  of  proceeding  is  to  set  up  a 
mass  flow-entropy,  enthalpy  table.  The  fact  that  S  and  H  are  con*- 
stant  along  streamlines  (and  that  we  can  compute  the  change  in  S 
through  a  shock)  means  that  we  can  set  up  functions  defined  point- 
wise  for  S  and  H  with  mass  flow  as  independent  variable .  Then  we 
replace  steps  4-5  by  another  process. 

4-5  a)  guess  S  and  H  ,  compute  p  ,  u  ,.  v  . 

V<  V  Vo 

b)  compute  the  mass  flow  through  AC, 
subtract  it  from  the  known  mass  flow 
below  A  to  get  the  mass  flow  at  C. 

c)  go  to  the  mass  flow  -  S,  H  table  and 

obtain  S„,  H 
C  C 

d)  compare  with  the  guessed  values  of  Sc 
and  H  and  iterate  until  convergence. 

Vo 

Then  continue  with  6)  in  the  indicated  manner. 

The  result  is  a  very  accurate  computation  with  the 
error  in  the  integrated  mass  flow  often  less  than  1%  in  the  compu¬ 
tation  of  two-dimensional  flows.  Of  course,  the  mass  flow-entropy 
table  must  be  modified  when  there  are  shocks. 

The  dependence  of  S  and  H  on  mass  flow  are 
determined  initially  from  the  known  conditions  on  the  initial  data 
line. 


This  technique  is  used  in  the  present  program. 
The  increment  of  mass  flow  across  AC  is  given  by 
pq  sin  [tan  ^  t  -  tan  ^  X1]  A  area  (20) 
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where  the  barred  quantities  are  averaged  between  A  and  C  and  Aarea 
is  the  area  between  A  and  C,  given  by, 


Aarea  = 


Aarea  = 


J  (*c-*ft>2  +  (vV2' 

/  2  2 
7  1  +  X  (ya  -  yc)  for 


for  2-diin.  flow  (j=0) 


axi-sym.  flow  (j=l) 


(21a) 

(21b) 


2.  Boundary  Points 

Boundary  points  are  computed  in  the  same  manner 
as  interior  points,  except  for  a  few  simplifications.  First,  since  in  an 
entirely  inviscid  flow  the  boundaries  are  streamlines,  we  know  a  priori 
what  S  and  K  are  along  the  boundary  (the  case  of  shocks  is  explained 
in  Part  D-g),  and  we  need  not  do  the  mass  flow  -  S,  H  iteration  required 
at  interior  points.  In  addition,  tc  =  vc/uc  =  dy/dx  on  the  boundaries 
is  a  specified  function  of  x.  To  compute  a  point  on  the  lower  wall  we 
use  Equation  17a  along  the  down  running  characteristic  to  get  pc.  To 
compute  a  point  on  the  upper  wall  we  use  17b  along  an  up  running 
characteristic.  Naturally,  ,erations  are  made  in  averaging  coefficients 
when  the  equations  are  solv<-  . 

3.  Centerli  j  Points 

When  the  flow  is  axisymmetric,  it  may  happen  that 
the  lower  wall  ends  at  the  axis  of  symmetry  and  then  the  lower  wall 
computation  must  be  replaced  by  a  symmetry  cojdition..^  The  only  effect 
this  has,  in  the  absence  of  shocks,  is  that  G  and  G1  ,  which  take  the 
indeterminate  form  0/0  at  y  =  0,  must  be  evaluated  by  L* Hospital's  rule. 


,1,  II 


_  .  *:,n.  t, 

y=0  -  <T-X  >  7 


=  I, II  dT 
y=0  dy 


y=0 


(22) 


The  program  requires  a  cusp  at  the  trailing  edge 


of  the  lower  wall. 


There  is  no  provision  to  compute  the  reflection  of 
shocks  at  the  axis  of  symmetry. 
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d.  Centered  Expansion  on  the  Boundaries 

Everytime  we  compute  the  location  of  an  upper  or  lower 
boundary  point,  a  test  is  made  to  determine  whether  or  not  there  is  a 
corner  between  the  current  point  and  the  last  point  on  that  boundary. 

If  a  corner  has  been  detected,  the  next  operation  is  to  find  out 
whether  or  not  it  is  an  expansion  or  a  compression  corner.  The  latter 
is  treated  in  Part  D-f. 

If  it  is  an  expansion  corner,  the  conditions  just  up¬ 
stream  of  the  corner  are  first  determined  (see  Figure  19) .  Then  the 
total  turning  angle  ,  A0<  of  the  flow  is  determined.  This  turning 
angle  is  divided  into  a  number  of  increments,  n,  of  size  6  ( A )  =A0/n , 
and  the  conditions  corresponding  to  the  flow  which  has  turned  from 
the  state  just  prior  to  the  corner,  say  0,  to  9=0  +  m  6  (A0) ,  m=l,n, 
are  determined.  This  is  done  by  determining  the  change  in  Mach  number 
to  go  from  0  to  9  +  m  6  (A0)  from  Equation  15) .  Then  the  rest  of  the 
conditions  can  be  determined  by  using  Equation  16. 

If  the  expansion  is  on  the  upper  boundary,  a  new  down 
running  characteristic  is  put  in  the  mesh  for  each  value  of  m,  and  the 
regular  upper  boundary  point  computation  is  bypassed  until  all  of  the 
down  running  characteristics  issuing  from  the  corner  have  been  computed 
(see  Figure  19a) . 

If  the  expansion  is  on  the  lower  wall,  a  number,  n+1, 
of  mesh  points  is  put  at  the  corner  and  the  current  down  running 
characteristic  is  finished  by  making  a  series  of  regular  interior 
point  computations  as  indicated  in  Figure  19b. 

e.  Computation  of  Interior  Shock  Points 

1.  Up  Running  Shock  Points 

In  Figure  20,  AB  is  a  down  running  characteristic 
along  which  we  know  everything,  including  the  location  and  slope  of 
any  up  running  shocks  crossing  AB.  Such  a  characteristic  shock  inter¬ 
section  is  designated  by  two  points  (H,I  of  Figure  20)  which  have 
the  same  physical  coordinates,  one  having  the  properties  of  the  flow 
just  in  front  of  the  shock  and  the  other  properties  just  behind  the 
shock.  We  want  to  compute  CD,  the  next  down  running  characteristic. 
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Supposing  the  shock  to  be  the  first  up  running 
shock  to  cross  CD,  all  points  through  E  on  CD  can  be  computed  as 
ordinary  points.  As  each  such  point  is  computed,  we  test  to  see  if 
the  down  characteristic  and  the  shock  intersect  before  we  can  com¬ 
pute  any  other  regular  points.  When  such  an  intersection  occurs, 
we  locate  it  (M  in  Figure  20)  by  intersecting  the  down  characteristic 
EM  from  E,  the  last  good  regular  point ,  with  the  up  characteristic, 
F?>M,  from  F*,  an  interpolated  point  on  AB.  The  conditions  at  M  are 
determined  by  a  regular  characteristic  computation  between  F*  and  E. 


Now  we  must  determine  the  new  shock  slope  at  M 
and  the  conditions  behind  the  shock  at  N.  To  do  this,  a  point  J* 
on  AB  is  found  by  interpolation  such  that  the  up  running  character¬ 
istic,  J*N,  from  J*,  intersects  the  shock  at  N.  Then  the  compatibility 
equation  (131)  along  J*N  is  solved  simultaneously  with  the  Rankine 
Hugniot  conditions  across  the  shock  to  determine  everything  including 
the  shock  slfcpe,  at  N.  In  outline,  the  process  of  solution  is: 


1)  guess  a,  the  shock  slope  at  M7 

2 

2)  compute,  p,  p,  t,  M  at  N  from  Rankine 
Hugniot  equations  (14)  and  the  known 
conditions  at  M, 

3)  compute  p',  at  N  from  the  known  conditions 
at  J*  using  the  compatibility  equation 
(17b)  along  J*N  and  the  t  from  2,  above, 

4)  compare  p'  and  p;  if  they  do  not  agree, 
guess  another  value  of  a  and  repeat  the 
process  until  convergence. 

When  the  iteration  has  converged,  an  additional 
test  is  made  on  If  an  expansion  has  caused  the  shock  to 

become  so  weak  that,  within  the  numerical  accuracy  of  the  computation, 
the  shock  no  longer  has  finite  strength,  it  is  desirable  to  eliminate 
it  from  the  flow  field.  Thus  if  £  1,  we  drop  the  shock. 

2.  Down  Running  Shock  Points 


The  intersection  of  a  down  running  characterist c 
and  a  down  running  rhock  were  outlined  in  Part  D-b-2,  although 
the  manner  in  which  the  iteration  is  carried  out  was  not  given.  The 
iteration  is  performed  exactly  as  in  the  case  of  a  down  character¬ 
istic  up  shock  intersection  detailed  in  Part  D-e-1,  except  that  the 
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compatibility  equation  is  17a  written  along  a  down  running  character¬ 
istic  (see  Figure  16  and  Part  D-b-2)  ,  There  is  presently  no  provision 
to  drop  a  down  shock  whic;  has  become  very  weak.  This  should  not  cause 
any  trouble  in  the  use  of  the  program. 

f.  Hov:  Are  Shocks  Started? 

Shocks  may  arise  in  two  fashions,  differing  only  in  their 
smoothness.  There  may  be  sharp  compressive  corners  in  the  flow  field 
which  will  give  rise  to  attached  oblique  shocks  (Figure  7) (the  case  of 
strong  shocks  is  not  considered  here).  Or,  a  continuous  smooth  compres¬ 
sion  may  coalesce  into  a  shock  as  illustrated  in  Figure  9. 

1.  Formation  of  Attached  Oblique  Shocks  by  Sharp  Corners 

a)  Lower  Boundary  -  When  a  corner  is  encountered  on 
the  lower  boundary,  the  conditions  just  in  front  of  the  shock are  found 
by  an  ordinary  lower  boundary  computation  (Part  D-c-2)  along  MA  of 
Figure  21  where  M  is  an  interpolated  point.  The  turning  angle  of  the 
flow  is  specified  and  the  conditions  behind  the  shock  at  B  and  the  shock 
slope  at^(A,B)  are  completely  determined  from  Equations  14  where 

6  =  tan  Tg-  tan  t_.  The  intersection  of  the  shock  with  the  down  running 
characteristic  (E  of  Figure  21)  is  found  and  the  conditions  in  front  of 
the  shock  are  determined  by  a  regular  characteristic  compuation  between 
G  and  H  where  G  is  an  interpolated  point  on  the  previous  down  character¬ 
istic  and  H  is  the  last  computed  point  in  front  of  the  shock  on  the 
current  characteristic.  C  is  the  point  on  the  wall  behind  the  shock 
whose  up  running  characteristic  intersects  the  shock  at  (E, F) .  The 
conditions  at  C  are  assumed  to  equal  those  at  B.  A  regular  shock  com¬ 
putation  between  E,  F,  and  C  is  then  performed  (see  Part  D-e)  determin¬ 
ing  the  conditions  at  F  and  the  slope  of  the  shock  at  (E,F) .  Then  D 
is  computed  as  a  regular  boundary  point. 

b)  Upper  Boundary  -  The  computation  of  the  for¬ 
mation  of  shocks  on  the  upper  boundary  has  been  outlined  in  Part  D-b-2 
except  that  no  indication  was  given  of  the  equations  used  to  determine 
the  conditions  behind  the  shock.  However,  the  equations  are  those  in¬ 
dicated  in  the  preceeding  paragraph  and  there  is  no  need  to  duplicate 
here  a  description  of  the  computational  procedure  outlined  in 

Parts  D-b-2  and  D-f-1. 

2.  Shocks  Formed  by  Coalescence  of  Smooth  Compressions 

The  formation  of  oblique  shocks  by  the  coalescence 
of  a  smooth  compression  is  manifested  numerically  by  the  crossing  of 
characteristics  of  the  same  family  as  illustrated  in  Figure  9  for  a 
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compression  from  the  lower  wall.  In  the  program  such  shocks  are 
detected  by  testing  each  regularly  computed  mesh  point  to  see  if 
there  has  been  a  crossing  of  either  characteristic  family.  In 
such  an  instance  a  shock  of  that  family  is  initiated. 

g.  Reflection  of  Oblique  Shocks  at  the  Walls 

1.  Up  Running  Shock  Reflected  at  the  Upper  Wall 

The  reflection  of  an  up  running  shock  at  the 
upper  wall  is  illustrated  in  Figure  22.  It  is  easy  to  detect  such 
a  reflection  since  the  counter  indicating  the  mesh  point  of  the 
shock  on  the  down  running  characteristic  approaches  one  (1)  as  the 
shock  approaches  the  boundary.  When  such  a  reflection  occurs,  we 
proceed  in  a  manner  similar  to  that  used  when  we  have  a  sharp  com¬ 
pressive  corner  on  the  upper  wall.  Referring  to  Figure  22,  AB  is 
assumed  to  have  been  completely  computed.  The  location,  (C,D,F)  of 
the  intersection  is  determined  and  the  conditions  at  C  in  front  of 
the  incident  shock  are  determined  by  a  regular  upper  wall  computation 
(H  is  an  interplated  point  on  AB) .  Then  the  slope  of  the  incident 
shock  at  the  point  of  intersection  and  the  conditions  at  D  (behind 
the  incident  shock  but  in  front  of  the  reflected  shock) ' are  deter¬ 
mined  by  a  regular  up  shock-down  characteristic  computation  (Part 
D-e-1) .  The  characteristic  DE  is  then  completed  in  the  usual  manner. 

Now  since  the  up  shock  reflects  as  a  down  shock, 
the  latter  will  separate  the  flow  field  into  two  strips,  say  the  I'th 
s  .  .p  in  front  of  the  reflected  shock  and  the  (1+1)  st  strip  behind 
it  The  incident  shock  turns  the  flow  toward  the  wall,  and  we  can 
think  of  there  being  an  effective  compressive  corner  at  the  point 
of  reflection  which  turns  the  flow  from  its  direction  at  D,  between 
the  two  shocks,  to  the  direction  parallel  to  the  wall  at  F,  behind 
the  reflected  shock.  Thus,  after  DE  is  computed  the  reflection  is 
treated  exactly  as  a  sharp  compressive  corner  (see  Part  D-f-'-b). 

2.  Down  Running  Shock  Reflected  at  the  Lower  Wall 

The  configuration  arising  when  a  down  running  shock 
reflects  at  the  lower  wall  is  illustrated  in  Figure  23.  While  working 
along  the  down  running  characteristic,  DEF,  i  ,  the  I'th  strip,  tests 
are  made  to  see  if  the  characteristic  intersects  the  lower  wall  before 
it  intersects  the  down  shock,  GA,  separating  the  I'th  and  (I-l)st 
strips.  When  this  happens,  the  point  of  intersection  (A,B,C)  of  the 
shock  with  the  wall  is  determined.  The  conditions  at  A  are  determined 
by  interpolating  on  the  (I-l)st  strip.  The  conditions  at  B  are 
determined  by  malting  a  regular  down  shock  computation  as  previously 

42 


outlined.  The  incident  shock  turns  the  flow  toward  the  wail  at  B, 
between  the  two  shocks.  The  slope  of  the  reflected  shock  and  the 
flow  at  C  is  then  determined  by  requiring  the  known  flow  at  D  to  turn 
back  parallel  to  the  wall  at  C.  By  extrapolating  the  reflected  shock 
we  find  E,  the  last  good  computed  point  on  DEF.  By  a  procedure 
identical  with  the  case  of  a  sharp  compressive  corner  on  the  lower 
wall  (Part  D-f-1)  the  rest  of  DEF  is  recomputed  as  indicated  in  the 
figure. 


h.  Intersection  of  Shocks  of  the  Same  Family 

In  general,  such  an  intersection  will  be  continued  as 
a  shock,  a  contact  discontinuity,  and  a  reflected  expansion  or  shock 
as  is  illustrated  in  Figure  4.  We  neglect  the  contact  discontinuity 
and  reflected  expansion  or  shock  because  they  should  be  weak  compared 
to  the  shock  OC. 

i.  Intersection  of  Shocks  of  Opposite  Family 

This  situation  is  illustrated  in  Figure  24.  Such  an 
intersection  will  result  in  the  diffraction  of  the  two  incident  shocks 
and  the  creation  of  a  contact  discontinuity  along  which  the  pressure 
and  velocity  direction  are  continuous.  While  working  along  the  down 
characteristic  AB  in  the  Ith  strip,  tests  are  made  to  determine  whether 
or  not  the  down  shock  COB  separating  the  two  strips  intersects  an  up 
shock  of  the  (1-1) st  strip.  When  this  happens,  the  pdiht  of  inter¬ 
section,  O,  is  determined.  The  conditions  at  0  (in  front  of  both 
intersecting  shocks)  are  obtained  by  interpolation  on  the  (I-l)st 
strip.  The  slopes  at  0  of  the  two  incident  shocks  and  the  conditions 
at  E  and  F  are  then  determined,  the  slopes  being  assumed  equal  to 
the  slope  at  the  nearest  previous  shock  point  and  the  conditions  at  E 
and  F  being  then  computed  directly  using  the  Rankine  Hugoniot  rela¬ 
tions  . 


The  slope  of  the  diffracted  up  shock,  OK,  is  then 
guessed  and  the  conditions  at  G  are  determined,  including  t  .  Then 
the  slope  of  the  diffracted  down  shock,  OB,  is  determined  w^ich  gives 
tH:=tg  •  Then  pQ  and  p-^  are  compared  and,  if  they  disagree,  the  slope 
of  TJR  is  changed  an  an  iteration  is  performed.  Then  the  intersection, 
K,  of  the  down  characteristic,  AB,  with  the  diffracted  up  shock,  OK, 
is  determined  as  are  the  other  intersections  indicated  in  the  figure. 
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The  next  down  running  characteristic  in  the  Ith  strip 
is  then  made  in  the  normal  fashion. 

j.  Intersection  of  Shocks  with  Contact  Discontinuities 

The  situation  which  arises  in  this  case  is  depicted 
in  Figure  6  where  it  is  seen  that  the  intersecting  shock  and  contact 
discontinuity  are  diffracted  at  the  point  of  intersection  with  a 
weak  centered  expansion  often  arising  at  the  point  of  intersection,  0. 

Currently,  contact  discontinuities  are  dropped  from 
the  computation  upon  intersecting  a  shock. 
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THREE-DIMENSIONAL  SUPERSONIC  FLOW  COMPUTER  PROGRAM  ELEMENTS 


A.  Summary  and  Introduction 

Computer  programs  have  been  written  for  machine  com¬ 
putation  of  the  basic  flow  elements  which  compose  a  large  variety 
of  three-dimensional  supersonic  inlet  flow  fields.  A  combin¬ 
ation  of  these  unit  problems  can  be  used  to  design  and  analyze 
three-dimensional  inlets  as  described  in  Section  III.  The  pro¬ 
grams  are  written  in  Fortran  IV  for  use  on  the  IBM  7094  high 
speed  digital  computer. 

The  overall  program  can  be  summarized  as  follows: 


General  -  Computer  calculation  of  the  basic  elements 
of  three-dimensional  supersonic  inlet  flow 
fields 


Input 


Output 


Plane  shock  -  conical  body 
Conical  shock  -  plane  body 
Conical  shock-conical  shock 
Delta  wing 

Uniform  free  stream  conditions 
Supersonic  leading  edges 
Any  combination  of  Mach  number,  body 
surface  angle  and  leading  edge  sweep  angle 
Perfect  gas 

Aerodynamic  flow  properties  along  the  axis 
of  symmetry  of  the  locus  of  the  three- 
dimensional  wave  intersections 
Aerodynamic  properties,  incident  and  re¬ 
flected  wave  angles  and  locations, 
immediately  behind  the  three-dimensional 
wave  intersections 

Complete  flow  field  behind  intersections 
calculated  for  plane  shock  -  conical  body 
and  conical  shock-  plane  body 


ii 


x 
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Complete  descriptions  of  these  program  elements  are  given 
in  this  section.  The  problems  are  divided  into  two  parts. 
Firstly,  the  swept  wedge  (Part  B)  which  is  considered  co  be 
infinite,  i.e.,  no  end  expansions.  This  problem  requires  the 
analytic  description  of  a  shock  envelope  and  flow  field  com¬ 
posed  of  a  central  conical  region  between  two  two-dimensional 
regions  which  extend  from  the  central  region  boundaries  to  the 
leading  edges.  Secondly,  for  the  remaining  problems  (Parts  C, 

D  and  E)  the  initial  shock  envelopes  and  flow  fields  are  either 
purely  conical  or  two-dimensional  hnd  therefore  =re  known  a 
priori. 

For  each  of  the  above  problems,  once  the ' initial  shock 
envelope  and  the  associated  flow  field  is  determined,  the  various 
three-dimensional  reflection  and  intersection  curves  are  com¬ 
puted  by  using  linearized  and  two-dimensional  supersonic  flow 
theory  and  working  in  planes  normal  to  the  tangent  of  the  inter¬ 
section  curve.  The  flow  v  elds  behind  the  intersection  curves 
are  then  constructed  from  two-dimensional  characteristics  with 
entropy  variations. 

It  is  suggested  that  the  reader  study  the  analysis  given 
in  Parts  B  through  E  before  using  these  programs,  in  order  to 
fully  understand  the  approximations  involved.  The  operating  in¬ 
structions,  including  the  input  formats,  definition  of  symbols, 
and  example  problems,  are  given  in  the  Appendices  to  this 
section. 
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S-l  Swept  Wedge  Flow  (Delta  Wing) 

According  to  linear  theory^  the  effect  of  the 
vertex  of  a  supersonic  swept  wedge  is  limited  to  a  region 
inside  the  free-stream  Mach  cone  from  the  vertex.  In  the 
analysis  presented  herein  the  two  dimensional  leading  edge 
shock  and  flow  field  (the  region  outside  this  central  Mach 
cone)  are  obtained  from  exact  oblique  shock  theory  applied 
in  a  plane  normal  to  the  wedge  leading  edge.  Since  an  ex¬ 
act  solution  for  the  region  influenced  by  the  vertex  Mach 
cone  is  nonexistent,  linear  theory  is  used  on  the  base 
plane  of  the  wedge  to  obtain  the  flow  properties  inside 
the  Mach  cone  corresponding  to  the  two  dimensional  Mach 

number.  This  cone  is  assumed  to  divide  the  two-  and  three- 

12 

dimensional  regions.  Similar  to  the  solution  of  Babaev 
and  Fowelll^  for  a  swept  flat  plate  at  an  angle  of  attack, 
a  second  order  curve  is  used  as  an  approximation  for  the 
shock  envelope  in  the  3-D  central  region.  The  flow  con¬ 
ditions  behind  the  reflection  of  the  wedge  shock  from  a 
plane  surface  are  then  determined  from  oblique  shock 
theory  in  planes  normal  to  the  tangent  of  the  intersection 
curve . 


The  supersonic  flow  over  a  swept  wedge  as  shown 
schematically  in  Figure  28  is  considered  here.  The  sweep 
angle  is  designated  by  and  the  wedge  angle  in  the  plane 

of  symmetry  is  denoted  .  The  leading  edges  are 

assumed  to  be  supersonic  ~ud  the  wedge  is  taken  to  be  of 
infinite  span. 

As  is  well  known^"\  the  influence  of  the  wedge 
apex  (O')  is  only  felt  inside  of  a  Mach  cone  whose  vertex 
is  at  O'.  Outside  of  this  cone  the  flow  field  behind  the 
leading  edge  shock  is  two-dimensional.  This  region  may 
be  analyzed  in  the  usual  way  by  breaking  the  free  stream 
Mach  number  into  two  components,  tangent  and  normal  to 
the  leading  edge.  Oblique  shock  theory  can  then  be  used 
in  the  plane  normal  to  the  leading  edge  to  determine  the 
flow  behind  the  two-dimensional  portion  of  the  shock - 
The  wedge  angle  in  this  normal  plane  is  determined  geo¬ 
metrically  to  be 


X 

n 


tan 


-1 


tan  X 
cos  x 


(23) 
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and  the  flow  deviation  on  the  wedge 

to  line  O'C  (Fig. 28),  is 


n  r  f  clC O  v*/ 1 1  h  2T6  sp(? c  t 


ivI 


a  --  cos  ^  (cos  X  sin  x)  •  tan 


n 


1 


,M  I 
1 


(30) 


The  general  shape  oi  the  swept  wedge  shock  consists  of 
the  two  plane  shocks  attached  to  the  leading  edges,  as 
described  above,  connected  through  the  central  region  by 
a  curved  section  which  is  determined  by  the  three- 
dimensional  r...  cure  oi  the  fl">w,.  In  under  to  analyze  this 
central  region  it  is  necessary  to  make  some  approximations 
concerning  the  shock  shape  and  the  distribution  of  the 
flow  properties.  The  local  Mach  cone  defined  by  Ug, 
with  the  vertex  at  O'  and  an  axis  along  O'C  is  assumed  to 
separate  the  two-dimons Lonal  flow  field  from  the  central 
portion,  which  is  conical  under  the  assumption  of  linear 
theory-^.  Due  to  the  conical  nature  of  the  central  flow, 
the  element  of  the  wedge  shock  in  the  vertical  plane  of 
symmetry  (x  =  0)  is  a  straight  line  and,  therefore,  a 
linear  function  of  z ..  In  a  plane  normal  to  the  z'  axis 
the  shock  is  assumed  to  be  a  parabola,  which  is  tangent 
to  the  two-dimensional  shocks  at  their  points  of  inter¬ 
section  with  the  local  Mach  cone  (Fig.  29).  A  parabola 
is  chosen  since  this  family  of  curves  have  asymptotes, 
which  is  desired  si>.r?  the  curve  must  match  with  the 
straight  two  h  imons  i  o  ial  v.hook .  The  angles  that  the  two- 
dimensional  shocks  make  with  the  base  plane  of  the  wedge 
in  the  y-z  ana  x-y  planes  are  respectively: 

ft* 

0 

S 

Using  these  r< 

the  plane  normal  to  the  z'  axis  can  be  found  as 


- 1 

-  can 

(can  0]a  cos  \) 

(31) 

.  -1 
~  tan 

(tan  8ln  sin  x) 

(32^ 

fat ions 

the  two-dimensional  shock  angle 

:  in 

'  -1 
0  -  tan 
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T 
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cos  X  tan  (j9*-X)  +  sin  x]  tan  Xj  (33) 

The  points  of  intersection  of  the  two-dimensional  shocks 
with  uhe  local  Mach  cone  in  the  plane  z'  =  cos  X  are 
given  by  x'.  =  ±  cos  0  cos  X  <tan(J S*-X)  sin  0^  + 

r  o  3  ft  c  l  c 

[tan  Ui-tan2()3  -X)cos^  0^J  J 


y!  = 

3 


tan  0  -  cos  X  tan  (/3  -X) 

c 

z !  -  cos  X 
3 


(34) 


55 


Parabolic  Surfac  . 


FIGURE  29  -  SHOCK  CONFIGURATION 
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Since  the  central  portion  of  the  wedge  shock  is 
represented  by  a  parabola  tangent  to  the  two-dimensional 
shocks  in  a  plane  z'  =  const,  and  is  linear  in  z',  the 
three-dimensional  shock  envelope  is  represented  by  a 
parabolic  sheet  whose  equation  is 

x' ^  =  Az'  (y '  +  hz' )  (35) 


The  constants  A  and  h  are  determined  by  requiring,  as 
stated  previously,  that  the  shocks, pass  through  the  point 
(*!,  4j)  and  have  slope  tan  9C  in  the  plane  z'  = 

const,  at  that  point.  The  resulting  values  of  these 
constants  are 


A  = 


iH 

tan0' 

c 


1 

cos  A 


h 


A  cos  A 


1 

cosA 


(36) 


The  constant  h,  as  can  be  seen  from  -(36),  is  the  tangent 
of  the  angle  0SW  that  the  three-dimensional  shock  makes 
with  the  z'  axis  in  the  plane  of  symmetry  (y-z) . 


Linear  theory  has  been  used  in  the  conical 
region  to  obtain  the  variation  of  the  flow  properties 
from  their  uniform  values  in  the  two-dimensional  region 
to  the  plane  of  symmetry.  For  this  present  work  the 
variations  of  the  properties  in  the  y  direction  has  been 
neglected.  Therefore  the  variations  used  are  truly  valid 
only  on  the  baseplane  (x-z*)  of  the  wedge.  However,  for 
small  wedge  angles  (<10  degrees  approximately)  which  are 
generally  used  in  supersonic  and  hypersonic  inlet  designs 
the  approximation  is  good  for  the  whole  region.  Since 
the  local  Mach  cone  has  been  used  previously  to  construct 
the  shock  and  linear  theory  uses  the  free  stream  Mach 
cone,  the  conical  rays  of  linear  theory  must  be  adjusted 
to  fit  into  the  local  Mach  cone.  Thus,  given  the  ratio 
(x/z)^,  which  represents  the  slope  of  a  conical  ray  .i, 
the  adjusted  parameter  used  in  linear  theory  is 


(37) 
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where  n  =  tan  %  v  M*  -  X  .  The  resulting  variations  of 
pressure  and  Mach  number  in  the  conical  region  are 


P  1  +  H  V  »T  C 

1  _  00  pi 


1  +  \  y  M  C 

‘  '  no  v 


M  .  M  / - — 

I"-? 


°pi)2  2 

L-  T  +  ui 


tan  X 


where  the  pressure  coefficients  from  linear  theory  for 
the  two-dimensional  and  conical  regions  are  respectively, 


2  tan  X 


'(m2  -  1)  (1  -  n2) 


C  .  =  C  1  -  —  sin 
px  po  ff 


2  J2 

.  -1  n  -Ci 

in  — r 

U-c: 


and  the  lateral  velocity  u^  (x-direction)  is 


n  tan  X  -1  n  ~Ci  -1  n  +Ci 

u.  »  / —  <cos  — ~ -  cos  — ~ — y\ 

1  wCJ  Lna-qJJ  L“d+t.a 


]1 

n  +C.  I 

i  (42) 


The  flow  deviation  with  respect  to  the  z'  axis  in  the 
conical  region  is 


a  -  o 
li  1 


tan  X 


C 

1  - 

2 


where  tx^  (Eq.  30).  represents  the  flow  deviation  in  the 
two-dimensional  region. 
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B-2  Swept  Wedge  Shock  Reflection 


A  swept  wedge  shock  reflecting  from  a  plane 
surface  is  considered  next.  The  three-dimensional  wedge 
shock  intersects  a  plane  body  along  a  curve  which  ends 
in  two  straight  lines  (Fig. 30),  corresponding  to  the 
intersection  of  the  two-dimensional  portions  of  the  shock 
with  the  plane.  For  the  case  under  consideration,  where 
the  plane  body  is  parallel  to  the  base  plane  of  the  wedge, 
the  two  straight  line  sections  are  parallel  to  the  lead¬ 
ing  edges  of  the  wedge.  The  central  curved  portion  of 
the  intersection  curve  is  represented  by  the  second  order 
equation 


x.  =  ± 
1 


A  (Bz.  +  2Cz .  +  D) 
l  l 


1/2 


(44) 


where 


B  =  cos  X  (sin  X  +  h  cos  X) 

2 

C  =  r  (h  sin  X  cos  X  +  sin  X  -  1/2) 

2 

D  =  r  sin  X  (h  sin  X  -  cos  X)  (45) 


and  where  r  is  the  distance  between  the  plane  and  the 
base  of  the  wedge.  Values  of  z.  are  obtained  by  in¬ 
crementing  from  the  value  of  z^_Q  obtained  from  Eq  (44) 

with  x.  =0,  i.e.,  z.=z  +i  Az. 

1=0  l  o 


For  each  point  on  the  curve  the  angle  between 
the  tangent  to  the  curve  and  the  z  axis  is  given  by 

C) 


=  tan 


-1 


rA(Bz±  + 


(46) 


When  the  value  of  reaches  (ff/2)-  x>  the  intersection 
curve  becomes  a  straight  line  and  the  incrementation  of 
z^  can  be  stopped. 


The  intersection  curve  represents  the  line 
along  which  the  shock  reflects  from  the  plane  body.  .The 
flow  field  between  the  incident  and  reflected  shocks  has 
been  assumed  to  be  the  same  as  the  flow  on  the  wedge  sur¬ 
face,  i.e.,  no  variation  in  the  y-direction  (Fig. 30).  In 
order  to  obtain  the  flow  conditions  behind  the  reflected 
shock,  it  is  necessary  to  use  the  oblique  shock  relations 
in  planes  normal  to  the  tangents  of  the  intersection 
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curve  for  each  point  i.  To  simplify  the  procedure  the 
components  of  the  Mach  number  between  the  incident  and 
reflected  shocks  in  the  x-y-z  directions  are  found  first. 


- 

M.  .  =  M.  .  c 

1  .  li 
xi 


cos  x  cos(x  "  an)  “  cos  sin  x  sin(x  - cr ,.) 


n 


'li' 


M  =  -Mu  sin  X  sin(X  -  0  > 

yi 


(47) 


*  ^ 

Mj_  =  Mi.  sin  X  cos(X  -  )  +  cos  X^  cos  x  sin(x  -cr^) 


where  the  angle  between  the  leading  edge  and  the  z'  axis  on 
the  wedge  surface  is 

*  -1 

X  =  cos  (cos  X  sinx) 

Then  the  Mach  number  components  normal  and  tangent  to 
the  intersection  curve  on  the  plane  surface  are  easily 
found  from 


M1  =  M1  sin  K  ~  M-  cos  K 

n.  z,  ‘Lx. 

ii  x 

=  M^  cos  +  M,  sin  x/ 


t .  z . 

x  x 


\l  1  , 
x. 

X 


(48) 


(49) 


The  total  Mach  number  in  the  plane  normal  to  the  inter¬ 
section  curve  is 


V  *  <  +  Mi  > 

1  1  Yi 


1/2 


(50) 


and  the  corresponding  flow  deflection  is 


6^  =  tan 

i 


-1 


M, 
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(51) 
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These  relations  give  the  flow  conditions  immediately  be¬ 
hind  the  curved  portion  of  the  intersection  curve  for 
each  point  i.  The  point  i  for  which  the  value  of  X.' 
reaches  (ff/2  -  gives  the  conditions  behind  the  two- 
dimensional  shock  reflection. 

The  flow  conditions  in  the  whole  region  of  the 
plane  surface  behind  the  shock  intersection  can  be  found 
by  using  two-dimensional  characteristics  which  take  the 
variations  of  entropy,  due  to  the  curvature  of  the  shock, 
into  account.  This  method  has  been  applied  to  the 
conical  flow  problems  of  Sections  C  and  E  and  has  been 
presented  in  such  a  manner  (separate  subroutine  in  the 
computer  program) ,  that  it  can  easily  be  used  in  this 
wedge  program. 

Other  methods  of  solution  for  supersonic  flow 
over  the  wedge,  which  are  more  complicated  than  the 
approximation  of  linear  theory  used  here,  include  an 
exact  solution  similar  to  the  indirect  method  for  a  swept 
flat  plate  at  an  angle  of  attack  of  Ref.  12.  Also,  the 
present  solutions  can  be  improved  by  using  three- 
dimensional  linear  theory  to  give  an  approximation  which 
allows  variation  of  the  flow  properties  in  the  y-direction 
A  method  for  determining  the  y-variation  of  the  flow, 
which  should  yield  a  better  approximation  than  three- 
dimensional  linearized  theory,  is  to  use  the  conical 
flow  variation  between  the  body  and  shock  of  an  equiva¬ 
lent  cone,  whose  surface  pressure  is  matched  to  the 
surface  pressure  of  the  wedge  at  the  same  free  stream 
Mcch  number.  In  any  case,  the  methods  of  the  present 
work  and  the  analytic  techniques  presented  above  can 
either  be  used  directly  in  an  inlet  design  and  analysis 
effort  or  as  a  starting  point  for  additional  investiga¬ 
tions. 


The  program  operation  and  example  problems 
are  described  completely  in  Appendices  V  and  VI. 
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C-l  Impingement  of  a  Conical  Shock  on  a  Plane  Surface 


The  configuration  examined  here  is  the  super¬ 
sonic  flow  over  a  cone  of. half  angle  yb  at  zero  angle 
of  attack,  in  the  presence  of  a  plane  surface  inclined 
by  an  angle  8^  with  respect  to  the  cone  axis  (Fig. 31). 


The  intersection  curve  of  the  conical  shock,  (vc) 
produced  by  the  cone,  with  the  plane  surface,  (x',  z' 
plane) ,  is  represented  by  the  second  order  equation 


x! 

1 


(az:  -2Bz: 
1  1 


+  c  ) 


(58) 


where 

A  =  sin2  0i  -  tan2y  cos20x 

c 

B  =  r  sin0i+q  tan2y  cosGj 

'c  (59) 

C  =  r3  -  q2  tan2y 
'c 

where  0 i  is  the  inclination  of  the  plane  surface  to  the 
free  stream  direction  and  where  r  and  q  are  the  coordinates 
of  the  origin  O',  which  is  assumed  to  be  in  the  y’,  z' 
plane  (Fig.  31) .  Values  of  z!  are  obtained  by  increment¬ 
ing  from  the  value  of  z!_  ,  obtained  from  Equation  (58) 
with  x^  =  0,  i.e.,  z^=  z^°+  i  Az' 


For  each  point  of  the  intersection  curve,  the 
angle  between  the  tangent  to  the  curve  and  the  z'  axis 
is  given  by 

B-Azj  I 


The  intersection  curve  is  cut  off  when  ^becomes  5°  or 
less  in  order  to  work  on  only  one  side  of  the  cone.  The 
conical  flow  field  of  a  cone  at  zero  angle  of  attack  is 
completely  defined  by  the  values  of  the  critical  Mach 
number  M*  and  the  flow  direction  angle  ^  for  each 
conical  ray  angle  y  ^^.The  Mach  number  components  along 
and  normal  to  a  conical  ray  of  angle  y  are,  from  Ref.  15. 


where 


M  =  u/a 
r 

Mfc  =  Mrtan  ($x 


-  y) 


_ _ l 

[l+tans  (0X 
64 


(61) 


and 


+  h.1  M*2  )  % 

+  y+1  M1  } 


(62) 


represent  the  velocity  along  the  conical  ray  and  the 
speed  of  sound  non-dimensionalized  with  respect  to  the 
limiting  velocity,  respectively.  The  Mach  number  com¬ 
ponents  immediately  behind  the  conical  shock,  which  is 
represented  by  the  conical  ray  y  =  yc,  are  found  from 
Eqs. (61)  and  (62)  by  setting  y  =  yc.  In  order  to  calculate 
the  flow  field  beind  the  reflected  shock,  the  oblique 
shock  relations  are  used  in  a  plane  normal  to  the  tangent 
of  the  intersection  curve  for  each  point  i.  The  components 
of  the  Mach  number  immediately  behind  the  conical  shock 
in  the  x',  y' ,  z'  directions  are  found  first; 


M2  ,  =  M  siny  cos <p.  +  M  cosy  simp. 

X  ,  3T  C  X  L  C  1 

X 


Ms  ,  =  M  (sir.y  cos  <p.  cos©!  +  cosy  sin0x) 
y^  r  l  c 

+  M^  (cosyc  coscp^  cos0x  -  sinyc  sin0x ) 

M3  ,  =  M  (cosy  COS0J.  -  siny  cos <p.  sin0x) 

Z  .  JC  C  *.  C  X 

1 

-  Mfc  (cosyc  COS0-  sinfij^  +  sinyc  eos©x) 


where 


<p.  =  tan  1  (xl/yl)  with  y!  =  r-z!  sin0x 

*i  l  i  Ji  l 


(63) 


Then  the  upstream  Mach  number  components  normal  and 
tangent  to  the  intersection  curve  on  the  plane  surface 
are  found  from 


M8  =  Ms  .  sinfl.  -  M3  .  cosQ. 

n.  z!  i  *x!  l 

ii  i 

Ma  =  -Ms  ,  cosO.  -  Ma  ,  sinQ, 

L  ,  Z.  X  X  .  X 

11  1 


(64) 


The  total  Mach  number  and  the  flow  deflection  in  the 
plane  normal  to  the  intersection  curve  are 

Ms^  =  (Mgn^  +  M 

68i  =  tan"1  (Ms  ,/Msn  ) 

Yi  i 

With  the  oblique  shock  relations  giving  the  change  in  flow 
properties  across  the  shock  as 
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(65) 
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N. 
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and  the  relation 


(a2/a4>i 


(67) 


the  flow  properties  behind  the  reflected  shock  can  be 
found  from 


(68) 


p4  /Poo  =  <VP2}i  pc/p- 

i 

a.  =  it/ 2-0.  +  tan-1  (M.  /M.  ) 
4 .  l  4  4 

l  t .  N. 

l  l 


e„  =  e.  -  6 

4.  4.  2. 

ill 

_  2y^2  sin2e4  "  ^"1) 

s  =  s  /R+  ~=z  lln  (  - Si - i -  ) 

4.  2y  2  2y  “  +  1 


(y+l)M^  sin  6^ 

__  N.  i 

-  y  ^n  (  - i -  )  ] 

(y-l)M2  sin264  +2 
N.  i 

l 


(69) 

(70) 

(71) 


(72) 


where  pc/pa>  and  S2/R  represent  the  pressure  ratio  and  the 
entropy  increase  across  the  conical  shock.  The  shock  angle, 
0^  is  referred  to  the  plane  surface, 
i 
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C-2  Characteristic  Subroutine 


The  flow  properties  on  the  plane  surface  behind 
the  reflection  curve  are  determined  by  using  two-dimensional 
Characteristics.  The  flow  conditions  immediately  behind 
the  reflected  shock  on  the  plane  surface  are  calculated 
by  the  methods  of  the  preceding  section.  Due  to  the 
conical  flow  in  front  of  the  reflected  shock,  the  properties 
behind  the  shock  vary  in  the  y-direction.  Therefore,  the 
flow  is  truly  three-dimensional.  However,  the  y-variations 
are  small  for  most  cases  of  interest  in  inlet  design  (  yt  <  10 
approximately)  and  thus  a  two-dimensional  characteristic 
calculation  on  the  plane  surface  has  been  assumed  appropriate. 

On  the  plane  surface,  whose  coordinates  are 
redefined  as  p,  ip  the  intersection  curve  is  considered  a 
data  line,  since  the  flow  conditions  immediately  behind 
it  are  known  (Fig. 32).  Due  to  the  flow  symmetry  about 
the  plane  x=x'=0,  the  p  atxis  must  be  a  streamline  and 
therefore  can  be  considered  as  a  wall  with  zero  flow 
deviation  (6=0).  The  data  line  terminates  when  either 
a  five  degree  slope  or  a  Mach  reflection  is  reached. 

The  program  computes  the  characteristic  mesh  between  the 
data  line,  the  wall  and  the  second  family  characteristic 
line  from  the  last  data  point.  The  characteristic  equa¬ 
tions  include  the  effect  of  entropy  due  to  the  varying 
strength  of  the  reflected  shock.  This  subroutine  can  also 
be  applied  to  find  the  flow  field  behind  the  intersection 
curves  of  any  of  the  problems  considered  in  this  report. 

The  data  required  for  the  characteristic  equa¬ 
tions  at  each  point  (  p^,ipi  )  on  the  data  line  are:  the 
flow  deviation,  0  ,  with  respect  to  the  p  axis,  the  local  ... 
entropy  fh=(S.-S  )/2c  ,  the  local  velocity,  w. ,  non- 
dimensionalized  with  respect  to  the  local  limiting  velocity 
and  the  local  Mach  angle  H..  The  values  of  w.  and  |i  . 
can  be  found  from  the  known  Mach  number  behind  the  inter¬ 
section  curve  from  the  relations: 


|ii=  sin 


-1 


(  1/M..) 


(73) 
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The  equations  which  determine  the  location  of 
a  typical  point  Qj  in  Fig.  33  are 

^-^l  =  (P3“P1)  [tan  (6  +  U)  31<3  (74) 

^3”^2  +  ^P3~P2^  ^tan  ^3 

where  the  double  subscript  indicates  the  average  value 
between  the  two  points.  The  relations  which  hold  along 
the  first  and  second  family  characteristic  lines  are, 
respectively 

<  -  -> !.  3  <VW1>-(V®1,“  '  FI  =°^>,,3  ■ 

■  ¥■  qTq  (75)  ■■ 

dn  ulw3 

,  cotan  |i,  ,  .  2  .  .  2  .  ds  - 

<  w  "  “»^2,3  S  «2«3 


where  ds/dn  represents  the  entropy  variation  between  the 
points  and  Q normal  to  the  streamlines. 

The  characteristic  grid  consists  of  first  and 
second  family  characteristics  designated  by  j=l, 2. . . . 2N-1 
and  i=l,2....N-l  respectively,  where  N  is  the  number  of 
data  points  (Fig. 32.)  If  the  Eqs. (74)  are  rearranged 
into  a  form  suitable  for  computation  there  results 

p.  -  jWj  ~  ♦i.j-rpm.jt“‘®^)+  Pj,j-i 

tan  (  e- ^  )  -  tan  (0+p) 


*i,j=  +  (pi,j  -  tan  <9-p> 

where 

tan ( 0+u)  =^5 [tan ( 0 .  . +P ,  .)+tan(0.  .+|i.  .)] 

i+l,3  1+1,3  1/3  i/3 

and 

tan(8-|i)=is[tan(  04  j  1  -p.  .  )  +  tan(0.  .-p.  .)] 

i/3-l  i/3  i/3 

For  the  first  approximation  of  these  averages  only  the 
values  at  the  known  point  are  used.  The  order  of  compu¬ 
tation  is  to  start  at  the  point  i=N,  j=N  and  to  move 
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along  second  family  characteristic  lines,  (i=const.), 


rrom  tne  aaca  line,  toward  the  wall.  This  can  be  repre¬ 
sented,  in  subscript  form,  by 


i  =  N,  N-l,  N-2, . . .1 
j  =  i,  i+1,  i+2,...2N-l 

A  wall  point  is  obtained  whenever  ,i=2N-j.  At  these 

points  0.  .  =  6  =0  and 

*  •  J  1  /  J 

+  otan 

where 

ctan ( 0-|i) -h [c tan ( 0 .  .  .  _)-ctan  p.  .] 

1,3-1  1,3-1  i,D 

Eqs.  (75)  ,  solved  for  w  and  0,  give 


w.  .  =  i  .  +  a  .  J 

1,3  A  1+1,3  i,  3-1 

e  -I[0  _0  (SUL 

i,:  a  w.+1  .  |  w1<j_1  )] 

where 

a  .  .=  -  (sin2d  .cosp  .)  ("•) L.  . 

1+1,3  1+1,3  i+l,3  dn  i+l » 3 


-0 .  .  . +w . 

1+1,3  1+1,3 


ctanp .  ,  . 
1+1,3 

w.  , 

1+1,3 


a  .  =  — —  (sin2p.  .  cosli.  .  )(••—)  L.  . 

i,j-l  —  ^  i,j-l  i,j-l  dn  i,j-l 

(ctanp.  ~  \ 

- 1 '  J  J 


ctanii.  ,  .  \  ctanU.  .  , 

A  = \+  - 2-Jzl 

w.  ,  .  |  w. 

,  i+l,3  /  1,3-1 


and  where 
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Li+1, jt(*i, j"*i+l, j}  +  (oi,j'pi+l,  j*  ^ 


L.  .  =[(#.  .-0.  .  J2+(p.  .-p.  .  J2]*5 

1,3-1  1,3  1,3-1  1,3  1,3-1 


■jds^_  1,3-1  1+1,3 

Adn  - 
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2  2 
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i,3  i,3 


jctanVi,i 
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\  -  V  1+1,3 

ctan|i.  . 
1,3 

l  v,3 

2  w 

/  1+1 ,  j 

+ 

w.  . 
i,3 
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The  Mach  angl-;  and  the  non-dimensionalized  entropy  are 

r  _  i-w2 

li.  =sin_1  ^ 

X'3  L  2  wi,j  ' 


(84) 


S.  .  =S .  .  -  ~  L.  (sinli.  ) 

1,3  1,3-1  dn  1,3-1  1,3-1 


(85) 


At  the  wall,  (i=2N-j),  the  flow  properties  are 

w  =v/  +  =dL=====  [0.  .  + 

i/j  i-j-1  | ctantij^  1  i,3"l  -_i 

_ ;  wj < j - 1  1 

X  (si^li^.^cosM.JC  ^  )  L._  ._J 


(86) 


(87) 


6.  .=0 
1  1 

where  '  J  S  . 


,ds.  =  ~i,i-l~SN,N 
W  Li(  j_i(sin^<  j_x) 


Since  the  wall  i_s  a  streamline,  the  entropy  has  the 
constant  value  S  and  the  pressure  distribution  can 

be  found  from 


1 


1 


2 

WN,N 


(88) 


For  this  particular  problem,  p^.  ^/p^  is  p 

The  program  operation  and  example  problems  are 
given  in  Appendices  v  and  VI . 
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D  Intersection  of  Two  Conical  Shocks  of 
Different  Strengths 

D-l  Conical  Flow  Fields 

The  problem  examined  here  is  the  flow 
field  immediately  behind  the  intersection  of  two  conical 
shocks  of  different  strength.  The  conical  shocks  are- 
generated  by  two  unequal  cones,  one  of  which  is  at  a 
small  angle  of  attack  with  respect  to  the  free-stream 
velocity.  The  cone  axes  are  assumed  to  be  in  a  plane. 
These  assumptions,  which  simplify  the  analytical  solu¬ 
tion  of  the  problem,  are  consistent  with  common  geo¬ 
metrical  configurations  of  inlets. 

The  analysis  includes  the  definition  of 
the  intersection  curve  between  the  two  conical  shocks, 
the  definition  of  an  orthogonal  reference  frame  attached 
to  the  intersection  curve,  and  the  use  of  the  two- 
dimensional  oblique  shock  relations  in  the  plane  defined 
by  the  normal  and  binormal  to  the  intersection  curve. 

In  the  analysis  that  follows,  the  cone 
with  the  axis  parallel  to  the  free-stream  direction  will 
be  denoted  as  the  "base  cone"  while  the  designation 
"second  cone"  will  be  used  for  the  cone  at  a  small  angle 
of  attack. 


The  flow  fields  of  the  two  conical  bodies 
are  found  first.  For  the  base  cone,  the  velocity  along 
a  conical  ray  line  (  y=  const.)  in  a  spherical  coordinate 
system,  and  the  speed  of  sound  nondimensionalized  with 
respect  to  the  limiting  velocity  are  respectively^5 


m; 


/ 


l+tan‘ 


(iii^-y) 


a  =  (1-  f  , 

2  <y+l 


(89) 
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where  M^*  and)/).*,  which  represent  the  critical  Mach 
number  and  the  flow  direction, are  tabulated  in  Ref. 15. 
The  Mach  number  components  along  and  normal  to  the 
conical  ray  line  y  =const.  become  (Fig. 34) 


M  =  u/a 
r 

M,  =  M  tan  U*  -  y) 
t  r  V1 


(90) 


The  effect  of  a  small  angle  of  attack  on 
the  conical  flow  field  can  be  considered  as  a  perturba¬ 
tion  of  the  flow  conditions  for  no  angle  of  attack.  It 
has  been  shown  in  Ref.  11  that  a  circular  core  at  a 
small  angle  of  attack  produces  a  circular  conical  shock 
having  the  same  half  angle  as  the  shock  which  would  be 
produced  by  the  same  concial  body  at  no  angle  of  attack, 
but  inclined  slightly  with  respect  to  the  free  stream 
direction  in  the  plane  containing  the  body  axes  (z*)  and 
the  free-stream  velocity  vector,  Fig. 35.  Therefore,  the 
total  nondimensional  velocity  components  in  a  cylindrical 
coordinate  system  u,  v,  and  w,  along  the  z  axis,  along 
the  r  axes  and  normal  to  the  meridian  p.lane  respectively, 
may  be  defined  by  (-6) 


u^  =  u^  cosy  -  v^  siny  +  *  u2  cost^'' 


Vi  =  Vl  cosv  +  ui  sinv  +  a  v2  cos'-°i 


(91) 


w.  =  fli  W^  simv 

where  is  the  meridian  angle  defined  in  Fig.  35,  which 
is  incorrectly  defined  in  Ref,  16.  The  velocity  components 
u^  and  v^  are  the  zero-angle  of  attack  values  along  and 
normal  to  the  conical  ray  y-=  y,  , where  y_  is  the  shock 
angle,  and  are  given  by 


u 


1 


1 


) 


(92) 


where  M*  0  and  i/jf  ?  are  the  critical  Mach  number  and  the 
flow  direction  o"f '  the  second  cone  flow  with  no  angle  of 
attack.  The  perturbation  velocities  derived  from  Ref. 16 
are 


where  the  flow  direction  and  the  critical  Mach  numbers 

M*  and  M*  „  are  tabulated  in  Ref.  16. 

2  wf  2 

To  complete  the  definition  of  the  flow 
field  cf  the  second  cone,  the  local  speed  of  sound  and 
the  second  cone  shock  axis  inclination  must  be  defined. 
The  speed  of  sound  non-dimensionalized  with  respect  to 
the  limiting  velocity,  is 


where  the  total  critical  Mach  nuriber  M+,  is 

l 

Mi  =  M*  2  +  01  M2  COS  <Ai 

The  inclination  of  the  shock  axis  with 
respect  to  the  free-stream  direction  is  given  by 

=  a  (l-rj/w) 

where  r\/&,  tabulated  in  Ref.  16,  represents  the  ratio 
between  the  deviation  of  the  shock  axes  from  the  cone 
axis  and  the  angle  of  attack  (Fig. 36). 


D-2  Shock  Intersection 


In  the  preceding  section,  the  conical  flow 
fields  have  been  defined..  In  order  to  calculate  the  flow 
field  behind  the  intersection  of  the  two  conical  shocks, 
it  is  first  necessary  to  determine  the  intersection  curve. 

On  the  base  cone  a  conical  coordinative  system,  p,  $  is 
defined  as  shown  in  Pig.  37.  The  intersection  curve  of 
the  base  cone  shock,  of  angle  and  the  second  cone  shock, 

of  angle  •/  ,  whose  axis  is  inclined  by  an  angle  0^  to  the 
free  stream  direction,  is  represented  in  this  coordinate 
system  by 


'  ^*1 


where 


A 

B] 

C] 

°] 

E1 

and  where 

b! 


2  1  ^ 

-  (A^  cos  ^  +  2E^  cos j/k  +  B^) 

■1  "  [(a2-D  sln2v1]/F1 

=  [(Vbj2  )  coS2Vl]  /Pt 
=  [(e2b22  sinfli  ~  d2  cosfli^  sinv-jJ  /P^ 

-  [(e2b22  cos0l  +  d2  sin*V  cosv1  j /Fl 


b?  = 


d„  = 


e„  = 


f„  = 


F1  = 


(f2  sinyj^  cosy^/Fj^ 


tan  v-] 
tan2v. 


.2 


cos  0^  -  b2  sm  -7^ 


sin“01  -  b^  cos^  0^ 
r  cos0^  +  q  sin  0^ 
r  sin0^  --  q  cos  0^ 
cos0^  sin0^  (l+b|) 


.  2  ,2  2 

d2  ”  h2  ®2 


(97) 


(98) 


(99) 
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T**6  curve  is  computed  by  inserting 


.1  ..A.  .1  _  _ 

'  ax  UC  a  KJ1.  yj  .  yivcii  Jjy  , 


j/)£=iAi/),i=0/l,2*  •  •i/)ma3/A(|)into  Eq. (97)  where  0  max.  represents 
the  furthest  point  out  on  the  intersection  curve,  before 
the  curve  bends  in  aejain 

D,  C,  -E,  /  D,C,-E,  2  D2  -B, 


1 1)  =  cos 

max 


-1 


1~1 


D^-E^2  _  ^ 
(  5  I  )  ~ 


Ci  "Ai 


Cr  AX 


'1  *1  V  CX  “*i 

In  order  to  determine  the  planes  normal  to  the  tangent  to 
the  intersection  curve,  the  direction  cosines  of  the 
tangent,  bi-normal  and  normal  with  respect  to  the  x-y-z 
directions  must  be  found.  Following  the  procedure  of  Ref. 17, 
the  direction  cosines  of  the  tangents  to  the  intersection 
curve  are  found  as 

cos  (^x )  i  =  sinyi  (p!  sln2^  -  pi  cosO^/A* 


cos  (^y).  =  siny^  sini|)i(pi' aosi/H-pp/ A| 
cos(4:)i  =  p^  cosy^  sin0^/  A* 


(100) 


where 


_i  —  _  — — — 

i  smi/j, 


IfiL 


„  (crAi)c°5»i+DicrEi 

cos20i+2  (D1C1-E1)  cosj/)i+D2-B1 


and 


A*  -  [tp^)2  sin^+P?  «in2Vl] 


The  bi-normal  to  the  curve  is  the  normal  to  the  plane 
tangent  to  the  cone  at  the  generic  point  i ,  and  the 
direction  cosines  are  given  by 

cos  (b^)  =  cosy^  sinjl)^ 


cos  (£y)  i  =  cosy^  cos il)^ 

cos(b/z)  .  =  -sinv, 
any  i  1 


(101) 

(102) 


(103) 


Finally,  the  vector  normal  to  the  curve  can  be  found 
from  the  vector  product  n  =  t  x  b  and  lies  on  the  plane 
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tangent  to  the  cone  at  the  generic  point  i.  Its  direction 
cosines  are 

cos(nx).  =  cos(l^y).  cos(l3z)  .  -  cos(tAz).  cos  (by) 

1  1  1 

cos (ny) ^  =  cos (£z) ^  cos  (6x)  ^  -  cos  (£x)^  COs(l3Z)^  (104) 

cos(nz)i  =  cos(^x)i  cos(l3y)i  -  cos(^y)i  cos(bAx)i 


The  Mach  number  components^  immediately  behind  the  base 
cone  shock  (y=y.  ) ,  in  the  t,  b  and  n  directions  are  from 
Eq.  (89)  and  (90) . 

M~  =  Mx  cos(^x)i+M  cos(^y)i+Mz  cos(4)  . 

^lj_  X^i  Y^i 

Mg  =  Mx  cos(A)i+M  cos(ny)i+Mz  cosfnz^  { 

1i  li  1  Yli  li 

Mg  -  Mx  cos(bAx).+M  cosfiSy^+Mj,  cosfbz^ 

*i  li  X  yli  *i 

where 

M  =  (M  siny  +  M  cosy.)  sinj/j. 

X  4  }T  1  tl  X  1 

±i 

M^l  =  (Mr  siny1  +  Mfc  cosy^)  COSj /k  ( 

M_  =  M  cosy.  -  M.  siny 
*■1  r  it:  l 

x 

The  Mach  number  component  and  the  flow  deviation  in  the 
plane  normal  to  the  tangent  to  the  intersection  curve  are 

M..  =  <M-  2  +  Mr  2  l*5  ( 

Nl.  nl-  “It  .. 

1  1  -1  %i 

\  *  ri  '  tan  (  5^“*  < 

where  ,  which  represents  the  angle  between  the  direction 
of  the  free-stream  Mach  number  and  the  normal  (n)  to  the 
intersection  curve  at  each  point  i  in  this  plane,  is 
given  by 


(105) 


(106) 


(107) 


(108) 


T\  =  tan 


cos  (bz)  ^ 


(109) 


84 


**%> 


i 

1. 


r* 


w 


m 

* 


For  the  second  cone ,  a  similar  procedure  is  i_ ar x. red  out 
to  find  the  Mach  number  components  along  the  E,  E,  and 
n  directions,,  when  the  value  of  <0^,  corresponding  to  the 
generic  point  i  of  the  intersection  curve,  is  known. 
The  coordinates  of  the  generic  point  p.,  j/j.  in  the  x',y', 
directions  (Fig.  36)  are 

x '  .  =  p  .  sin>j) .  siny, 

111  1 

y1^  =  (p^  cosjji^  siny^-r)  cos«  +  (p^  cosy^-q)  sina 


(110) 


from  which 


<0.  =  tan 
1 


-1 


x '  . 
1 


(111) 


For  the  second  cone,  the  relations  corresponding  to 
equations  (105)  become  from  equations  (91),  (94),  and  (111) 
with  y  =  y 


Mr  =  M 

fc2.  x2 . 

1  1 


cos(£x).+  M  cos(^y).+  M  cos(i^z) 


J2. 

1 


M-  =  M  cos  (nx)  .  +  M  cos  (ny)  .  +  M  cos  (riz)  . 

n2.  X2.  1  y2.  1  Z2. 

IX  1  1 


“E. 


=  M  cos(bx).+  M  cos  (by).  +  M  cos  (^z)  . 


2 . 
1 


2. 

1 


(112) 


where 

v .  w. 

M  =  —  sin<0.  +  —  cos(0. 

x2  a2  1  a2  1 

i  i  i 

v,  w.  u. 

M  =( - cos  to.  +  —  sinw.)  cost*  -  —  sina 

y2.  a2.  1  a2.  1  a2. 

111  1 

V.  w.  u. 

M  =( — —  cos <0.  +  —  sin^o.)  sina  +  —  ccsa 

Z2.  a2.  1  a2.  1  a2. 

111  1 

The  Mach  number  component  and  the  flow  deviation  in  the 
plane  normal  to  the  tangent  to  the  intersection  curve  are 
represented,  similar  to  Eqs.(107)  and  (108)  by 


(113) 


2. 

1 


(M-  )  +  (Mr  )' 

2.  2 . 
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(114) 


Me  \ 


Ti  "  tan’  (  iT^  ) 

n2. 

.  i 


(115) 


The  oblique  shock  relations  in  the  plane 
normal  to  the  tangent  to  the  intersection  curve  are  used 
to  find  the  flow  conditions  at  each  point  i,  immediately 
behind  the  intersection  of  the  two  shocks,  which  due  to 
the  conical  nature  of  the  flow  field,  remain  valid  only 
in  the  immediate  vicinity  of  the  point.  In  this  present 
work,  the  flows  of  the  base  cone  and  the  second  cone  are 
considered  to  be  deflected  by  the  angles  62  and  6^ 

l  i 

respectively.  Using  the  oblique  shock  relations  with 

V 

and  62  /  the  Mach  number  behind  the  shock  Mg,  the  shock 

angle  A  and  the  pressure  and  speed  of  sound  ratios  P  /P-  * 
A2/A^  are  found  for  the  lower  portion  of  the  intersection. 

The  total  Mach  number  behind  the  shock  intersection  is 

M3  =  (M  2  +M  2)* 

1  t  .  1  . 

1  1 

where  M_  .  =  Mr  /(A./A.) 

3t±  tl.  2  1  1 


(116) 


(117) 


The  total  pressure  ratio  and  the  total  entropy  change  are 

P2  P1 
'  w  i  =  <  F7>x  r 

1  ( ®  r  _  2  .2  i 

-  ,  S  ^  I  ^vMn^ .  Sin  03i-(y-l)  j 

T  +  <n  . - T77 - J 


i  2y 


-  y  In 


y  +  1 


(y+ljM^sin^  «3 

1.  1 

1 


sm2  A_  +  2 
3. 

1 


(118) 


(119) 


where  p./p  and  S./k  are  the  values  immediately  behind  the 

JL  CO  X 


base  cone  shock. 


Using  and  6-^  ,  the  upper  portion  of  the 
2i  1 

flow  behind  the  shock  intersection  is  calculated  in  a  similar 
manner  from  the  oblique  shock  relations.  The  flow  properties 
immediately  behind  the  shock  are 


*v  =  Mt, /(Wi 


(120) 


M4.=  {M4.  +  %  5 

i  '-i 


P4i  P2  P' 
- —  =  (-2-  )  - 

P  P,  i  P 


(y+l)-(y-l)  (M*  2+  M*  cos»5^) 


(v+1) - (r-1) (mj  2} 


(121) 
y— 

y-l 

(122) 


e  -  —  a  cos-^ 


y-l  ,S1, 2  S2 

=  “  (-^—  +  «  —  COStfi.) 


4i  2? 


(123) 


1 

+  2y 


The  last  two  expressions  contain  factors  which  represent  the 
contribution  of  the  angle  of  attack  of  the  second  cone.  In 
the  present  analysis,  the  shock  strengths  have  been  assumed  to 
be  constant  through  the  intersections.  An  iterative  procedure 
can  be  used  to  improve  the  results.  The  iteration  consists  of 
varying  the  flow  deviations  6,  and  6o  until  the  pressures  in 

1  i 

the  lower  and  upper  region  behind  the  shock  intersection  are 
matched . 


f 

2y  mJ  sin2  fl  -(y-l)  ‘ 
2.  i 

i 

(y+ljM2  sin2 

2 .  i 

l 

|  In 

...  # 

rH 

+ 

1  ^ 

-y  -on 

(>-l)M^2  sin2  04  +  2 

2i  1 

i 

In  order  to  define  the  flow  field  behind  the 
shock  intersection,  the  characteristic  subroutine  can  be  used 
in  a  plane  which  contains  the  projection  of  the  intersection 
curve.  This  curve  becomes  the  data  line  and  the  flow  properties 
immediately  behind  it  are  obtained  by  averaging  the  values  in 
the  upper  and  lower  regions  and  finding  the  flow  deviation  in 
this  plane  from  the  projections  of  the  normal  and  tangent  Mach 
number  components. 

The  program  operation  and  example  problems  are 
given  in  Appendix  V  and  VI.  87 


E.  Impingement  of  a  Plane  Shock  On  a  Conical  Surface 

The  impingement  of  a  plane  shock  from  a  two- 
dimensional  wedge  of  angle  6^  on  a  conical  body  is  treated 
here.  In  order  to  treat  themost  general  case,  the  shock 
of  the  conical  body  is  included  ^Pig.  38).  The  intersec¬ 
tion  of  the  conical  and  plane  shocks  is  found  first.  After 
the  plane  shock  passes  through  the  conical  shock  it  is  no 
Iona'.  due  to  the  effects  of  the  conical  flow  field. 

The  ioc,v‘  !.  lination  of  this  shock  can  be  found  by  using 
varir  *  .ons  in  the  conical  field  between  the  conical 

shock  ‘  i  cone  surface  as  shown  by  dashed  lines  in 

Figure  •  In  this  analysis  the  computation  is  carried  out 
only  at  the  cone  shock  and  at  the  cone  surface.  The  proper¬ 
ties  obtained  by  crossing  the  "plane"  shock  at  the  cone  sur¬ 
face  are  used  in  the  manner  discussed  in  the  preceding  sections 
to  obtain  the  shock  reflected  from  the  cone  surface.  The  flow 
field  behind  the  reflected  shock  is  found  by  using  the  Charac¬ 
teristic  Subroutine  after  transforming  the  reflection  curve, 
which  lies  on  the  cone  body,  into  a  curve  on  a  plane. 

The  computation  starts  by  using  oblique  shock  theory 
for  the  given  wedge  angle,  6^,  and  free  stream  Mach  number  to 
obtain 

Ml'  0r  Pj/P.  and  ai/aw 


which  represent  the  uniform  flow  properties  behind  the  plane 

shock.  As  in  the  previous  sections,  the  conical  field  of  the 

cone  at  zero-angle  of  attack  is  defined  by  imputing  the  para-, 

meters,  y  ,  M*  ,  j/)1  ,  P  /p  T  /T  and  S„/R  behind  the  conical 
C  1  1  C®,Ce°  2 

C  C 

shock  and  y^ ,  "It-  *ib\  p./p^  on  the  body  (see  Ref.  15).  This 
data  allows  the  “conical  flow  properties  to  be  computed.  The 
Mach  number  components  along  and  normal  to  the  conical  ray 
y  =  const. are 


M  =  u/a 
r 

Mt  =  Mr  tan  ^i”y) 


(124) 
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The  intersection  of  the  conical  shock  and  the 
the  z'-x'  plane,  (Pig.  38),  is  given  by 


2Bz + 
1 


where 


A  =  sin29^  -  b2  cos2fi^ 

B  =  r  sinfl^  +  b2q  cosQ^^ 

C  =  r2  -  b2q2 
b2=  tan2yc 

and  where  z is  obtained  by  incrementing  from 
corresponding  to  x'  =0  (z'q  )  and  is  given  by 

z  =  z  ’  +  i  A  z  1 
i  o 


(125) 

plane  shock  on 

(126) 

(127) 

the  value  of  z ' 

(128) 
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The  angle  between  the  tangent  to  the  intersection  curve  and 
the  z'-axis  is 


B-Az  1  . 


1 


i 


The  computationof  the  intersection  is  terminated  when 

Q.  <5°  in  order  to  work  on  one  side  of  the  cone.  The  conical 
i 

flow  Mach  number  components  immediately  behind  the  shock, 
(y-v  ) ,  in  the  x'-y'-z'  directions  are 


+ 


M 


2 


sinvc  sirufi^  +  Mfc  cosyc  sinw^ 

M^  (sinyc  costa^  cos0^  +  cosyc  sin0^ ) 

M,_  (cosy^  cos^  cosP^  -  siny^  sin^) 

M  (cosy  cos0_  -  siny  cos®,  sin©..) 

r  C  X  C  X  X 

M  (cosy  cos  sinV  +  siny  cosfl  ) 

u  O  -L  X  C  X 


where  the  meridian  angle 


=  tan 


is 


with 


x .  =  x 

x  x 


y.  =  q  +  z sin0. 
x  x  1 


Equations  (13$  projected  along  the  normal  and  the  tangent 
to  the  intersection  curve  with  z',  x'  plane  become 


sinO. 

x 


M  cosQ. 

^  i  X 

x  . 

1 


M 


2 


t. 

l 


-M„  cosQ. 

2z'. 

1 


sinO. 

l 


(129) 


(130) 


(131) 


(132) 


(133) 
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With  the  total  Mach  number  component  in  the  generic  plane 
normal  to  the  intersection  curve  at  point  i  given  by 


M2  =  (M2  ^  +  Mo  2  ^ 

\  y'i  n. 

1  2-1 


(134) 


and  the  strength  of  the  plane  shock  given  by, 


tan9  .i.  tanf^-6^ 

6i*  ■ tan  1  -tan  Tirol — 


(135) 


the  oblique  shock  relations  can  be  used  to  find  the  flow 
properties  behind  the  lower  region  of  the  shock  intersection. 
The  properties  thus  found  are  the  Mach  number  component 
M  ,  the  shock  angle,  A  ,  in  the  plane  normal  to  the  inter - 
4n.  l 

section  curve  and  the  pressure  and  speed  of  sound  ratios 
(P4/P2)^  and  (a4/a2)i,  respectively.  The  tangent  Mach  number 
component  ,  the  total  Mach  number  M4  ,  the  pressure  ratio 

t .  i 

l 

(p^/P^)^*  the  speed  of  sound  ratio  (a4/a^)^  and  the  change  of 
entropy  are  respectively 


V  ■  m2  '  r  »i 

t.  t.  4 

l  i 


\  =  V  *  +  (\f 


<VpJi  =  <Wi  (pc/p.  1 
<vaJi  -  <a4/a2>i  /tX 


(136) 


r  2vM0  1  Sin  A  -  (y-1) 
!  .  8,  |  \  4i 

s4i=  3?  H - Thl - 

(y+l)  M22  Sin2  A4  "1  ’J 


-v  tn 


(v-1)  M  2  sin2  A  +  2 
N.  -  i 

i 


The  flow  deviation  with  respect  to  the  z'  axis  is 

m4 

-1 

Oa.  =  -  tan  pr - 


a4.  =  Qi 
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(137) 


For  the  eventual 
it  will  be  useful 
which  is 


62*  =  tan 
i 


computation  of  the  upper  intersection  region 
to  know  the  strength  Of  the  conical  shock, 


M. 


-1 


'M, 


tan 


-1 


tanA 

(HTn.> 


(138) 


n . 

i 


In  each  plane  normal  to  rhe  intersection  curve  a  new 
orthogonal  coordinate  system,  (£,rj),  is  defined,  in  which  £  is 
in  the  x'-z’  plane  and  is  coincident  with  the  normal  direction 
previously  defined,  (Fig.  39).  The  angle  of  the  lower  shock, 
referred  to  the  f  axis,  is 

m2 

®  4 =  e4.  -  tan"1  (mT^)  (139). 


For  each  normal  plane  the  point  *i^,  of  intersection 

between  the  shock  and  the  conical  body  is  found  from  the 
following; 


A2_  n 

VP;  P 


(140) 


=  -r?.  ootan 
1  1 

where 

P  =  L  +  ex  (E  eL  +  2  G) 

Q  =  e1  F  +  H 
and  where 

2  2  2  2  2 
E  =  cos  Oi  +  sin  Q^(sin  0^-b  cos  4^) 

2 

F  =  -x.'  cos  0.  -  (F  sin©  +b  F  cos0  )  sinO. 

1  X  X  X  ^  X  X 

2 

G  =  -sin0^  cosA^  sinJh  (1+b  ) 

2 

H  =  F^  cos0^-b  F2  sin0^ 

2  2  2 
L  =  cos  0^-b  sin  0^ 


(141) 


(142) 
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N  =  (x..1  )2  +  F12-b2  F22 

F  =  r-  z!  sin0, 
111 

F  =  q  +  z cos 6, 

2  1  1 

e, =  -cotan  I  A '  | 


(143) 


The  coordinates  of  these  intersection  points  on  the  body  are 
transformed  into  a  conical  coordinate  system,  p,  tr>,  which 
represents  the  polar  distance  from  the  conical  body  apex  and 
the  meridian  angle  in  the  z  =  constant  plane,  respectively, 
(see  Fig.  40) .  This  transformation  is  given  by 


j.  =  +  { x .  2  +  y .  2  +  z .  2} 

i  l  i,  i,J 


<o.  =  tan 


V\ 


(144) 


where 


x.  =  x.'  -  £.  cos  Q. 

n  i  -i 


y.  =  r  +  r».  cos  0.  -  sin  0  (z.*  +  .£.  sin  a) 

X  i  1  • 

z.  =  q  +  n.  sin  0,  +  cos  0,  (z.'  +  f.  .  . 

i,  n  i,  1  1  l  si,  sm  a. 


(145) 


In  order  to  analytically  represent  the  intersection 
of  the  shock  on  the  conical  body,  a  second  order  curve  of  the 
form 


ak  <p2  +  bk  p2  +  dk  p  +  1  -  0 

is  fitted  to  three  selected  points 


(146) 


f>o<V  V  ®1!  °2 '  ®2 

which  are  the  center  line  point,  a  point  in  the  region  near 
the  center  line  and  one  of  the  last  points,  respectively.  The 
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coefficients  of  Eq.  (146)  are 


*k  *  (prpo)  (p2-po)/Ak 

bk  =  [<012  (p2“po)“<022  (prpo)]/Ak 

dk  ~  ^22  (pl2-po  ^1  ^p2  -po2)^/Ak 

where  Ak  =  PQ[o12  P2  fo2-DQ)-  *2  P1(o1-PQ)  3 


Similar  to  Section  D  a  set  of  coordinates  tancjent,  binormal 
and  normal  to  the  curve  represented  by  Eq.  (146)  are  defined. 
Their  direction  cosines  with  respect  to  the  x-y-z  axes  are 


cos(i^x).  =  (<«.  sin  >n.  sin  y,  -  AA.  cos  ro.)/BB. 

l  'b  l  x  x 

cos  (t^y)  t  a  cos  >oi  sin  yfa  +  AAi  sin  (148) 

cosrfz^  =  (<fii  cos  yb)/BBi 


and 


>(&). 


.  =  sm  (Q.  cos  y, 
l  l  b 


cos  rfyh  =  COS  (p.  cos 


cc  .  :£z)  .  =  -  sin  y, 
i  b 


cos (nx) .  =  NA. 

i  i 

cos (ny) .  =  NB . 
J  x  1 

/  A  . 

cos(nz)^  =  AA^ 


(149) 


(150) 
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where 


Mi  =  <bk  pi  +  T 


f  2  2  a  *5 

BB.  =  -I  to.  +  AA.  }> 
l  iTl  l  J 


p .  sm  y 
Hi  'b 


“k 


AA.^  =  cos(£x)^  cos  (lay)  ^  -  cos(^y)i  cos  (lax)  ^ 
NA^  =  cos(^y)i  cos(^z)i  -  cos^z).^  cos  (lay)  ^ 

i 


NB.  =  cos(^z).  cos(bx).  -  cos  (l^x)  .  cos  (laz) 


Next  it  is  desired  to  find  the  flow  field 
behind  the  "plane"  shock  near  the  conical  surface.  The  Mach 
number  components  of  the  conical  flow  field  in  front  of 
the  "plane"  shock  on  the  conical  body  are 
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where  is  found  from  Eqs. (124)and  (125)with  y  =  v^, 

and  i ||).  =  0.  .  The  components  normal  and  tangent  to  the 
1  ib  . 

original  intersection  curve  are 
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(153) 


The  flow  field  behind  the  shock  near  the  conical  body  can 
be  found  by  using 


(154) 


and  6  *  (Eq.  135  ),  in  the  oblique  shock  relations.  The 
i 

shock  strength  has  been  assumed  to  be  constant;  however, 

the  flow  behind  the  shock  near  the  body  differs  from  that 

of  the  lower  region  of  the  shock  intersection  since  the 

flow  field  in  front  of  the  shock  is  conical.  In  the 

usual  manner  the  properties  obtained  are:  the  normal  Mach 

number  M.  ,  the  shock  angle  0.  ,  the  pressure  and  speed 
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number  and  the  total  change  in  entropy  are 
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The  Mach. number  components  normal  and  binormal  to  the 
original  intersection  curve  are 
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To  obtain  the  flow  field  behind  the  reflected 
shock  it  is  necessary  to  find  the  upstream  Mach  number  com¬ 
ponents  in  the  directions  tangent,  normal,  and  binormal  to 
the  reflection  curve  which  are 
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Then  the  total  Mach  number  component  in  a  plane  normal  to  the 
reflection  curve 
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and  the  flow  deflection  across  the  reflected  shock  in  this 
plane 

6  =  tan 

4 . 

1 

can  be  used  in  the  oblique  shock  relations  to  give  the  flow 
properties  behind  the  reflected  shock  on  the  cone  surface. 
The  results  are  finally  given  by  the  following  relations: 
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where 


is 


obtained  from  the  oblique  shock  relations  and 
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The  shock  angle  with  respect  to  the  body  is 


The  pressure  ratio  is 
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The  flow  deviation  in  the  p ,  <p  plane  with  respect  to  the 
p  axes  is 
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and  the  tofal  <~>v>  Anna  l  r>  nr*4*v nr\» *  4  _ 
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To  find  the  flow  field  on  the  body,  the  characteristic  sub¬ 
routine  can  be  used  in  the  p,  <p  plane,  with  the  above  condi¬ 
tions  which  are  valid  immediately  behind  the  reflected  shock. 


The  program  operation  and  example  problems 
are  given  in  Appendices  v  and  VI. 
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APPENDIX  I 


HOW  TO  USE  THE  PROGRAM 

(Two-Dimensional  and  Axisymmetric  Flow  Computer  Program) 

1.  General  Comments 

Very  often  instructions  for  operating  a  computer  program  are 
more  complicated  than  even  the  program  itself.  Many  quite  useful 
programs  have  been  rendered  almost  useless  because  the  answers  to 
prospective  users  questions  are  either  not  provided  by  the  in¬ 
struction  manual  or  are  hidden  in  a  plethora  of  jargon.  The  im¬ 
portant  questions  which  are  answered  here  are: 

a.  What  does  the  program  do? 

b.  What  are  its  restrictions  and  limitations? 

c.  What  rnput  data  are  required  by  physics?' 

d.  What  additional  input  data  are  required  to  make  the 
program  work? 

e.  What  is  the  required  form  of  the  input  data? 

f.  What  is  given  in  return  (i.e.  what  are  the  outputs)? 

g.  How  is  the  output  to  be  interpreted? 

h.  What  might  go  wrong?  What  can  be  done  about  it? 

Answers  and  parts  of  answers  to  some  of  these  questions  have 
already  been  given.  In  such  instances,  the  appropriate  sections 
will  be  designated  and  that  information  will  not  be  repeated  in 
this  section.  The  questions  and  answers  are: 

a.  What  does  the  program  do? 

The  program  computes  the  inviscid  supersonic  internal  flow 
of  a  perfect  gas  of  constant  specific  heat  ratio  for  configurations 
typical  of  supersonic  and  hypersonic  inlets.  (Also  see  Section  IV-A,B) . 
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b.  What  are  its  restrictions  and  limitations? 

In  standard  notation: 

y  =  constant  (constant  ratio  of  specific  heats) 
p  =  pRT  (perfect  gas  equation  of  state) 
ju  =  0  (inviscid) 

M  >  1  (supersonic) 

c.  What  input  data  are  required  by  physics? 

Referring  to  Figure  11,  the  physics  of  these  flows  requires 
one  to  specify  the  flow  conditions  completely  on  some  initial  data  line 
(which  is  not  characteristic)  (AA1  of  Figure  11)  and  the  geometry 
(y(x),  y 1 (x) )  of  both  the  walls.  When  this  data  has  been  specified 
(maximum  ordinate  of  the  initial  data  must  be  great  enough  so  that  a 
down  characteristic  from  that  point  intersects  the  leading  edge  of  the 
upper  wall)  the  entire  inviscid  flow  in  the  inlet  can  be  determined,  so 
long  as  it  remains  everywhere  supersonic. 

Needless  to  say,  one  must  also  specify  y,  the  ratio  of  specific 

heats . 


d.  What  additional  input  data  is  required  to  make  the  program  work? 

In  addition  to  the  physically  required  data,  the  program  needs 
the  following  information:  abscissa  of  upper  and  lower  leading  edges, 
maximum  abscissa  of  the  upper  wall  to  be  computed,  angle  of  attack  of 
the  free  stream  flow  with  respect  to  the  horizontal  axis,  spacing  of  mesh 
points  on  the  initial  data  line  (uniform  incident  flow) ,  a  code  indicating 
whether  data  along  the  initial  line  i"  supplied  or  the  simplified  input  is 
used,  a  code  indicating  whether  the  flow  is  two-dimensional  or  axi- 
symmetrical,  and  other  codes  related  to  the  wall  geometry.  In  addition, 
for  the  convenience  of  the  user,  provision  is  made  for  inputing  and 
outputing  the  date  and  a  run  indentification  number. 

e.  What  is  the  required  input  data? 

Physical  Data 

As  mentioned  before  all  input  data  must  be  nondimensionalized 
in  the  manner  specified  in  Part  C-a,  Non  Dimensionalization  and  General 
Equations  of  Motion.  Read  Part  C-a  before  continuing  with  this  section. 
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If  the  incident  flow  is  uniform,  the  user  may  specify 
the  incident  Mach  number,  M,  the  specific  heat  ratio,  y,  and,  in 
the  case  of  two-dimensional  flow,  the  angle  of  attack.  In  this 
case  the  reference  pressure  and  temperature  will  be  assumed  to 
be  the  incident  pressure  and  temperature. 

In  the  case  of  uniform  incident  flow,  the  leading  edge 
of  the  lower  wall  must  be  upstream  of  the  leading  edge  of  the 
upper  wall.  However,  if  it  is  not,  this  requirement  can  be  over¬ 
come  by  extending  the  lower  wall  upstream  of  its  leading  edge 
parallel  to  the  incident  flow.  This  can  not  be  done  in  the  case 
of  axisymmetric  flow  when  the  leading  edge  of  the  lower  wall  has 
a  sharp  corner  on  the  axis  of  symmetry.  If  it  is  desired  to  have 
the  output  referred  to  another  pressure  and  temperature,  it  is 
necessary  to  input  data  on  the  initial  data  line  point  by  point, 
the  data  being  referred  to  the  desired  reference  pressure  and 
temperature. 

If  the  incident  flow  is  non-uniform,  it  is  necessary  to 
specify  the  data  point  by  point,  starting  at  the  leading  edge  of 
the  lower  wall.  A  reference  pressure  and  temperature  (and  length, 
if  desired)  must  be  chosen.  Then  the  coordinates  of  each  data 
point  must  be  specified  along  with  the  Mach  number,  the  angle  the 
streamline  makes  wi*-1-'  the  x  axis,  and  either  the  static  pressure 
and  temperature  or  the  total  enthalpy  and  total  pressure  (all 
non-dimensionalize J) . 

It  is  important  to  realize  that  the  initial  data  line  may 
not  be  arbitrarily  specified.  The  first  point  must  be  at  the  lead¬ 
ing  edge  of  the  lower  wall  and  the  subsequent  points  must  be  in 
order  of  increasing  ordinate.  Also,  the  initial  data  line  must 
be  oriented  so  that  the  characteristics  emanating  there-from  point 
downstream  from  the  data  line  (see  Figure  26) .  In  addition,  the 
inital  data  line  must  be  upstream  of  the  leading  edge  of  the  upper 
wall  (cowl) . 

The  required  data  are: 

1)  Free  stream  Mach  number,  M  (simplified  input,  uniform 
incident  flow) or  x,y, M,  0,  either  p  and  T  or  pfc  and  H  on  the 
initial  data. 

2)  Specific  heat  ratio,  y. 

3)  Abscissae  of  the  two  leading  edges  and  the  maximum 
abscissa  to  be  computed  on  the  upper  wall.  All  lengths  must  be 
dimensionally  the  same  or  non-dimensionalized  consistently.  If 
non-dimensionalized,  the  program  does  not  need  to  know  how  this 
was  done,  nor  is  there  any  way  of  telling  it. 
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4)  Angle-of-attack  in  degrees  measured  positive  counter¬ 
clockwise  from  the  horizontal  axis.  (Two-dimensional,  uniform  incident 
flow,  only.) 

5)  The  code  for  two-dimensional/axisymmetrical  flow  is 
an  integer  which  is  zero  (0)  if  the  flow  is  two-dimensional  and  one 
(1)  if  it  is  axisymmetrical. 

6)  The  run  identification  number  and  the  date  are  read  in 
as  integers,  the  identification  number  having  a  maximum  of  5  digits. 

7)  Miscellaneous  control  integers. 

8)  y  (x)  on  the  lower  and  upper  walls  in  the  form  of 
second  order  polynomials.  There  are  a  maximum  of  9  second  order 
polynomials  on  each  wall.  The  polynomials  are  of  the  form 

2 

y(x)  =  a.  +  b.  (x-x.)  +  c.  (x-x.- )  ,  1=1,  9 

1  1  A  i  1 

for  top  and  bottom,  a.,  b.,  c.,  and  x.  must  be  specified  for  each 
polynomial.  i 2 3 4 5  “ 

Note:  -  A  more  complicated,  but  extremely  versatile 
way  of  specifying  body  ge'  '’etries  is  given  in  Appendix  A.  It  can 
be  used  to  specify  body  c  ometries  in  virtually  any  manner. 

9)  The  x  coo.  inate  of  every  sharp  corner  on  the  upper 
and  lower  wall,  but  not  ^  eluding  the  leading  edges.  In  addition, 
the  user  must  supply  dy/dx  before  and  after  che  corner.  There  are 
a  maximum  of  six  such  corners  on  each  wall. 

f.  What  is  given  in  return?  (What  are  the  output?) 

Data  at  each  mesh  point,  including  the  location  of  the  point, 
are  written  on  scratch  tape  in  the  following  form: 

1)  abscissa  and  ordinate  of  the  point  (in  the  same  units 
or  non-dimensionalized  in  the  same  manner  as  the  input  lengths) . 

2)  pressure  divided  by  the  reference  pressure 

3)  temperature  divided  by  the  reference  temperature 

4)  the  streamline  angle,  0 

5)  Mach  number 


107 


6)  (S-SqVv]v  where  S  is  dimensional  entropy  (SQ  is  the 
entropy  corresponding  to  the  place  where  p^  and  Tq  are) . 

7)  total  enthalpy  divided  by  the  ratio  of  reference 
pressure  to  reference  density. 

8)  the  mass  flow  between  the  mesh  point  and  the  lower 
wall  divided  by  [square  root  of  (reference  pressure  x  reference 
density)  times  the  square  of  the  length  used  to  non-dimensionalize 
lengths  (if  they  are  non-dimensionalized) ] . 

.  mass  flow  p  =  o  R  T 

i.e.  T==?t  .  v  o  Ko  o 

*  PQ  Po  4 

In  addition,  in  axisymmetric  flow  the  mass  flow  is 
also  divided  by  n- 

9)  Pitot  pressure  divided  by  the  reference  pressure, 
g.  How  is  the  output  to  be  interpreted? 

The  output  is  given  as  a  function  of  mesh  points  on  each 
down  running  characteristic.  In  each  strip  the  down  running  character¬ 
istics  are  identified  integrally  beginning  with  1  for  the  first  point 
to  be  computed  (see  Part  D-b  on  the  order  of  progression  of  the  com¬ 
putation)  in  that  strip.  Similarly,  the  strips  are  identified  integ¬ 
rally  beginning  with  1  for  the  strip  containing  the  initial  data  line. 

These  should  be  no  problem  so  far  as  reference  quantities 
are  concerned. 

There  is  a  little  problem  arising  from  the  way  in  which 
the  computation  proceeds.  It  will  be  recalled  that  (see  Part  D-b-2,3) 
when  a  down  shock  is  initiated,  the  region  between  (see  Figure  16) 
the  last  down  characteristic  BC,  in  the  previous  strip  (Ith  strip)  and 
the  down  shock,  BD,  separating  the  two  strips  (i.e.,  the  region  between 
BC  and  BD)  is  computed  first  as  a  part  of  the  Ith  strip  without  con¬ 
sidering  the  influence  of  the  compression  at  B,  and  then  is  recomputed 
as  a  part  of  the  (I+l)st  strip  taking  into  account  the  compression. 
This  is  necessary  since  we  have  no  way  of  knowing  a  priori  just  where 
the  shock,  BD,  will  go.  We  only  know  that  it  will  always  be  upstream 
of  BC. 
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For  this  reason,  in  computing  the  Ith  strip,  we  write  out 
all  of  the  data  in  that  strip  up  to  and  including  all  of  the  data  on 
BC.  Then,  in  recomputing  the  region  common  to  both  strips,  we  out¬ 
put  the  new  "good"  data  between  BD  and  BC.  Thus,  one  must  be  careful 
to  pick  the  "good"  data  for  the  common  region. 

But  this  is  made  simple  by  the  following  device. 

Recall  that  the  various  strips  are  separated  by  down 
running  shocks.  In  general,  the  Itn  down  shock  will  separate 
the  Ith  and  the  (I+l)st  strip.  Neglecting,  for  the  moment,  the 
effect  of  any  other  down  shocks,  the  data  ambiguity  is  a  result  of 
the  fact  that  some  of  the  data  computed  in  the  Ith  strip  is  not  valid 
because  it  is  downstream  of  the  Ith  shock.  So  the  proper  procedure 
is  to  plot  first  the  down  running  shocks  and  to  use  their  locations 
as  a  guide  in  determining  which  data  has  become  invalid. 

Note  that,  if  two  down  running  shocks  coalesce,  being 
continued  as  one  shock,  downstream  of  the  coalescence  the  difference 
between  the  strip  number  on  either  side  of  the  shock  will  be  greater 
than  one  (Figure  25)  . 

h.  How  long  will  a  computation  take? 

The  program  is  written  in  Fortran  IV  for  use  with  IBSYS  13 
on  the  IBM  7094.  A  typical  case  will  take  about  1-4  minutes  execution 
time. 


On  the  CDC  6600,  the  same  computation  would  be  completed,  in 
less  than  2  minutes  including  compilation  time. 

i.  What  might  go  wrong  and  what  can  be  done  about  it? 

Any  computer  program  developed  to  compute  phenomena  as  com¬ 
plicated  as  those  under  consideration  here  is  bound  to  have  trouble 
from  time  to  time.  Prospective  users  should  be  neither  surprised  nor 
frightened  by  this  fact.  A  good  program  will  have  difficulties,  but 
they  should  not  occur  too  frequently  under  normal  use  nor  should  they 
occur  when  relatively  "simple"  cases  are  being  computed.  But,  when 
troubles  arise,  the  user  should  be  given  every  opportunity  to  get 
around  the  difficulties  in  order  to  obr.ain  the  information  (or  at  least 
some  of  it)  desired. 


The  difficulties  of  this  flow  field  program  will  ordinarily 
fit  into  two  categories,  those  resulting  from  computational  inaccuracies 
and  those  resulting  from  logical  complications  caused  by  extremely  com¬ 
plicated  flow  field  patterns. 

1)  Problems  resulting  from  computational  inaccuracies  - 
The  usual  reason  for  problems  to  arise  because  of  computation  in¬ 
accuracies  is  that  dependent  variables  have  large  variations  on  the 
scale  of  a  mesh  size  or  less.  For  instance,  if  one  has  an  expansion 
corner,  the  pressure,  density  and  streamline  orientation  are  dis¬ 
continuous  at  the  corner.  Hence,  we  make  provisions  to  compute  the 
flow  near  the  corner  in  special  ways  in  order  to  retain  the  desired 
accuracy.  If  we  did  not  make  such  provision,  we  not  only  would  lose 
accuracy,  but  we  would  usually  find  that  the  errors  made  would  be  so 
large  that  the  computation  would  soon  fail  to  converge  in  one  of  the 
many  subsequent  iterations. 

Suppose  that  in  some  region  one  of  the  boundaries  changes 
so  rapidly,  relative  to  the  local  mesh  scale,  that  there  are  large 
changes  in  the  pressure,  density, and  streamline  orientation  during  one 
mesh  interval  on  the  body.  Then,  even  though  things  are  smooth,  because 
of  the  rapid  changes  it  is  not  unlikely  that,  in  addition  to  inaccuracies, 
there  will  be  difficulties  and/or  failures  in  attempting  to  converge  in 
some  iteration. 


One  obvious  solution  is  to  decrease  the  mesh  size  by  de¬ 
creasing  the  initial  mesh  size.  This  can  be  done  within  limits  (the 
computer  has  limited  storage  capacity) ,  though  a  price  is  paid  in  in¬ 
creased  computer  time.  If  it  is  impossible  to  refine  the  mesh  enough, 
an  alternative  is  to  replace  the  rapid  smooth  turning  with  a  sharp 
corner  strategically  placed  and  oriented.  Then,  locally,  near  the 
corner,  the  solution  will  not  describe  well  the  flow  over  the  original 
body,  but  the  engineer  can  usually  correct  the  computed  flow  to  give 
him  the  required  information. 

In  certain  cases  the  user  is  justified  in  relaxing  certain 
tolerances  in  order  to  allow  an  iteration  to  converge.  Such  action  must 
be  done  with  an  appreciation  for  the  effect  chat  this  will  have  on  the 
computation.  Some  tolerances  can  be  relaxed  by  an  order  of  magnitude 
under  certain  conditions  without  having  a  marked  effect  on  the  accuracy 
of  the  computation.  Others  should  not  be  relaxed  by  even  a  factor  of 
two.  Such  decisions  require  a  quite  intimate  knowledge  of  the  program 
and  should  not  usually  be  attempted  by  infrequent  users. 
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?)  T.oaical  nroblems  resnlhina  from  exf-remfilv  romnl i rated 

•  —  -t  x  _/  4  a. 

flow  field  patterns  -  For  such  problems  all  of  the  prescriptions 
applicable  to  inaccuracy  problems  also  apply,  in  addition  to  a  few 
others.  The  problems  we  are  now  discussing  result  from  very  complicated 
flow  field  patterns.  For  example,  if  three  down  running  shocks  con¬ 
verge  near  the  lower  wall  (perhaps  they  converge  just  before,  or  the 
reflected  shocks  converge  just  after,  or  some  combination  happens)  near 
a  sharp  corner,  it  would  not  be  too  surprising  to  find  that  the  program 
gets  a  little  confused  trying  to  keep  track  of  what  happened,  what 
is  happening,  and  what  might  happen.  Sometimes  the  program  will  fail 
to  work  perfectly  for  such  reasons.  This  does  not  mean  that  iu  is 
useless  or  that  the  user  is  lost.  There  are  many  alternatives  avail¬ 
able,  including  those  already  discussed  under  a) ,  above. 

Another  "corrective"  measure  is  best  explained  by  an 
illustrative  example.  Suppose  that  for  a  given  geometry  and  Mach 
number,  such  a  complicated  pattern  develops  in  a  particular  region  and 
the  program  gets  confused  and  stops.  Very  often  a  change  in  Mach 
number  (in  both  directions)  will  spread  things  out  so  that  they  are 
not  quite  so  concentrated.  One  can  use  this  to  infer  his  desired 
results . 

j.  Special  control  cards;  Tape  Assignments 

The  program  is  written  in  FORTRAN  IV  for  use  on  the  IBM 
7094  or  CDC  6600  computer. 

Approximately  72,000  octal  locations  are  needed  on  the 
CDC  6600.  No  special  control  cards  are  needed  ether  than  specifi¬ 
cations  of  the  use  of  logical  tape  5  for  input,  6  for  output,  and 
3  and  4  as  scratch  tapes. 

On  the  IBM  7094,  the  same  tape  assignments  are  used.  Be¬ 
cause  of  storage  requirements,  it  is  necessary  to  use  the  alternate 
input-  output  option,  the  ALTIO  routines.  This  is  dome  by  specifying 
ALTIO  as  one  of  the  options  on  the  $IBJ0B  control  card. 


Ill 


■Wi 


Input 


£.  . 


Input  are  read  on  logical  tape  5  and  the  results  written  on 
logical  tape  6.  In  addition,  the  results  are  stored  on  logical 
tape  4  for  use  during  the  computation.  But  that  tape  can  be  used 
later  also.  Logical  tape  3  is  also  used  as  a  scratch  tape. 

The  data  for  several  cases  can  be  stacked. 


INPUT  DATA 


All  Integers  Must  Be  Right  Adjusted 


Card  Columns 
No. 


Datvm- 
Fortr an 
Symbols 


Format 


Must  be 

Right 

Adjusted 


Recom¬ 

mended 

Value 


Description 


6-10 

11-15 

16-20 


21-25 


26-30 


31-35 


36-40 


NRUN 


MONTH 

MDAY 


MYEAR 


IUNFLO 


ILWFTS 


IUPFTS 


Run  number 

Number  of  the  month 

Number  of  the  day 

Last  two  digits 
of  the  year 

Must  be  0-If  initial  data  is 

sO  uniform  and  to  be 

computed  from 
given  Mach  number 
and  y.  In  this 
case  the  initial 
data  line  is  the 
line  x=XBEG  (see 
card  #2).  If  the 
initial  data  is  to 
be  read  for  each 
mesh  point,  then 

IUNFLO  equals  the 
number  of  such 
points. 

Must  be  0  for  two-dimen- 

0  or  1  sional  flow,  1  for 
axisymmetric  flow. 

<;9  Number  of  quadratic 

fits  on  the  lower 
wall 

[y  (x)=ai+bi(x-xi)  + 
c. (x-x. ) 2] 

i  l 

<.3  Number  of  quad¬ 

ratic  fits  on  the 
upper  wall 
[y(x)=ai+bi(x-xi)  + 

c .  (x-x . ) ^ ] 

i  l  J 
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Card 

Columns  Datum- 

Format  Must  be 

Recom¬ 

Description 

No. 

Fortran 

Right 

mended 

Symbols 

Adjusted 

Value 

1  41-45  INPT  15  Yes 


46-50  ICORNL 

51-55  IC5RNU 

±  i 

2  1-10  MACH  E10.8  Optional 

11-20  GAMMA 

21-30  ALPHA 

31-40  DY 

41-50  XAU 

"1-60  XEND 


61-70  XBEG 


+1  INPT=+1  if  IUNFLO  /  0 

and  p  a:.d  T  are  to  be 
specified . 

INPT=-1  if  IUNFLO 
and  pT  and  H  are  to 
be  specified. 

s6  Number  of  corners  on 

lower  wall  (not  count¬ 
ing  lower  leading 
edge) 

6  Number  of  corners  on 

upper  wall  (not  count¬ 
ing  upper  leading  edge 

y=c  /c  ,  the  ratio  of 
p  v 

specific  heats.  Angle 
of  attack  of  free 
stream  with  respect 
to  x  axis,  positive 
counter-clockwise 

Ay  interval  for 
initial  line  x  co¬ 
ordinate  of  the  lead¬ 
ing  edge  of  the  upper 
wall 

Maximum  abscissa  of 
the  upper  wall.  The 
computation  will  be 
stopped  when  this 
value  is  reached  and 
the  next  case  will  be 
begun. 

Abscissa  of  the  lead¬ 
ing  edge  of  the  lower 
wall . 


Card  Columns 

No. 


Datum- 

Fortran 

O  *  »  wVi  r\ "!  c* 
Wjr  « 4. 


Format 


Must  be 
Right 

7k/^  *i  n  r<  r»  rl 
nu  J  U O 


Recom¬ 

mended 

TT*i  1  ii/i 
V  C4.U.  VX  V5 


Description 


3  1-10 

11-20 
21-30 

51-60 


XFLAGL(l)  E10.8  Optional 

(2) 

(3) 

(6) 


1-10 

TFLAGL  (1,1) 

11-20 

(1,2) 

51-60 

(1.6) 

1-10 

11-20 

« 

TFLAGL  (2,]) 
(2,2) 

« 

«< 

51-60 

(2,6) 

6  1-10 
11-20 


51-60 


XFLAGTJ { 1 ) 

(2) 

(6) 


*  If  there  are  no  corner s,  insert  blank  data  cards. 


x  coordinate  of 
the  sharp  corners 
on  the  lower  wall . 

There  are  a  maxi¬ 
mum  of  6  such 
corners.  They 
must  be  in  order 
beginning  at  the 
first  sharp  corner 
after  the  leading 
edge.  There  are 
ICORNIi  corners 
(See  card  #1)  * 

T=dy/dx  on  the 
upstream  edge  of 
the  sharp  corners 
on  the  lower  wall, 
in  order .  * 

T=dy/dx  on  the 
downstream  side 
of  the  sharp 
corners  on  the 
lower  wall,  in 
order,  .* 

x  coordinate  of 
the  sharp  corners 
on  the  upper  wall 
(maximum  of  6)  . 
They  must  be  in 
order  beginning 
with  the  1st  corner 
after  the  leading 
upper  edge.  There 
are  ICORNU  corners 
(See  card  #  1)  * 
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Card  Columns  Datum-  Format  Must  be  Recom-  Description 

No.  Fortran  Right  mended 

Symbols  Adjusted  Value 


7  1-10  TFLAGU  (1,1)  E10.8  Optional 

11-20  (1,2)  I  I 


51-60  (1,6) 

8  1-10  -TFLAGU  (2,1) 

11-70  (7,2) 


51-6C  (2,6, 

9  1-10  YLBC (1,1) 

11-20  (1,2) 

21-30  (i/3) 

31-40  (1/4) 


T=dy/dx  on  the  up¬ 
stream  side  of  the 
sharp  corners  on 
the  upper  wall,  in 
order.  * 


T=dy/dx  on  the 
downstream  side  of 
the  sharp  corners 
on  the  upper  wall, 
in  order.  * 


Coefficients  in  the 
first  of  the  second 
order  polynomial 
fits  for  y(x)  on 
the  lower  wall. 
y(x)=a  +b  (x-x1)  + 

c1(x-x1)2  x1s:xsx2 

where  a^  is  in 

columns  1-10,  b^  is 

in  columns  11-20, 
c^  is  in  columns 

21-30,  Xj^  is  in 

columns  31-40. 

There  are  a  maximum 
of  9  such  fits. 

The  actual  number 
of  such  fits  is 
ILWFTS  which  is  a 
datum  on  the  first 
card . 
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Cord  Columns  Datum--  Fux. maL  Must  be  Recum- 

No.  Fortran  Right'  mended 

Symbols  Adjusted  </alue 


Description 


10  1-10 
11-20 
21-30 
31-40 


YLBC  (2,1)  E10.8  Optional 

(2,2) 

(2.3) 

(2.4) 


9+ILWFTS-l  1-10 
(ILWFTSS9) 


9+ILWFTS  1-10 
11-20 
21-30 
31-40 


YI£C  ((ILWFTS,!) 


YDBC  (ILWFTS,  4) 
YUBC (1,1) 

(1,2) 

(1.3) 

(1.4) 


9+ILWFTS+  1-10 
IUPFTS-1  11-2- 
21-30 
31-40 


YUBC (IUPFTS , 1) 

•  (IUPFTS, 4) 


Coefficients  of  the 
second  of  the  second 
order  polynomial  fits 
for  y(x)  on  the  lower 
wall.  y(x)=a2+b^ 

(x-x2)+c2(x-x2)2  , 

x2^xsx3  where  a2  is 

in  columns  1-10  , 
b2  is  in  columns 

11-20, c2  is  in  coltxms 

2i-30,x  is  in  columns 
31-40  .. 

Coefficients  of  the 
last  of  the  second 
order  fits  for  y(x) 
on  the  lower  wall. 


Coefficients  of  the 
1st  of  the  2nd  order 
polynomial  fits  for 
y(x)  on  the  upper 
wall. 

y(x)=a1+b1 (x-x1)+ 

^(x-Xi)2 

x1^xsx2  where  a^^  is 

in  columns  1-10, 
bn  is  in  columns  11-20, 

c^  is  in  columns  11-30, 

x^  is  in  columns  31-40  , 

Coefficients  of  the 
last  of  the  2nd  order 
polynomial  fits  for 
y(x)  on  the  upper 
wall. 
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If  the  specified  initial  data  is  uniform  with  given  Mach  number,  there 
are  no  more  data  cards.  If  the  initial  data  is  non— uniform,  it  is  read 
on  the  following  cards.  The  first  point  is  the  point  corresponding  to 
the  minimum  value  of  y  and  will  usually  correspond  to  the  leading  edge 
of  the  lower  wall  (the  axis  point,  for  axially  symmetric  flow) .  Recall 
that  the  coordinate  system  is  right  handed,  with  x,  the  abscissa, 
positive  in  the  general  direction  of  the  flow.  The  characteristics 


issuing  from  the 

data  line 

must 

point 

downstream 

of  the 

initial  line. 

Card  Columns 

Datum 

Format 

Must  be 

Recom¬ 

Description 

No. 

Fortran 

Right 

mended 

Symbols 

Adiusted 

value 

9+ILWPTS+  1-10 

X  (1,1) 

■  E10.8 

Optional 

x,  abscissa  of  the 

IUPFTS 

1st  data  point 
(must  equal  XBEG, 
card  number  2) . 

11-20 

Y(l,l) 

y,  ordinate  of  the 
first  data  point. 

21-30 

TAU (1,1) 

Angle  (in  degrees) 
made  by  the  stream' 
line  through  the 
first  data  point. 

31-40 

SM (1,1) 

Mach  number  of  the 
first  data  point. 

41-50 

P  (1,1) 

Static  pressure 
(INPT=+1)  or  total 
pressure  (INPT=-1 
of  the  first  data 
point  (non-dimen- 
sionalized) . 

51-60 

H(l, 1) 

Static  temperature 
(INPT=+1)  or  total 
enthalpy  (INPT=-1) 
of  the  first  data 
point  (non-dimen**‘ 
sionalized) . 

9+ILWFTS+ 

Second  data  card 

IUPFTS+1 

for  initial  data 

• 

• 

line. 

• 

9+ILWFTS+ 

Last  ,  data  card 

IUPFTS+IUNFLO 

for  initial  data 

-1 

line. 

'  '-I 
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3 ,  Output 


The  output  is  done  by  writing  on  logical  tape  6  as  the  com¬ 
putation  proceeds.  After  each  down  running  characteristic  has  been 
computed,  the  data  along  that  characteristic  is  output.  The  data 
is  nondimensionalized  as  indicated  in  Part  c-a. 

Sample  output  sheets  are  given  in  Appendix  II. (pp,  126-129). 

At  the  beginning  of  the  computation,  the  pertinent  input  data 
are  output  and  identified. 

At  the  top  of  the  data  of  each  down  characteristic,  the  strip 
number  and  the  number  of  the  characteristic  in  that  strip  are 
printed  out.  These  are  foil -v'-id  by  a  block  with  important  data 
concerning  the  location  of  c^./tact  discontinuities  and  up  running 
shocks.  (See  page  127,  Appei  dix  II.)  The  contact  discontinuities  and 
the  shocks  are  numbered  from  one  (1)  to  ten  (10)  beginning  with  those 
nearest  the  upper  wall  as  one  proceeds  along  the  characteristic. 

Thus,  on  Page  128,  Appendix  II,  one  sees  that  there  is  a  contact  dis¬ 
continuity  at  the  20th  mesh  point  on  the  18th  characteristic  of  the 
3rd  strip.  The  data  on  the  'upper'  side  of  the  discontinuity  is 
at  mesh  point  20  #  while  the  data  at  the  'lower'  side  of  the  dis¬ 
continuity  is  at  mesh  point  21  .  Similarly,  there  are  up-running 

shocks  at  mesh  points  16  and  25  ,  the  data  behind  the  shocks 

being  stored  in  points  17  and  26  respectively.  The  column  under 

SIGMA  givesthe  local  slope  of  the  shock  (dy/dx  of  the  shock)  while 
1/SJGMA  is  the  inverse  of  the  slope. 

Following  is  the  data  along  the  characteristic,  beginning  at 
the  upper  wall  and  continuing  to  either  the  lower  wall  or  a  down 
running  shock,  as  the  case  may  be.  The  column  of  integers  on  the 
extreme  left  and  right  identify  the  mesh  point  number  on  that  character¬ 
istic.  The  various  quantities  output  are  the  mesh  point  coordinates, 
the  non-dimensional  pressure  and  temperature,  the  angle  (in  degrees) 
the  streamline  makes  with  the  x  axis,  the  Mach  number,  and  the  non- 
dimensional  entropy,  enthalpy,  mass  flow  (between  the  lower  wall  and 
that  mesh  point),  and  total  pressure. 

The  only  thing  to  be  cleared  up  is  how  to  be  able  to  distinguish 
whether  the  down  characteristic  ends  at  a  down  running  shock  or  bottom 
wall.  But  this  is  really  very  easy  because  the  initiation  of  a  down 
running  shock  is  always  indicated  by  one  of  the  following  three 
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messages  s 


1)  SHOCK  ON  UPPER  WALL  AT  X  =  ,  y= 

2)  REFLECTION  ON  UPPER  WALL  AT  X=  ,  Y= 

3)  DOWN  SHOCK  FORMED  BY  COALESENCE  OF  SMOOTH 
COMPRESSION  WAVES 

The  first  message  is  written  when  a  down  shock  is  formed 
by  a  sharp  compression  corner  on  the  upper  wall. 

The  second  message  is  written  when  a  down  shock  is  formed 
by  the  reflection  of  an  up  running  shock  on  the  upper  wall. 

The  third  message  is  written  when  smooth  compression  waves 
coalesce  to  form  an  'envelope  shock’. 

So,  after  any  one  of  the  three  messages,  the  last  point  on 
the  down  running  characteristics  will  be  on  the  downstream  side 
of  the  down  running  shock,  until  that  shock  reflects  at  the  lower 
wall.  This  is  indicated  by  the  message 

4)  REFLECTION  ON  LOWER  WALL  AT  X  =  ,  Y= 

indicating  that  a  down  running  shock  has  reflected  off  the  lower 
wall,  the  reflected  shock  being  an  up  running  shock.  All  the 
down  characteristics  printed  after  message  4) ,  will  terminate  on 
the  lower  wall  until  one  of  the  messages  1)  -  3)  is  again  printed. 

Sample  output  are  shown  in  Appendix  II. 
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APPENDIX  II 


SAMPLE  INPUT  AND  OUTPUT 


1 .  Sample  Input 

Sample  input  cards  for  a  case  with  the  simplified  input 
(uniform  incident  flow)  are  given  on  the  following  pages. 

The  data  for  this  case  are: 


Two-dimensional 
y  =  1.4 
M^=  12.3 

Uniform  incident  flow 

Zero  degrees  angle  of  attack 

Ay  "  0.1  on  the  initial  data  line 

lower  wall  leading  edge  at  x  =  0.0,  y  =  0.0 

upper  wall  leading  edge  at  x  =  2-75,  y  =  2.75. 

maximum  upper  wall  abscissa  at  x  =  9.5 

The  geometry  of  the  lower  wall  is: 

y  *  0.17633X  ;  0gxs:5.5 

y  =  0.96962  +  0.23087  (x-5.5)  5.5<xs8.1122 

y  =  1.5729+0.23087  (x-8. 1122) -0.097186  (x-.l . 1122) 1  8.1122<x^9.3 
y  =  1.71  x>9 . 3 


The 

y  = 
y  = 
y  = 
y  = 


geometry  of  the  upper  wall  is: 

2.75  -  0.08749  (x-2.75)  2.75^xs5.0 

2.5527  -  0.17633  (x-5.0)  5.0<xs8.782 

1.8907  -  0.17633  (x-8 . 782) +0 . 17021  (x-8. 782)^  8.782<X*9.3 
1.845  9 . 3<x 
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SAMPLE  xlnPUT  CARD 


SAMPLE  INPUT  CARDS  (cont'd) 


2.  Sample  Output 

Sample  output  sheets  for  the  case  of  the  preceding  pages 
is  given  on  the  following  pages  (up.  126-129). 

The  first  page,  Table  I,  gives  the  input  data.  If  the 
simplified  input  is  not  used,  just  disregard  those  parameters 
which  are  used  only  for  the  simplified  input.  They  are: 


MACH 

ANGLE  OF  ATTACK 

DY  Ay  on  initial  lines 

For  completeness,  a  review  of  the  meaning  of  Table  I  follows: 


The  first  two  lines  of  output  are  self-explanatory. 


IUNFLO  =  0 


means  the  incident  flow  is  uniform  and  the 
simplified  input  is  being  used 


MACH  *  12.3  M*  ■  12.3 

GAMMA  -  1.4  Y  =  1.4 

ANGLE  OF  ATTACK  ■  0  is  self-explanatory 


DY  ■  0.1  Ay  -  0.1  on  the  initial  data  line  which  is 

computed  internally 


JA  «  1 
ILWFTS  -  4 

IUPFTS  -  4 

INPT  -  0 

ICORNL  =  1 


means  the  flow  is  axieymmetric 

means  there  are  four  geometry  fits  on  the 
lower  wall 

means  there  are  four  geometry  fits  on  the 
upper  wall 

means  the  flow  is  uniform  and  the  simplified 
input  is  being  used 

means  there  is  one  sharp  corner  on  the  lower 
wall  after  the  leading  edge 
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.<  "7" 


ICORNU  =  1  means  there  is  one  sharp  corner  on  the 

upper  wail  after  the  leading  edge 

X  UPPER  =2.75  is  the  abscissa  of  the  leading  edge  of  the 

upper  wall  (cowl) 

X  FINAL  =9.5  is  the  maximum  abscissa  to  be  computed  on  the 

upper  wall 

X  IOWER  =0.0  is  the  abscissa  of  the  leading  edge  of  the 

lower  wall 

Following  are  the  abscissae  and  slopes  (dy/dx)  of  each  of  the 
corners  on  the  upper  and  lower  walls.  TATI  LEFT  is  dy/dx  on  the 
upstream  side  of  the  corner.  TAU  RIGHT  is  dy/dx  on  the  downstream 
side-  Remember,  the  leading  edges  are  not  counted  as  corners  in 
this  part  of  the  input.  ^ 

Following  are  the  coefficients  of  the  polynomial  fits  for 
the  walls.  This  data  is  self-explanatory. 

Table  II  shows  the  data  on  the  eleventh  down  running  character¬ 
istic  of  the  first  strip.  Note  that  the  leading  edge  of  the  upper 
cowl  has  not  yet  been  intersected.  Also  note  that  since  the  data 
at  x  =  0,  y  =  0,  is  on  the  first  characteristic,  one  would  expect 
to  have  the  eleventn  characteristic  having  an  ordinate  of  y  =  1.0 
since  Ay  =  0.1  was  input.  The  fact  that  y  =  0.9  here  indicates 
that  twice  a  failure  has  occurred  in  an  iteration  and  the  initial 
Ay  was  temporarily  halved.  In  this  example  this  happened  at  the 
very  beginning  near  the  symmetry  axis  of  the  conical  centerbody. 

The  lower  wal|l  has  a  leading  edge  of  10°  yielding  a  conical  shock 
of  11.9°.  For  a  more  extensive  discussion  of  how  to  read  this  table, 
see  Appendix  I,  Part  c. 

Table  III  gives  the  date,  along  a  down  running  characteristic 
further  downstream. 

The  shock  wave  pattern  for  this  example  is  given  in  Appendix 
III  along  with  some  comments  about  the  computation. 
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DPCTTT.U1C 


The  shock  wave  pattern  for  the  example  case  given  in  Appendix 
II  is  shown  in  Figure  41.  The  strip  numbers  are  shown  in  the 
Figure  in  order  to  illustrate  once  again  that  feature  of  the  program. 

Another  sample  case  is  shown  in  Figure  42.  This  is  an  in- 
ternal  flow  problem  with  many  sharp  compressions  and  expansions  and 
illustrates  very  well  the  interaction  of  the  various  wave  patterns. 
The  flow  is  two-dimensional  and  at  an  incident  Mach  number  of  3, 

P  is  the  ratio  of  the  local  pressure  to  free  stream  pressure  and  M 
is  the  local  Mach  number.  The  values  of  P  and  M  are  indicated  for 
various  points  in  the  flow  field.  Thus,  at  point  D  we  have 
P  -  0.109  and  M  =  4.63,  Point  D  is  just  in  front  of  the  intersection 
of  the  two  shocks.  Do  not  interpret  the  letters  as  indicating  that 
an  entire  area  has  the  specified  properties.  For  instance,  the 
properties  of  F  are  at  a  point  at  the  wall  slightly  in  front  of  the 
expansion  corner.  The  properties  vary  along  the  wall  between  E  and 
F. 
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APPENDIX  IV 


DISCUSSION  ON  THE  GEOMETRY  SPECIFICATION  AND 
SUGGESTIONS  ON  HOW  TO  ACHIEVE  GREATER  FLEXIBILITY 


It  has  already  been  stated  that  the  geometry  of  the  upper  and 
lower  walls  is  specified  through  a  number  of  second  order  polynomial 
fits  for  y  as  a  function  of  x.  Presently,  this  number  is  nine,  but 
it  can  be  increased  by  simply  increasing  the  dimension  of  the 
appropriate  variables.  In  addition,  there  is  provision  for  six 
corners  on  each  wall.  That  number  can  also  be  increased.  The 
dimension  of  the  appropriate  variables  must  be  one  more  than 
the  number  of  fits  or  corners,  respectively. 

It  would  not  be  difficult  to  increase  the  order  of  the  poly¬ 
nomials.  The  limitations  are  primarily  machine  storage  space  and 
convenience  in  input.  However,  the  polynomials  could  be  increased 
to  sixth  order  while  still  retaining  the  pattern  of  one  data  card 
per  polynomial.  It  is  only  necessary  to  modify  the  dimension  state¬ 
ments  for  YLBC  and  YUBC  from  (10,  4)  to  (10,  7) ,  to  change  the  cor¬ 
responding  input  and  output  statements,  and  to  modify  the  coding 
in  functions  YU'B,  YLB,  TUB,  TLB,  which  compute  y(x)  and  dy/dx  on 
the  upper  and  lower  boundaries. 

These  functions  can  be  made  as  versatile  as  desired,  the  only 
penalty  being  that  a  certain  amount  of  programming  needs  to  be  per¬ 
formed.  Each  of  these  functions  has  an  argument,  XIN,  which  is  the 
value  of  the  abscissa  at  which  one  would  like  either  y  or  dy/dx. 

That  is, 

FUNCTION  YLB  (XIN)  computes  y (XIN)  on  the  lower  boundary, 

FUNCTION  TLB  (XIN)  computes  T=dy/dx  at  x=XIN  on  the  lower  boundary 

FUNCTION  YUB  (XIN)  computes  y  (XIN)  on  the  upper  boundary, 

FUNCTION  TUB  (XIN)  computes  T=dy/dx  at  X=XIN  on  the  upper  boundary 

An  example  can  best  illustrate  the  method  for  changing  the  geometry 
routines.  Suppose  that  the  second  order  fits  are  satisfactory  on  the 
upper  wall  but  not  on  the  lower  wall  where  it  is  desired  to  specify 
y  (x)  asj 


y(x)  =0 

Osxsl 

y  (x)  =tanh  (x-1) 

I<xs2 

y(x)=0. 76159 

x>2 
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then 

y'  (x)-0  0«jx<1 

2 

y'(x)=sech  (x-1)  l<x<2 

y 1 (x) =0  x>2 

The  procedure  is  to  recompile  FUNCTION  YLB(XIN)  and 
FUNCTION  TLB(XIN)  as  follows: 

$IBFTC  YL 

FUNCTION  YLB(XIN) 

IF (XIN. LE . 1 . ) ) YLB=0 . 0 

IF (XIN . GT . 1 . ) . AND . XIN . LE . 2 . 0 ) YLB-TANH (XIN- 1.0) 

IF (XIN«GT» 2 .0) YLB=0. 76159 

RETURN 

END 

$IBFTC  TL 

FUNCTION  TLB (XIN) 

IF (XIN f LE . 1 . 0 ) TLB=0 . J 

IF  (XrN.GT .  1 . 0 . AND .XIN.LE .  2 . 0)  TLB=1  ,/COSH (XIN-1 . 0) )  **2 

ID (XIN . GT . 2 . 0 ) TLB=0 . 0 

RETURN 

END 

Do  not  forget  that  it  is  necessary  to  input  the  proper  data 
indicating  that  there  is  a  compressive  corner  at  x=1.0  and  an  ex¬ 
pansion  corner  at  x=2.0. 


133 


APPENDIX  V 

PPOGRAM  OPERATION  (THREE-DIMENSIONAL) 

These  programs  are  coded  in  the  IBM  Fortran  IV 
language.  It  is  designed  for  use  with  the  IBSYS  Monitor 
system  for  the  IBM  7090/94  digital  computer.  The  only 
tapes  utilized  are  the  system  input  tape  mounted  on  unit 
Bl  (logical  tape  5) ,  and  the  system  BCD  output  tape 
mounted  on  unit  B2  (logical  tape  6) . 

No  sense  switch  settings  are  interrogated,  nor  any 
program  stops. anticipated.  This  program  may  be  "backed" 
with  other  runs  and  may  be  pre-processed  via  the  IBM1401, 
say,  in  the  usual  manner. 

Execution  time  for  the  longer  programs  is  approxi¬ 
mately  1  minute  per  case. 


1 .  Wedge  Program 

The  wedge  program  computes  the  shock  and  inviscid 
flow  field  about  a  swept  wedge  with  supersonic  leading 
edges  including  the  effects  of  an  angle  of  attack.  The 
conditions  behind  the  reflection  of  this  shock  from  a 
plane  surface  parallel  to  the  base  of  the  wedge  as  well 
as  the  reflected  shock  geometry  are  also  computed.  In 
its  present  form  the  program  uses  a  ratio  of  specific 
heats,  y  ,  equal  to  1.4. 


Input 


The  card  input  to  the  program  is  bn  a  single  card 
and  is  read  as  follows: 

v,  M  ,  X,  ft,  Az ,  r 
00 

FORMAT  (6E10.4) 

The  symbols  in  the  input  statement  have  the  follow¬ 
ing  meanings: 

CHI  (y) :  Wedge  sweep  angle  (deg.) 

M-INF(M  )  i  Free  stream  Mach  number 
LAMBDA  (* x)tWedge  angle  in  y-z  plane  (deg.) 

ALPHA  («): Angle  of  attack  with  respect  to  wedge  base 
plane  (deg.) 

DELTA Z  ( Az) step  size  in  z  for  determining  intersection 
curve.  A  typical  value  is  0.05r. 

R(r):  Distance  between  reflection  plane  and  the  base 

plane  of  the  wedge  (normalized  in  any  convenient 
manner) . 

Output 

The  output  consists  of  the  conditions  in  the  two- 
dimensional  region  and  the  constants  that  define  the 
shock  shape,  followed  by  tables  of  conditions  in  the 
conical  region  for  each  computed  point  on  the  shock  re¬ 
flection  curve.  It  should  be  noted  that  the  conditions 
corresponding  to  the  last  point  in  these  tables  represent 
the  two-dimensional  shock  reflection. 

The  meanings  of  the  symbols  in  the  first  part  of 
the  output  are: 

LAMBDA »N(Xn) :  Wedge  angle  in  a  plane  normal  to  the 
leading  edge  (radians). 

ALPHA-N  (»h):  Sum  of  wedge  angle  and  angle  of  attack  in 
plane  normal  to  leading  edge,  i.e.,  flow 
deflection  across  two-dimensional  shock 
(radians) . 


MN-INF(m_  ): 

11 

00 

MT-INF(Mt  ): 

00 

THETA  (9^)  . 

Pl/P-XNF'fyfcJ 
MTl (Mti) : 

THETA-N  (0  )j 
In 

M-l  (ML)i 
CHI-STAR  (x*) • 

SIGMA-1  (or  )  : 

THETA-S (0  ) : 

s 

BETA-STAR (8*): 

THETA-CP ( 0 '  ): 

c 

A  j 

XJP  (x.'  )  : 

3 


Component  of  free  stream  Mach  number  in  a 
plane  normal  to  leading  edge. 

Component  of  free  stream  Mach  number 
tangent  to  leading  edge. 

Two-dimensional  shock  angle  with  respect 
to  Mr^  in  a  plane  normal  to  leading  edge 
(radians)  . 

: Pressure  jump  across  two-dimensional  shock. 

Componen  t  of  Mach  number  behind  two- 
dimensional  shock,  tanghnt  to  leading  edge. 

Shock  angle  with  respect  to  base  plane  of 
wedge  in  a  plane  normal  to  leading  edge 
(radians) . 

Total  Mach  number  on  wedge  surface. 

Angle  between  z 1 -axis  and  leading  edge 
(radians)  . 

Flow  deviation  on  wedge  surface  with 
respect  to  z'  ^xes  (radians). 

Two-dimensional  shock  angle  with  r.^ect 
to  the  base  plane  of  wedge  in  a  plane 
parallel  to  the  x-y  plane  (radians) . 

Two-dimensional  shook  angle  with  respect 
to  base  plane  of  wedge  in  a  plane  parallel 
to  the  y-z  plane  (radians) . 

Two-dimensional  shock  angle  with  respect 
to  base  plane  of  wedge  in  a  plane  normal 
to  the  z’  axis  (radians). 

Constant  in  equation  of  parabolic  sheet  (35) 

The  x  coordinate  of  the  point  of  inter¬ 
section  of  the  two-dimensional  shock  and 
the  local  Mach  cone  in  the  plane  z*  -  cos  X 


H  (h) : 


Constant  in  equation  of  parabolic  sheet  (36) 


THETA-SW  (A  )  s  Angle  betwen  the  three  dimensional  shock 
sw 

and  the  z'  axes  in  the  y-z  plane  (radians). 

The  meanings  of  the  symbols  in  the  tables  of  conditions 
on  e^ch  conical  ray  i  are: 


I(i): 


Designates  conical  ray  in  central  region. 


Z(z.) : 


X  (xi)  : 


CHI  (X  y  : 


Normalized  z  coordinate  of  the  intersection 
point  of  conical  ray  i  and  the  shock  re¬ 
flection  curve  on  plane  surface. 

Normalized  x  coordinate  of  the  intersection 
point  of  conical  ray  i  ind  the  shock  re¬ 
flection  curve  on  plane  surface. 

Angle  between  the  tangent  to  the  inter¬ 
section  curve  and  the  z-axes  (radians) . 


SIGMA-1  (  )  : 

i 

M-l (Mx  ) : 
i 

MT  (M  )  : 

it . 

1 

MN  (M1  )  : 


MKN 


<V 


Flow  deviation  on  conical  ray  i  (radians) . 

Total  Mach  number  on  conical  ray  i. 

Upstream  Mach  number  component  tangent 
to  intersection  curve. 

Upstream  Mach  number  component  normal  to 
intersection  curve,  on  plane  surface. 

Total  Mach  number  component  in  a  plane 
normal  to  intersection  curve. 


DELTA-1  (6l  ): 

i 


Flow  deflection  across  reflected  shock 
(radians) . 


M2T  (M-  ) 


Component  of  Mach  number  behind  reflected 
shock,  tangent  to  intersection  curve. 


P2/P-INF 

(p2i/pj 

M2  (M2i  ) : 


Pressure  ratio  behind  reflected  shock. 
Total  Mach  number  behind  reflected  shock. 


SIGMA-2  (<r2i)  : 


Flow  deviation  with  respect  to  z  axes 
behind  reflected  shock  (radians) . 


THETA-P  Reflected  shock  angle  with  respect  to  the 

(q  ,).  reflection  plane  in  a  plane  normal  to  the 

2ni  intersection  curve  (radians) . 

The  subroutine  which  contains  the  oblique  shock 
relations  is  equipped  with  four  stops  printed  out  as 
follows: 


1)  M  =  1:  If  the  total  Mach  number  becomes  sonic. 

oo 

2)  Upstream  normal  Mach  number  <1:  If  the  Mach  number 
component  in  a  plane  normal  to  the  leading  edge  or 
to  the  intersection  curve  becomes  sonic. 


3)  Mach  Reflection:  If  the  flow  deflection  used  in  the 
shock  relations  exceeds  the  maximum  value  given  by 


4 

3  vT" (yf r) 


%2-  n3/2 

V 


where  M„  is  the  Mach  number  component  in  a  plane  normal 
to  the  leading  edge  or  to  the  intersection  curve.  The 
computation  of  the  intersection  curve  is  terminated  but 
the  rest  of  the  computation  is  continued  for  the  points 
already  obtained. 


4)  Imaginary  roots:  If  a  pair  of  roots  of  the  cubic 
equation  for  the  shock  angle  are  imaginary. 


;& 
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2  Conical  Shock-Plane  Surface  Program  and  Characteristic 
Subroutine 

This  program  starts  with  the  shock  and  flow  field 
about  a  cone  at  zero  angle  of  attack  and  computes  the 
reflection  of  this  shock  from  a  plane  inclined  at  an 
angle  0^  to  the  cone  axes.  The  characteristic  subroutine 
computes  the  flow  field  behind  the  reflection  curve  on 
the  plane  surface  by  using  two-dimensional  characteristics 
and  may  be  applied  to  any  of  the  programs  of  the  report. 

The  ratio  of  specific  heats,  y  ,  is  an  input  quantity  in 
the  program,  while  the  characteristic  subroutine  auto¬ 
matically  uses  the  value  which  the  main  program  has  used. 

Main  Program 


Input 


The  card  input  to  the  program  is  on  two  cards  and 
is  read  as  follows: 


CARD  1: 
CARD  2: 


M»'  V  V  V'  r>  q>  Az‘ 

FORMAT  (7E  10.4) 

Pc/Peo'  M*1C,  01c  i  S^/R 
FORMAT  (4E  10.4) 


The  symbols  in  the  input  statement  have  the  follow- 


ing  meanings: 

M-INF  (M  ): 

03 

Free  stream  Mach  number. 

THETA-1 (0^) . 

Angle  between  plane  surface  and 
stream  direction  (deg.) 

free 

GAMMA-C  (v  )  j 
c 

Conical  shock  angle  from  Reference  15  (rad.) 

R  (r): 

y  distance  between  origins  0  and  O'  (Fig. 31) 
normalized  in  any  convenient  manner. 

Q  (q): 

Nondimens ionalized  z  distance  between- 0 
and  O',  usually  zero. 

DELTA-Z(  Aa"1 )  : 

Step  size  in  z'  for  determining 
section  curve.  A  typical  value 

inter- 
is  0.05  r. 
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PC/P-INF 
(pc/p.) > 


Pressiirp  ratio  across  conical  shock,  from 
Reference  15. 


Ml-STAR  (M*  )  : 

c 


Critical  Mach  number  behind  conical  shock, 
from  Reference  15. 


PSI  ($lc) :  Angle  of  flow  behind  conical  shock  with 

respect  to  the  cone  axes  (rad.). 

S2/R(S2/R) :  Entropy  change  across  conical  shock 

normalized  with  respect  to  the  gas  constant. 


Output 

The  output  of  the  program  consists  of  tables  of  con¬ 
ditions  in  front  of  and  behind  the  reflected  shock,  near 
the  plane  surface,  for  each  point  i  on  the  reflection 
curve . 

The  meanings  of  the  symbols  in  the  output  are: 

I(i):  Designates  point  on  reflection  curve. 

ZP(z'.):  Normalized  z'  coordinate  of  point  i  on 

reflection  curve. 


XP  (x1.);  Normalized  x*  coordinate  of  point  i  on 

reflection  curve. 


OMEGA  (0i)  :  Angle  between  the  tangent  of  the  reflection 

curve  at  point  i  and  the  z'  axes.  The 
computation  terminates  whenQ^<5  (rad.). 


M2XP 

(M2x'.)! 

1 

M2YP 

<v.): 

J  X 

M2ZP 

(M,  )  : 

z '  . 

l 

M2N 

(M2ni>  ! 

Mach  number  component  behind  conical  shock 
in  x'  direction. 

Mach  number  component  behind  conical  shock 
in  y'  direction. 

Mach  number  component  behind  conical  shock 
in  z'  direction. 

Mach  number  component  in  front  of  the  re¬ 
flected  shock,  normal  to  the  reflection 
curve,  on  the  plane  surface. 
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M?T  (M2  )  : 


Mach  number  component  in  front  of  reflected 
shock,  tangent  to  the  reflection  curve. 


M2KN  (M2n  ): 

i 


Total  Mach. number  component  in  front  of 
the  reflected  shock,  in  a  plane  normal  to 
the  reflection  curve. 


DELTA- 2 


<52i>‘ 


Flow  deflection  across  the  reflected  shock 
in  a  plane  normal  to  the  reflection  curve 
(radians) . 


P4/P2  (p^/p2)  ^-s  Pres  sure  jump  across  reflected  shock. 


A4/A2  (a./a.).  •  Speed  of  sound  ratio  across  reflected 
1  shock. 


M4N  (M4n  )= 
i 


Mach  number  component  behind  reflected 
shock  in  a  plane  normal  to  the  reflection 
curve . 


THETA  4  (A  )  . 

4 . 

i 

M4T  (M4t  ): 
i 


Reflected  shock  angle,  in  a  plane  normal 

to  the  reflection  curve,  with  respect  to 

the  direction  of  M.,  ,  (radians) . 

2N . 

X 

Mach  number  component  behind  reflected 
shock,  tangent  to  reflection  curve. 


P4/P-INF 

(p4  /p.  ) :  Pressure  ratio  behind  reflected  shock, 

i 

SIGMA-4(cr4^)  .  Flow  deviation  behind  reflected  shock, 
"L  with  respect  to  z*  axes  (radians). 


THETA-P-4  ( 0'  )  ; 

4 . 

l 


Reflected  shock  angle  with  respect 
to  plane  body  (radians) . 


S(S4  ): 
i 


Change  'in  entropy  from  the  free  stream 
vakie,  nondimensionalized  by  twice  the 
specific  heat  at  constant  pressure  (2c)* 


Characteristic  Subroutine 


The  Characteristic  Subroutine  is  an  integral  part 
of  the  progran  of  this  section,  Conical-Shock-Plane  Body, 
and  the  program  of  Section  V,  Plane  Shock-Conical  Body. 

In  these  cases  no  additional  input  is  needed.  The  sub¬ 
routine  also  exists  as  a  separate  program  for  which  the 
input  is  given  below.  The  output  of  the  subroutine  is 
the  same  for  all  cases  and  is  also  given  below. 

Input 

When  used  as  a  separate  program,  the  input  is  on 
N+l  cards  and  is  read  as  follows: 

CARD  1:  y,  N 

FORMAT  (2E10. 4,  HO) 

CARD  2:  p ^  0^,  SN  (Data  at  last  point  of 

main  program) w 

.  (One  card  for  each  data  point) 


CARD  N+l:  M^,  p^,  0^,  (Data  at  first 

point  of  main  program) 

The  symbols  in  the  input  statement  have  the  follow¬ 
ing  meanings: 


GAMMA  (y)  :  Ratio  of  specific  heats 

P/P  (PN  jj/p^)  : Pressure  ratio  at  point  N, N  (Fig.32)  . 

N:  Number  of  data  points. 


M  (M.  .): 
i.D 


Total  Mach  number  at  data  point,  i=j. 


RHO(p  ): 

X  /  J 


Nondimensionalized  coordinate  of  data- 
point  (Fig.32). 


PSI  (ib.  .)  : 


Nondimensionalized  coordinate  of  data 
point  (Fig.32.) 


THETA  (8. 

1.3 


Plow  deviation  with  respect  to  p  axes, 
at  data  point  (radians) , 


(S.  .): 
1.3 


Entropy  change  from  free  stream  at  data 
point,  normalized  with  respect  to  twice 
the  specific  heat  at  constant  pressure. 


Output 


The  output  consists  of  the  Mach  angle  and  velocity 
at  the  data  points,  which  is  printed  next  to  the  data, 
and  a  page  of  output  for  each  second  family  characteristic 
line,  i=constant,  which  gives  the  flow  conditions  at  each 
point  along  the  line. 

The  output  symbols,  not  defined  under  input,  have 
the  following  meanings: 


MU  (p.  .)  : 

1.3 

W  (W.  .)  : 

i.3 


Mach  angle  (radians) . 

Flow  velocity,  normalized  with  respect  to 
the  limiting  velocity. 


I  (i): 


Second  family  characteristic  line.  The 
computation  starts  at  I=N,  Fig.  32. 


P/P-INF (p .  ./p  ) : 

1  9  J  " 


Pressure  ratio  at  the  wall  points, 
i.e.,  last  point  on  each  i  line. 


J  (j):  First  family  characteristic  line,  Fig. 32. 
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3 .  Double  Conical  Shock  program 

The  double  conical  shock  program  computes  the  flow 
field  behind  the  intersection  of  two  conical  shocks  of 
different  strength  when  given  the  flow  conditions  behind 
the  shocks.  An  option  is  included  which  allows  one  of 
the  cones  to  be  oriented  by  an  angle  attack,  a,  with  re¬ 
spect  to  the  free  stream  direction.  The  properties 
immediately  behind  the  upper  and  lower  shocks  behind  the 
intersection  are  found,  along  with  the  new  shock  angles. 
This  data  can  be  used  in  the  characteristic  subroutine, 
(Section  C) ,  to  find  the  flow  field  on  a  plane  passing 
through  the  shock  intersection. 

Input 


The  card  input  consists  of  three  cards  if  both  cones 
are  at  zero  angle  of  attack  and  four  cards  if  one  cone 
is  at  an  angle  of  attack,  and  is  read  as  follows: 


CARD 

1: 

Ml*,  *!*,  vr  Pj/P.,  S/R 

FORMAT  (5E10.4) 

CARD 

2: 

m*  o--  </>?  o  '  p->/p  /  s,  0/R 

FORMAT  (5E10.4) 

CARD 

3 : 

r,q,A$,y,  S2/R,KALF 

FORMAT  (5E10.4,  IlO) 


CARD  4:  a,  ij/a,  M2*,  M*  2>>P2* 

FORMAT  (5E1Q.4) 

(Only  required  if  the  angle  of  attack  a  /  0) 

The  symbols  in  the  input  have  the  following  mean¬ 
ings: 

Ml*  (M^*) :  Critical  Mach  number  behind  the  shock  of 

the  first  cone,  ("base  cone"),  which  has 
zero  angle  of  attack. 

PSI  1*  (j/j^*)  Flow  direction  angle  behind  the  shock  of 
the  base  cone,  (radians). 
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GAMMA  1  (y^) .  Base  cone  shock  angle,  (radians). 

Pl/P-INF  (Pq/P^ ) 1  Pressure  jump  across  base  cone  shock. 


Si/K  (5.^/R) 


tncropy  increase  across  base  co.'.e  shock 
normalized  with  respect  to  the  universal 
gas  constant,  Reference  15. 


M12*  (M*  2) :  Critical  Mach  number  behind  the  shock  of 

the  second  cone,  for  zero  angle  of  attack. 


PS12  (jfj*  ) 


Flow  direction  angle  behind  the  shock  of 
the  second  cone,  for  zero  angle  of  attack, 
(radians) . 


GAMMA  2  (v2) :  Shock  angle  of  second  cone  at  zero  anqle 
of  attack,  (radians)  . 


P2/P-INF (p^/p  ):  Pressure  ratio  across  the  second  cone 

shock,  for  zero  angle  of  attack. 

S12/R(S^  / R ) :  Entropy  increase  across  second  cone  shock, 

'  for  zero  angle  of  attack. 


R(r) : 


Vertical,  (y) ,  distance  between  cone 
verticies,  normalized  in  any  convenient 


manner . 


Q  (q) : 


Normalized  z  distance  between  cone  verticies. 


DELTA-PSI  (A0)  :  Step  size  in  ip  for  determining  the  inter¬ 
section  curve.  A  typical  value  is  5  deg. (deg.) 

GAMMA-BAR  (y)  :  Ratio  of  specific  heats. 


S2/R  (S2/R) 


Entropy  perturbation,  due  to  angle  of 
attack  of  second  cone^. 


HALF: 


Index  for  including  angle  of  attack  of 
second  cone;  =0,  zero  angle  of  attack;  ^  0, 
nonzero  angle  of  attack. 


ALPHA  (a)  ; 


Angle  of  attack  of  second  cone,  (deg.) 


ETA/ALPHA  (  v/et)  :  Ratio  of  tilt  angle  of  second  cone  shock 
to  angle'  of  attack  1®. 


M2*  (M2*): 


Critical  Mach  number  perturbation,  due  to 
angle  of  attack,  behind  second  cone  shock. 
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£ 

v> 


MW2*(M*  Critical  Mach  number  normal  to  meridian 

"*  plane,  behind  second  cone  shock. 


PSI2*(02*) # 


Plow  direction  angle  perturbation,  due  to 
angle  of  attack,  behind  second  cone  shock, 
(radians)  . 


Output 

The  output  consists  of  the  maximum  value  of  ip,  three 
quantities  which  are  invariant  along  the  intersection  curve 
and  tables  which  describe  the  intersection  curve,  the  base 
and  second  cone  flows  and  the  flows  in  the  upper  and  lower 
intersection  regions,  for  each  point  i  on  the  intersection 
curve . 


The  meanings  of  the  symbols  in  the  output  are: 


PSI-MAX (i/jmax . ) :  Maximum  value  of  intersection  curve 

coordinate,  ip,  before  reaching  double 
values  of  p,  (degrees) . 

THETA-1  (9^) :  Angle  between  second  cone  shock  axis  and 

free  stream  direction,  (radians) . 

COS(BZ)  (cos(bz)):  Direction  cosine  of  binormal  to  intersection 

curve  and  z  axes. 


MZi  (M  ): 
Z1 


Mach  number  component  in  z  direction,  behind 
lower  cone  shock . 


PSI  (ip±)  : 


Meridian  coordinate  of  point  on  intersection 
curve,  (radians). 


COS(TX)  (cos(tx).)  sDirection 

curve  and 


cosine  of  tangent  to  intersection 
x  axes. 


COS (TY)  (cos (ty ^)  sDirection  cosine  of  tangent  to  intersection 

curve  and  y  axes. 

COS(TZ)  (c$os  (tz).)  sDirection  cosine  of  tangent  to  intersection 

curve  and  z  axes. 


COS  (3X)  (cos(bx).)  sDirection 

curve  and 


cosine  of  binormal  to  intersection 
x  axes. 
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COS  (BY)  (cos  (by).)  : Direction  cosine  of  binormal  to  intersection 

curve  and  y  axes. 

COS  (NX)  (cos  (nx).)  : Direction  cosine  of  normal  to  intersection 

curve  and  x  axes. 

COS  (NY)  (cos  (ny)^)  {Direction  cosine  of  normal  to  intersection 

curve  and  y  axes. 

COS  (NZ)  (cos  (nz).)  {Direction  cosine  of  normal  to  intersection 

curve  and  z  axes. 

RHO  (p.){  Nondimensionalized  radial  coordinate  of 

i 

point  on  intersection  curve. 

GAMMA  (F. )  :  Angle  between  the  normal  to  the  intersection 

curve  and  the  free  stream  direction,  (radians) . 

Mach  number  component  in  y  direction  behind 
base  cone  shock. 

Mach  number  component  in  v.  direction  behind 
base  cone  shock. 

MT-BAR  1  (M£  ) :  Mach  number  component  tangent  to  the  inter - 

^i  section  curve,  behind  the  base  cone  shock. 

MN-BAR  1  (M-  ) :  Mach  number  component  along  the  normal  to 

nii  the  intersection  curve,  behind  the  base 
cone  shock. 

MB-BAR  1  (Mg  ) :  Mach  number  component  along  the  binormal  to 

the  intersection  curve,  behind  the  base 
cone  shock. 

) {  Total  Mach  number  component  in  a  plane  normal 

to  intersection  curve,  behind  the  base  cone 
shock  upstream  of  the  intersection. 

DELTA  1  (6^):  Strength  of  the  base  cone  shock  in  a  plane 

1  normal  to  the  intersection  curve,  and  which 
is  assumed  to  be  constant  through  the  inter¬ 
section,  (radians)  . 


MNl  (Mn^ 


MYl  (My  ) : 


MXl  (KL.  )  : 


147 


X  (x  ' )  : 


Nondimensionalized  x'  coordinate  of  inter¬ 
section  curve  point  i. 


Y  (y  ' )  : 


nondimensionalized  y*  coordinate  of  inter¬ 
section  curve  point  i. 


Z  (z‘ )  i 


nondimensionalized  z'  coordinate  of  inter¬ 
section  curve  point  i. 


PHI  (tfi.)i 

i 

U-BAR  (u. ) s 

i 

V-BAR  (v . ) : 

i 

W-BAR  (w.): 

l 

M*(Mi*) : 


Meridian  angle  of  intersection  curve  point, 
(second  cone  coordinates) ,  (radians) . 

Velocity  component  behind  the  second  cone 
shock  along  the  z  axes,  normalized  with 
respect  to  limiting  velocity. 

Normalized  velocity  component  behind  the 
second  cone  shock  along  the  r  axes. 

Normalized  velocity  component  behind  the 
second  cone  shock,  normal  to  the  meridian 
plane. 

Total  critical  Mach  number  behind  the  second 
cone  shock. 


A2 

^a2.  ^  : 
l 

MX2 

(M 

2. 

l 

MY2 

(My 

y2. 

i 

MZ2 

): 

): 

) 


Speed  of  sound  behind  the  second  cone  shock, 
normalized  by  the  limiting  velocity. 

Mach  number  component  in  x  direction,  behind 
second  cone  shock 

Mach  number  component  in  y  direction,  behind 
second  cone  shock . 

Mach  number  component  in  z  direction  behind 
second  cone  shock. 


MT-BAR  2 


Mach  number  component  tangent  to  the  inter¬ 
section  curve,  behind  the  second  cone  shock. 


MN-BAR  2  (M-  ) : 

n2. 

l 


Mach  number  component  normal  to  the  inter¬ 
section  curve,  behind  the  second  cone  shock. 


i- 
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MB -BAR  2 


Mach  number  component  along  the  binormal 
to  the  intersection  curve,  behind  the 
second  cone  shock. 


MN2  (Mn 


2 


Total  Mach  number  component  in  a  plane 
normal  to  the  intersection  curve,  behind 
the  second  cone  shock. 


DELTA  2  ($2  ) :  Strength  of  the  second  cone  shock  in  a 

i  plane  normal  to  the  intersection  curve, 

and  which  is  assumed  constant  through  the 
intersection,  (radians)  , 


MBi  (M  )  : 

Bl. 

l 


Mach  number  component  behind  lower  portion 
of  intersection,  in  a  plane  non;  al  to  the 
intersection  curve. 


THETA  (A  )  : 

I 


Second  cone  shock  angle  with  respect  to  6-^  , 

i 

in  a  plane  normal  to  intersection  curve,  be¬ 
hind  the  intersection,  (radians)  . 


M3T  (M^  ) i  Mach  number  component  behind  lower  portion  of 

t.  intersection,  tangent  to  the  intersection 

x  curve. 


M3  (M  )  : 

l 


Total  Mach  number  behind  lower  portion  of 
intersection. 


P3/P-INF (p3  /pj) :  Pressure  ratio  immediately  behind  lower 
i  portion  of  intersection. 


S3  (S3  ) s 
i 


MB2  (Mfi  ) 


THETA  (f>4  )  : 

i 


Entropy  change  from  free  stream,  behind 
lower  portion  of  intersection,  normalized' 
by  twice  the  specific  heat  at  constant 
pressure. 

Mach  number  component  behind  upper  portion 
of  intersection,  in  a  plane  normal  to  the 
intersection  curve. 

Base  cone  shock  angle,  with  respect  to  62  , 

in  a  plane  normal  to  the  intersection  x 
curve,  behind  the  intersection,  (radians) . 


149  • 


150 


M4T  (M  ) 

4. 

1 


Mach  number  component  behind  upper  portion 
of  intersection,  tangent  to  the  inter¬ 
section  curve. 


M4  (M4  ): 
i 


Total  Mach  number  behind  upper  portion  of 
intersection. 


P4/P-INF  (p^ /p^) ; Pressure  ratio  behind  upper  portion  of 
i  intersection. 


S4  (S4  ).  : 
i 


Entropy  change  from  free  stream  behind  upper 
portion  of  intersection,  normalized  by  twice 
the  specific  heat  at  constant  pressure. 


4. 


Plane  Shock  -  Conical  Body 


When  given  a  conical  flow  field  and  a  free  stream  Mach 
number,  this  program  computes  the  flow  behind  the  shock 
from  a  two-dimensional  wedge,  both  before  and  after  its 
intersection  with  the  conical  shock.  The  reflection 
curve  of  the  modified  plane  shock  on  the  conical  body  is 
then  determined,  together  with  the  flow  field  behind 
the  reflected  shock  on  the  cone  surface.  The  characteristic 
subroutine,  described  in  Section  C,  is  incorporated  in  the 
program  to  compute  the  rest  of  the  flow  field  on  the 
conical  body. 

Input 

The  card  input  to  the  program  consists  of  four  cards 
and  it  is  read  as  follows: 


CARD  Is 


CARD  2: 


CARD  3: 


M  ,  6. ,  r ,  q,  AZ,  y 

CO  i 

FORMAT  (6E10.4) 

PC/P„.  VP-'  VT.»'  S2/R 
FORMAT  (4E10.4) 

V  M*lc'  *lc 

FORMAT  (3E10.4) 


CARD  4s  Vjj,  M*ib' 

FORMAT  (3E10.4) 

The  meanings  of  the  symbols  in  the  input  statement  are: 


M-INF  (M  ):  Free  stream  Mach  number. 

00 

DELTA-1  (6^):  Two-dimensional  wedge  angle,  (deg.) 

R  (r) :  y  distance  of  the  wedge  vertex  from  the  cone 

vertex,  normalized  in  any  convenient  manner, 
(Figure  38)  . 


Q  (q) : 


Normalized  z  distance  of  the  wedge  vertex 
from  the  cone  vertex. 


DELTA-Z  (Az1):  Step  size  in  s'  for  determining  inter¬ 

section  curve. 

GAMMA  (y)  :  Ratio  of  specific  heats. 

PC/P-INF  (pc : Pressure  ratio  behind  conical  shock. 

PB/P-INF  :Pressure  ratio  at  cone  surface. 

TC/T-INF  (T iTemperature  ratio  behind  conical  shock. 

S2/R  (S^/R) :  Entropy  increase  across  conical  shock, 

normalized  with  respect  to  the  universal 
gas  constant. 

GAMMA  (y  ,  y^) :  Conical  shock  and  cone  angles,  respec¬ 
tively,  (radians). 

Ml-STAR  (M*^  M*^) :  Critical  Mach  numbers  behind  conical  shock 
'  and  on  cone,  respectively. 

PSI-1  j/j^)  j  Flow  deviation  behind  conical  shock  and 
'  on  cone,  respectively,  (radians)  . 


Output 

The  output  consists  of  many  parts.  The  first  includes 
some  quantities  in  the  two-dimensional  and  conical  fields, 
which  are  independent  of  the  shock  intersection.  Second, 
are  tables  describing  the  intersection  curve  of  the  plane 
and  conical  shocks,  and  the  flow  conditions  ahead  and 
behind  the  modified  plane  shock  near  the  intersection.  Next, 
the  reflection  curve  of  the  modified  plane  shock  on  the  cone 
surface  is  described,  along  with  the  flow  conditions  in 
front  and  behind  this  shock  near  the  body.  Finally,  the 
conditions  behind  the  reflected  shock  on  the  cone  surface 
are  given.  The  characteristic  subroutine  is  included  in 
the  program  and  its  output  is  described  in  Appendix  V-2. 

The  meanings  of  the  symbols  in  the  output  are: 

GAMMA-C  (y  ) :  Cone  shock  angle,  (radians). 

MR  (M  ) :  Mach  number  component  behind  the  conical 

shock,  parallel  to  the  shock. 
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Mx  (M^) : 


Mach  number  component  behind  the  conical 
shock,  normal  to  the  shock. 


M-l  (M^) :  Mach  number  behind  the  two-dimensional  shock 

THETA-1  .  Two-dimensional  shock  angle,  with  respect 

to  the  free  stream  direction,  (radians) . 


Pl/P-INF(p./p  ): 

-L  oo 

Al/A-INF (a, /a  ): 

1  « 


ZP  (2 ' ±)  : 


XP  (x 1 ±)  : 


OMEGA  (Qi)  : 


M2XP  (M  )  : 

x'  . 

i 

M2YP  (M  )  : 

y'i 

M2ZP (M  )  : 
z’  . 

l 

M2NP  (M2  )  : 

n . 

i 


M2TP  (M2  )  : 


Pressure  jump  across  two-dimensional  shock. 

Speed  of  sound  increase  across  two-dimen¬ 
sional  shock. 

Normalized  z'  coordinate  of  point  on  two- 
dimensional  and  conical  shock  intersection, 
on  plane  of  two-dimensional  shock. 

Normalized  x'  coordinate  of  point  on  shock 
intersection  curve. 

Angle  between  the  tangent  to  the  inter¬ 
section  curve  and  the  z'  axes.  The  com¬ 
putation  of  the  curve  terminates  when 
CL<5°  ,  (radians)  . 

Mach  number  component  in  x'  direction, 
behind  conical  shock. 

Mach  number  component  in  y'  direction, 
behind  conical  shock. 

Mach  number  component  in  z*  direction, 
behind  conical  shock. 

Mach  number  component  normal  to  the 
intersection  curve  in  the  plane  of  the 
two-dimensional  shock,  behind  the  conical 
shock . 

Mach  number  component  tangent  to  the  inter¬ 
section  curve,  behind  the  conical  shock. 


M2KN  (M  ) : 

2N. 

i 


Total  Mach  number  component  in  a  plane 
normal  to  the  intersection  curve,  behind 
the  conical  shock. 


DELTA  (2) s 


OMEGA- STAR: 
GAMMA: 


Plow  deviation  behind  conical  shock  with  respect 
to  plane  of  two-dimensional  shock,  in  a 
plane  normal  to  the  intersection.  (First 
term  of  Eq.  138  )  ,  (radians)  . 

v/2-0  ,  (radians). 

Free  stream  flow  deviation  with  respect  to 
the  plane  of  the  two-dimensional  shock,  in 
a  plane  normal  to  the  intersection  curve. 
(Second  term  of  equation  138  )  ,  (radians) . 


DELTA  (2)  STAR  (6*,^.):  Strength  of  the  conical  shock  in  a  plane 

normal  to  the  intersection  curve,  (radians) . 


DELTA  (1) :  Flow  deviation  behind  two-dimensional  shock 

with  respect  to  two-dimensional  shock  plane, 
in  a  plane  normal  to  the  intersection  curve. 
(Second  term  of  equation  135  ')  (radians)  . 


DELTA (1) STAR( 6*  ): Strength  of  the  plane  shock  in  a  plane  normal 
l  to  the  intersection  curve,  (radians). 


P4/P2(p4/p2)iS 
A4/A2  (a4/a2)  i! 


Pressure  jump  across  modified  plane  shock. 

Speed  of  sound  increases  across  modified 
plane  shock. 


M4N  (M.  ) : 


N. 

l 


Mach  number  component  in  a  plane  normal  to 
the  intersection  curve,  in  the  lower  inter¬ 
section  region. 


THETA (4)  (A..):  Modified  plane  shock  angle  behind  inter- 
1  section,  with  respect  to  upstream  conical 
flow  direction,  in  a  plane  normal  to  the 
intersection  curve,  (radians) . 


M4T  (M.  ) : 

t. 

i 


Mach  number  component  tanget  to  the  inter¬ 
section  curve,  in  the  lower  intersection 
region. 
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M4  (M,  ) !  Total  Mach  number  in  the  lower  intersection 

^ . 

1  region. 

P4/P-INF (p,/p  ),  jPressure  ratio  in  the  lower  intersection 

4  oo  i 

region . 

A4/A-INF  (a  . /a  )  j  Speed  of  sound  ratio  in  the  lower  inter - 

i  section  region. 

SIGMA(4)(ct^  ):  Plow  deviation  in  the  lower  intersection 

i  ■  region  with  respect  to  the  z'-axes,  (radians). 


THETAP(4)  (4*  ):  Lower  shock  angle  with  respect  to  the  £  axes, 

i  Figure  39,  (radians). 


S4  (S4  ): 
i 


Entropy  increase  from  free  stream,  in  the  lower 
intersection  region,  normalized  by  twice  the 
specific  heat  at  constant  pressure. 


GAMMA  (K)  (y  )  :  Conical  body  half  angle,  (radians). 

I  (i) s  Shock  -  conical  body  intersection  points. 


ETA  (r).  )  :  Normalized  coordinate,  normal --to  plane  of 

L1  the  two-dimensional  shock,  of  intersection 

point  of  modified  plane  shock  and  conical 
body,  (Figure  39) . 


ZETA  (£  .  ) : 
11 


Normalized  coordinate  of  intersection  point, 
in  plane  of  two-dimensional  shock  normal  to 
intersection  curve. 


X  (x.  ): 

11 


Normalized  x  coordinate  of  intersection 
point. 


Y 


): 


Normalized  y  coordinate  of  intersection 
point . 


Z 


)  s 


Normalized  z  coordinate  of  intersection 
point . 


RHO  (Pi): 


Normalized  radial,  conical  coordinate  of 
intersection  point,  (Figure  40) . 
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PHI  (<pi)  : 


Meridian,  conical  coordinate  of  inter¬ 
section  point,  (radians)  . 


CHECK: 


Checks  if  x.  ,  y.  ,  z.  ,  lie  on  the 

al  ~  1 1  1  o  2 

conical  body  (x^,  +  y  .  -(tan  v.)  z^  =0) 

1  \  b  1 


THREE  POINT 
SELECTION  : 


I=a,  I=b  ,  I=c  :  Intersection  points 

chosen  for  determining  reflection  curve. 
(The  program  presently  uses  1=1,  1=4 
1=10) . 


RH0-0,RH0-1 ,RH0-2,  Normalized  coordinates  of  the  three  sel- 
{pQ,  p.,  p2)  •.  ected  points. 

PHI-1,  PHI-2,  Cooordinates  of  last  two  selected  points. 

(tp.  ,  /o_)  :  The  coordinate  w  =  0,  (radians). 

o 

COS  (TX)  (cos  (^x).): Direction  cosine  of  tangent  to  reflec¬ 
tion  curve  and  x  axes. 


COS(TY)  (cos  (ty).)  •  Direction  cosine 

tion  curve  and  y 


of  tangent  to  reflec- 
axes. 


COS  (TZ)  (cos  (£z)^)  ‘Direction  cosine 

tion  curve  and  z 


of  tangent  to  reflec- 
axes. 


COS(BX)  (cos  (lax).)*  Direction  cosine 
1  tion  curve  and  x 


of  binormal  to  reflec 
axes . 


COS  (BY)  (cos  (lay  )^):  Dir  action.,  cosine 

tion  curve  and  y 


of  binormal  to  reflec- 
axes . 


COS(BZ)  (cos  (^z)^) ;  Direction  cosine 

tion  curve  and  z 


of  binormal  to  reflec- 
axes . 


COS  (NX)  (cos  (roc).). Direction  cosine 

tion  curve  and  x 


of  normal 
axes . 


to  reflec- 


COS  (NY)  (cos  (r^).  ): Direction  cosine 

tion  curve  and  y 


of  normal  to  reflec- 
axes . 
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A 

COS (NZ) (cos (nz) . )  :  Direction  cosine  of  normal  to  reflection 

curve  and  z  axes. 


M2XP (M_  ) : 

Mach  number  component, 

in 

x 1 

direction, 

on 

x ' 

bx 

conical  surface. 

M2YP (M-  ) s 

'~.ach  number  component , 

in 

y* 

direction, 

on 

y  bi 

■'•'nical  surface. 

M2ZP (M-  ) : 

Mach  number  component, 

in 

z1 

direction, 

on 

conical  surface. 

M2T  (M. 


'bi 


Mach  number  component  tangent  to  intersection 
curve,  on  conical  surface. 


M2N(M2  ): 
"bi 


Mach  number  component  normal  to  intersection 
curve,  on  conical  surface. 


M2KN (M  ) 


N, 


bi 


Total  Mach  number  component  in  a  plane  normal 
to  the  intersection  curve,  on  the  conical 
body . 


P4/P2(p4/p2)b.: 


Pressure  jump  across  the  modified  plane  shock 
near  the  cone. 


A4/A2(a4/a2)bi: 


Speed  of  sound  inciease  across  the  modified 
plane  shock  near  the  cone. 


M4CN(M4 


Total  Mach  number  component  in  a  plane 
normal  to  the  intersection  curve,  in  front 
of  the  reflected  shock,  near  the  cone. 


THETA  (4)  (6  ): 


Modified  plane  shock  angle  near  the  cone, 
with  respect  to  the  direction  of  M2  , 
(radians)  . 


S4(S  ): 

b. 

1 

M4T (M4  )  : 

fcbi 


Total  entropy  change  behind  the  modified 
plane  shock  near  the  cone,  normalized  by 
twice  the  specific  heat  at  constant 
pressure . 

Mach  number  component  tangent  to  the  inter¬ 
section  curve,  in  front  of  the  reflected 
shock,  near  the  cone. 


157  ■ 


M4N(M4  ): 

K» 

“bi 


Mach  number  component  along  normal  to 

1  i  An  i  n  Pv*r»w4-  n -F  — 

•Vi*  VV.1.  UV-V^-kVin  U4.  V  W  ,  ^11  X.  Am  V11W  U  JL  JL.  'wX  XUU  ~ 


ted  shock,  near  the  cone. 


M4B (M.  ) : 


bi 


Mach  number  component  along  binormal  to 
intersection  curve,  in  front  of  reflec¬ 
ted  shock,  near  the  cone. 


M4TB (M4_  ) : 


Mach  number  component  tangent  to  the 
reflection  curve,  in  front  of  the  reflec¬ 
ted  shock,  near  the  cone. 


M4VB(M4J  : 
n. 

l 


Mach  number  component  along  the  normal 
to  the  reflection  curve,  in  front  of  the 
reflected  shock, near  the  cone. 


(v  -i 


Mach  number  component  along  the  binormal 
to  the  reflection  curve,  in  front  of  the 
reflected  shock  near  the  cone. 


DELTA-4 (64  )s 
i 


Plow  deflection  across  the  reflected  shock 
in  a  plane  normal  to  the  reflection  curve, 
(radians)  . 


M4KN(M4_  ) : 


Total  Mach  number  component  in  a  plane 
normal  to  the  reflection  curve,  in  front 
of  the  reflected  shock. 


THETAP-6 (0  g  )  : 

i 


Reflected  shock  angle  with  respect  to  the 
cone  surface,  in  a  plane  normal  to  the  re¬ 
flection  curve  (radians)  . 


P6/P4(p6/p4)is 


Pressure  jump  across  reflected  shock  on  cone 
surface . 


A6/ A4 ( a^/a4  X : 


Speed  of  sound  increase  across  the  reflected 
shock  on  cone  surface. 


M6N  (M.  )  : 

D„ 


Total  surface  Mach  number  component  in  a 
plane  normal  to  the  reflection  curve  behind 
the  reflected  shock. 
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M6T(M,  ): 

6t. 


P6/P-INF (p 


SIGMA- 6  (<7 

D 


S6(S6  )i 


Surface  Mach  number  component  tangent  to 
the  reflection  curve,  behind  the  reflec¬ 
ted  shock . 

/p^) ^ : Pres sure  ratio  on  cone  surface  behind  re¬ 
flected  shock. 

) :  Surface  flow  deviation  behind  the  reflected 

shock,  with  respect  to  the  fl  axes. 

Total  surface  entropy  change  behind  the  re¬ 
flected  shock,  normalized  by  twice  the 
specific  heat  at  constant  pressure. 
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APPENDIX  VI 


(SAMPLE  PROBLEMS  -  THREE-DIMENSIONAL) 


1.  Wedge  Program  Example 

To  illustrate  the  use  of  the  wedge  program  a  particular 
example  was  carried  out  for  the  following  input  data: 

X  =  55  deg. 

M  =  4.23 
00 

X  =  8  deg. 

*  =  0 

Az=  0.05 

All  lengths  have  been  normalized  with  respect  to  r  by  in¬ 
putting  r=l.  The  program  output  is  included  in  this  appendix 
and  the  symbols  cure  defined  in  Appendix  V-l. 

Each  set  of  non-dimensionalized  coordinates  x  and  z 
represent  a  point  on  the  intersection  curve  of  the  shock  on 
the  plane  surface.  However,  if  the  ratio  x/z  is  considered 
the  variation  of  flow  conditions  on  the  wedge  surface  is 
determined.  Each  set  of  conditions,  i,  corresponds  to  the 
flow  on  a  ray  on  the  base  plane  of  the  wedge  passing  through 
the  wedge  vertex  and  making  a  angle  tan"*- (x/z)  ^  with  the  z 
axes,.  Due  to  the  assumption  of  linear  theory  the  flow  on  the 
wedge  surface  is  the  same  as  that  on  the  base  plane  of  the 
wedge.  The  last  set  of  conditions  (i  =  22)  corresponds  to  those 
on  the  edge  of  the  local  Mach  cone  and  therefore  also  to  the 
constant  conditions  in  the  two-dimensional  region.  Note  that 
the  step  size  has  been  reduced  by  a  factor  of  5  for  the  last 
set  of  points  in  order  to  approach  the  two-dimensional  values 
more  closely.  Figure  43  presencs  the  variation  of  Mach 
number,  pressure  ratio  and  flow  deviation  across  the  wedge 
surface.  Only  half  the  wedge  is  shown  due  to  its  symmetry. 

These  conditions  are  also  valid  in  front  of  the  reflected 
shock.  The  strength  of  the  reflected  shock  at  any  point  i  on 
the  intersection  curve  depends  on  the  upstream  Mach  number  in  a 
plane  normal  to  the  intersection  curve  at  that  point  and  on  the 
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FIGURE  43  -  WEDGE  SURFACE  CONDITIONS 


deflection  of  the  flow  necessary  to  turn  the  flow  parallel 
to  the  plane  surface.  These  quantities  are  plotted  for  each 
point  on  the  intersection  curve  i..  Figure  43.  At  first, 
one  would  expect  that  the  shock  strength  would  decrease  as 
points  further  and  further  from  the  z  axes  are  considered 
since  the  normal  Mach  number  decreases.  However,  it  is  also 
necessary  to  consider  the  fact  that  the  flow  deflection  in¬ 
creases  as  the  two-dimensional  region  is  approached.  There¬ 
fore,  the  shock  strength  may  increase  or  decrease  depending 
on  the  relative  effects  of  these  two  quantities.  It  is  also 
possible  that  as  the  two-dimensional  region  is  approached  the 
required  deflection  will  become  too  large  for  the  corresponding 
normal  Mach  number,  resulting  in  a  Mach  reflection.  This  did 
not  occur  in  this  example  but  would  result  in  a  program  stop  as 
described  in  Appendix  V. 

The  conditions  immediately  behind  the  intersection  curve 
on  the  plane  surface  are  plotted  in  Figure  45  .  Figures  44 
and  45  do  not  represent  a  cross  section,  but  points  on  the 
intersection  curve  (Figure  3u  of  the  main  text) .  It  can  be 
seen  that  the  shock  strength  decreases  a  little  and  then  in¬ 
creases  as  the  intersection  curve  is  traversed  from  the  z  axes 
to  the  two-dimensional  region.  It  should  be  noted  that  the 
values  of  02n£  can  used  to  construct  the  reflected  shock. 

In  particular  A  is  the  angle  between  the  reflected  shock 

i=l 

and  the  z  axes  in  the  y-z  plane. 
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FIGURE  44  -  CONDITIONS  AFFECTING  REFLECTED  SHOCK  STRENGTH 


FIGURE  45  -  CONDITIONS  BEHIND  INTERSECTION  CURVE 


2.  Conical  Shock-Plane  surface  Program  Example 

This  program  was  run  together  with  the  characteristic 
subroutine  program  for  the  following  input  data. 


M  =4.0 
00 

9^  =  10  deg. 

y  =  0.30918  rad.  (10  deg.  cone) 
c 

y  =1.4 
r  =  2.5 
q  =2.0 
Az  =  0.01 
p  /p  =  1.562 

C  ao 

M*  =  2.09123 
lc 

j h  =  0.0803353  rad. 

rlc 

S2/R=  0.089929 

This  program  output  is  included  in  the  appendix  and  the 
symbols  are  defined  in  Appendix  V.2. 

Corresponding  to  each  value  of  i,  the  set  of  coordinates 
x'  and  z'  represent  a  point  on  the  reflection  curve  of  the 
conical  shock  from  the  plane  surface.  This  curve  also  represents 
the  data  line  for  the  characteristic  subroutine  and  is  shown  in 
Figure  46.  .  The  curve  is  terminated  after  computing  a  maximum 
number  of  25  points.  If  it  is  necessary  to  find  the  flow 
further  out  on  the  curve  the  step  size,  Az,  would  have  to  be  in¬ 
creased.  If  a  large  enough  step  size  were  used  the  reflection 
curve  computation  would  terminate  either  at  a  Mach  reflection  or 
when  the  slope  of  the  curve  with  respect  to  the  z'  axes  becomes 
5°  or  less. 


The  variation  of  the  reflected  shock  strength  can  be 

determined  from  Figure  47  where  the  total  upstream  Mach  number 

component  in  a  plane  normal  to  the  reflection  curve  (M  )  and 

2N. 

i 

and  the  flow  deflection  necessary  to  turn  the  flow  parallel  to 
the  plane  surface  (62  )  are  plotted  vs.  the  slope 

i 
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FIGURE  46  -  INTERSECTION  CURVE  OF  CONICAL  SHOCK  ON  PLANE  SURFACE 


FIGURE 


47  - 


CONDITIONS  AFFECTING  STRENGTH  OF 
REFLECTED  CONICAL  SHOCK 


of  the  reflection  curve  .  The  effect  of  these  parameters 
is  described  in  Appendix  VI . 1  where  it  is  also  pointed  out 
that  this  figure  does  not  represent  a  cross-section  but  points 
along  the  reflection  curve. 

The  conditions  on  the  plane  surface  behind  the  reflection 
are  given  in  Figure  48.  The  first  part  of  the  figure  shows 
the  pressure  and  Mach  number  variations  immediately  behind  the 
reflection  curve.  It  can  be  seen  that  a  Mach  relection  must 
occur  quite  soon  since  the  pressure  is  increasing  so  rapidly. 
The  second  part  of  the  figure  shows  the  variation  of  the  pres¬ 
sure  along  the  line  of  symmetry  (p  or  z'  axes)  on  the  plane 
surface  as  obtained  from  the  characteristic  subroutine  program 

Something  should  be  said  here  about  the  output  of  the 
characteristic  subroutine.  The  first  page  gives  the  value  of 
the  pressure  behind  the  reflection  curve  on  the  centerline  and 
the  flow  conditions  along  the  reflection  curve,  which  is  now 
called  the  data  line.  The  data  is  given  in  reverse  order,  i.e 
the  point  furthest  out  on  the  curve  is  listed  first  and  the 
centerline  point  is  given  last.  The  next  page  gives  the  condi¬ 
tions  along  the  second  family  characteristic  line  closest  to 
the  centerline  data  point.  This  line  is  designated  by  the 
value  of  i,  and  the  value  of  p/pw  given  at  the  top  is  the 
pressure  ratio  at  the  intersection  of  this  second  family 
char acteristic  line  and  the  centerline  (wall) .  The  procedure 
continues  until  the  second  family  line  from  the  data  point 
furtherst  out  on  the  reflection  curve  has  been  considered. 
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FIGURE  48  -  CONDITIONS  BEHIND  REFLECTED  SHOCK  ON  PLANE  SURFACE 
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3.  Double  Conical  Shock  Program  Example 


This  nronran!  was  carried  out  for  the  case  where  one 
cone  is  at  an  angle  of  attack  for  the  following  input  data* 

M*  =  2.03183  (M  =4.0,  15  deg.  cone). 

ij)*  =  0.165935  rad. 

y^  =  0.38032  rad. 

p/p  =  2.4056 
1  ® 

Sj/R  =  0.067733 

M*  =  2.06346 

X  t  z 

if)*  =  0.12262  rad. 

/  z 

y  =  0.342957  rad.  (M  =  4.0,  12.5  deg.  cr  , 
p  /p  =  1.9441 

2  CO 

S  /R  =  0.029667 
*  t  z 

r  =  2.0 

q  =  1.0 

A>[;  =  5.0  deg. 

y  =  1.4 

S2/R  =  0.52932 

KALF  =  1 

O'  =  2.0  deg. 

Tj/a  =  0.654926 

M2*  =  -0.408226 

M*  ,,  —  1.86308 

w,2 

0*  =-0.505885  rad. 
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The  program  output  is  included  in  this  appendix  and  the 
symbols  are  defined  in  Appendix  V.3. 

The  intersection  curve  of  the  two  conical  shocks  is 
described  in  terms  of  conical  coordinates  p  and  ip,  which 
lie  on  the  base  cone  surface.  The  computation  of  this 
curve  is  terminated  when  the  value  of  ipmax  >  which  repre¬ 
sents  the  point  where  the  intersection  curve  starts  to 
bend  back  in,  is  reached. 

The  pressure,  Mach  number  and  entropy  immediately 
behind  the  intersection  curve  in  the  lower  and  upper  regions 
are  given  by  Figures  49  and  50.  These  quantities  are 
plotted  versus  the  meridian  angle  tp  and  therefore  are  at 
points  along  the  intersection  curve.  As  shown  in  the  out¬ 
put  these  computations  are  terminated  when  a  so-called  Mach 
reflection  is  reached.  This  means  that  the  normal  Mach 
number  components  are  too  small  for  the  necessary  deflections 
and  a  cusped  intersection  with  a  normal  shock  region  occurs. 
Comparison  of  Figures  49  and  50,  shows  that  the  pressure 
ratios  in  the  lower  and  upper  regions,  at  each  point  on  the 
intersection  curve,  differ  by  less  than  2%  even  though  the 
assumption  of  constant  shock  strengths  has  been  used  instead 
of  an  iteration  procedure. 
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FIGURE  49  -  LOWER  CONICAL-CONICAL  INTERSECTION  REGION  PROPERTIES 


UPPER  CONICAL-CONICAL  INTERSECTION  REGION  PROPERTIES 
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TABLE  VI  (Cont'd) 


4. 


The  intersection  curve  of  the  plane  and  conical  shocks 
in  coordinates  on  the  plane  surface  is  given  in  Figure 
51.  As  noted  on  the  first  page  of  the  output,  this  com¬ 
putation  terminates  when  a  so-called  Mach  reflection  is 
reached,  which  means  a  cusped-normal  shock  intersection 
exists  in  the  plane  normal  to  the  intersection  curve. 
Figure  52  shows  the  variation  of  Mach  number,  pressure 
ratio  and  entropy  in  the  lower  intersection  region,  immedi¬ 
ately  behind  the  intersection  curve  (conditions  4)  versus 
the  slope  of  the  curve  with  respect  to  the  z'  axes  (Q.). 


The  program  next  computes  the  intersection  curve  of  the 
modified  plane  shock  and  the  cone  surface.  To  do  this  it 
is  necessary  to  pick  three  points  through  which  to  construct 
the  second  order  curve  (Section  E  ) .  The  points  chosen  in 
this  case  are  i  =  1,4,  and  10.  This  curve  is  terminated  due 
to  a  Mach  reflection.  The  flow  properties  M  ,  p  /p 

bi  4bi  °° 

and  immediately  behind  the  modified  plane  shock  near 

"U  ^ 

the  cone  surface  are  shown  in  Figure  53  and  differ  some¬ 
what  from  those  in  Figure  52  since  the  flow  upstream  of 
the  modified  plane  shock  is  conical.  The  flow  properties 
immediately  behind  the  reflected  shock  on  the  cone  surface 
are  plotted  versus  the  slope  of  the  reflection  curve  on  the 
conical  surface,  /3,  given  by 


+  tan 


-1 


in  Figure  54.  The  pressure  variation  along  the  centerline 
on  the  cone  surface  obtained  from  the  characteristic  subroutine 
is  given  in  Figure  55.  The  output  of  the  characteristic  sub¬ 
routine  is  described  in  Appendix  VI-2. 


FIGURE  52  -  CONDITIONS  IN  LOWER  INTERSECTION  REGION  OF  PLANE  AND  CONICAL  SHOCKS 


FIGURE  53.  -  FLOW  PROPERTIES  BEHIND  THE  MODIFIED  PLANE  SHOCK 

NEAR  THE  CONE  SURFACE 
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